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Relativistic interaction of short-pulse lasers with underdense plasmas has recently led 
to the emergence of a novel generation of femtosecond x-ray sources. Based on radi- 
ation from electrons accelerated in plasma, these sources have the common properties 
to be compact and to deliver collimated, incoherent and femtosecond radiation. In this 
article we review, within a unified formalism, the betatron radiation of trapped and 
accelerated electrons in the so-called bubble regime, the synchrotron radiation of laser- 
accelerated electrons in usual meter-scale undulators, the nonlinear Thomson scattering 
from relativistic electrons oscillating in an intense laser field, and the Thomson backscat- 
tered radiation of a laser beam by laser-accelerated electrons. The underlying physics 
is presented using ideal models, the relevant parameters are defined, and analytical ex- 
pressions providing the features of the sources are given. Numerical simulations and a 
summary of recent experimental results on the different mechanisms are also presented. 
Each section ends with the foreseen development of each scheme. Finally, one of the 
most promising applications of laser-plasma accelerators is discussed: the realization of 
a compact free-electron laser in the x-ray range of the spectrum. In the conclusion, the 
relevant parameters characterizing each sources are summarized. Considering typical 
laser-plasma interaction parameters obtained with currently available lasers, examples 
of the source features are given. The sources are then compared to each other in order 
to define their field of applications. 
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I. INTRODUCTION 

X-ray radiation has been, ever since its discovery over 
a century ago, one of the most effective tools to explore 
the properties of matter for a broad range of scientific re- 
search. Successive generations of radiation sources have 
been developed providing radiation with always higher 
brightness, shorter wavelength and shorter pulse dura- 
tion (Koch, 1983). Despite remarkable progress on x-ray 
generation methods, there is still a need for light sources 
delivering femtosecond pulses of bright high-energy x-ray 
and gamma-ray radiation, emitted from source size of 
the order of a micron (Pfeifcr et ai, 2006; Service, 2002). 
Indeed, the intense activity on the production of such 
radiation is motivated by countless applications in fun- 
damental science, industry or medicine. (Bloembergen, 
1999; Martin et at, 1992; Rousse et ai, 2001b; Zewail, 
1997). For example, in the studies of structural dynam- 
ics of matter, the ultimate time scale of the vibrational 
period of atoms is a few tens of femtoseconds. Fun- 
damental processes such as dissociation, isomerization, 
phonons, and charge transfer evolve at this time scale. 
High-energy radiation is used to radiograph dense ob- 
jects that are opaque for low-energy x rays, while micron 
source size allows one to obtain high-resolution images 
and makes possible phase contrast imaging to see what is 
invisible with absorption radiography. Several techniques 
are being developed to produce femtosecond x rays. In 
the accelerator community, large-scale free-electron laser 
facilities can now deliver the brightest x-ray beams ever, 
with unprecedented novel possibilities (Barty et al., 2008; 
Brock, 2007; Chapman et ai, 2006; Fritz et al, 2007; 
Gaffney and Chapman, 2007; Marchcsini et ai, 2008; 
Neutze et ai, 2000). The slicing technique, combining 
a conventional accelerator with a femtosecond laser to 
isolate short electron slices, allows synchrotrons to pro- 
duce radiation pulses with duration of the order of 100 fs 
(Schoenlein et ai, 2000). High-energy radiation can be 
delivered by radioactive sources, x-ray tube, and Comp- 
ton scattering sources based on a conventional accelera- 
tor. However, even if widely used these high-energy radi- 
ation sources have limitations in terms of storage, pulse 
duration, spectrum tunability, energy range, and source 
size. 



In parallel, alternative and complementary methods 
based on laser-produced plasmas have been developed to 
produce ultrashort compact radiation sources covering a 
wide spectral range from the extreme ultraviolet (XUV) 
to the gamma rays. While several laser-based source 
schemes were proposed in the early 1970s, this field of 
research has seen rapid development when lasers have 
been able to produce intense femtosecond pulses (Perry 
and Mourou, 1994; Strickland and Mourou, 1985). At 
laser intensities on the order of 10 14 W/cm 2 , XUV radi- 
ation, in the few tens of electronvolts (eV) energy range, 
can be produced using the mechanism of high-order har- 
monics generation from gas targets (Brabec and Krausz, 
2000; Corkum, 1993; Krausz and Ivanov, 2009; Protopa- 
pas et ai, 1997) or by XUV laser amplification in a laser- 
produced plasma (Daido, 2002). These sources can de- 
liver, in most recent configurations, up to a microjoule 
of radiation within a beam of a few milliradians diver- 
gence. At laser intensities on the order of 10 16 W/cm 2 , 
x-ray sources from laser solid target interaction can pro- 
duce a short pulse of K a line emission, emitted within 
4-7T steradians (Kicffer et al., 1993; Murnane et al., 1991; 
Rousse et al, 1994). Discovered more than a decade ago, 
these sources have been widely developed and have led 
to the first structural dynamics experiments at the fem- 
tosecond time scale (Cavalleri et al, 2001; Rischel et ai, 
1997; Rose-Petruck et ai, 1999; Rousse et ai, 2001a; 
Siders et ai, 1999; Sokolowski-Tinten et ai, 2003, 2001). 
With recent developments, laser systems can deliver fo- 
cused intensities above 10 18 W/cm 2 and the laser-plasma 
interaction has entered the relativistic regime (Mourou 
et ai, 2006; Umstadter, 2003). At this light intensity, 
relativistic effects become significant. Electrons can be 
accelerated within the laser field or in the wakefield of 
the laser up to relativistic energies (Esarey et ai, 2009, 
1996b; Everett et ai, 1994; Joshi et ai, 1984; Malka et ai, 
2002; Modena et ai, 1995; Patel, 2007; Tajima and Daw- 
son, 1979; Umstadter et ai, 1996a). In particular, laser 
wakefield acceleration has led to the production of high- 
quality femtosecond relativistic electron bunches (Faure 
et ai, 2004, 2006; Geddes et al, 2004; Mangles et ai, 
2004) created and accelerated up to the gigaelectron- 
volt level (Hafz et ai, 2008; Kneip et ai, 2009; Lecmans 
et ai, 2006) within only a few millimeters or centimeters 
plasma. Using these relativistic electrons, several novel 
x-ray source schemes have been proposed over the past 
decades to produce collimated and femtosecond radiation 
in a spectrum ranging from the soft x rays to gamma 
rays. Most of these schemes are based on the wiggling 
of relativistic electrons accelerated in a laser wakefield. 
In this article, the physics of these sources is reviewed, 
and the opportunities offered by these relativistic elec- 
trons to generate ultrashort x-ray radiation (Catravas 
et ai, 2001; Fritzlcr et ai, 2003; Griincr et al, 2007; 
Hartemann et ai, 2007; Jaroszynski et ai, 2006; Lee- 
mans et ai, 2005; Malka et ai, 2008; Nakajima, 2008) 
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are highlighted. These sources can deliver x rays or 
gamma rays as short as a few femtoseconds, as they in- 
herit the temporal profile of the laser-plasma electron 
bunch, whose few-femtosecond duration was recently ex- 
perimentally demonstrated (Lundh et ai, 2011). 

The aim of this article is to review the novel x-ray 
sources based on relativistic laser and underdense plasma 
interaction and to highlight their similitude by using a 
common formalism for their description. The paper is 
organized as follows. In Sec. II, the general formalism of 
radiation from an accelerated relativistic electron is pre- 
sented, which provides a framework for the description 
of the sources discussed throughout the paper. From 
this formalism, the relevant parameters describing the 
properties of the radiation, such as its spectrum, diver- 
gence, number of emitted photons and duration, can be 
extracted. As the x-ray sources presented here are based 
on electrons accelerated by laser wakefields [i.e., by the 
laser wakefield accelerator (LWFA)], a description of the 
most efficient laser-based electron accelerator to date is 
given in Sec. III. 

In Sees. IV, V and VI, different methods for the pro- 
duction of incoherent x rays from relativistic electrons are 
reviewed; the objective is to define the relevant regimes 
to accelerate and wiggle electrons in such a way that they 
emit x rays. In Sec. IV, betatron radiation is described. 
In that plasma cavity created in the wake of an 

intense laser pulse acts as both an electron accelerator 
and a wiggler (referred to as a plasma undulator). The 
scheme presented in Sec. V relies on the use of laser wake- 
field accelerated electrons, transported and wiggled in a 
meter-scale periodic arrangement of permanent magnets 
(conventional undulator). This method is the closest to 
synchrotron technology. In Sec. VI, the nonlinear Thom- 
son scattering and the Thomson backscattering sources 
are reviewed. In nonlinear Thomson scattering, electrons 
are directly accelerated and wiggled in an intense laser 
field. For the Thomson backscattering case, the plasma 
is used to accelerate electrons which are then wiggled in 
a counterpropagating electromagnetic wave (EM undula- 
tor). 

In Sec. VII, an introduction to the topic of the free- 
electron laser is given and the conditions for realizing 
such an ultrahigh-brightness and coherent source of ra- 
diation with wavelength down to the angstrom (hard x 
rays) from laser-accelerated electrons are discussed for 
different types of undulators (plasma, conventional, and 
EM). 

Other radiation sources based on the laser-plasma in- 
teraction have been developed and could provide pho- 
tons in the keV range. For example, high-order harmon- 
ics from gas (Brabec and Krausz, 2000; Corkum, 1993; 
Krausz and Ivanov, 2009; Protopapas et al., 1997) or solid 
targets (Dromcy et ai, 2006; Tarasevitch et at, 2000; 
Teubncr and Gibbon, 2009; Thaury et ai, 2007), sources 



based on the flying mirror concept (Bulanov et ai, 2003; 
Esirkepov et ai, 2009; Kando et at, 2007; Mourou et ai, 
2006), or the K a source (Kicffcr et ai, 1993; Murnane 
et ai, 1991; Rousse et ai, 1994) can produce femtosec- 
ond XUV or x-ray radiation. However, they are not based 
on the same physical principle of acceleration and wig- 
gling of relativistic electrons, and will therefore not be 
reviewed here. 



II. GENERAL FORMALISM: RADIATION FROM 
RELATIVISTIC ELECTRONS 

In this section, the radiation from relativistic electrons 
is introduced. A qualitative understanding of the phe- 
nomenon is highlighted and analytical results for the rel- 
evant general parameters determining the radiation fea- 
tures are given. The formalism described here is general 
and will be common to all the sources presented through- 
out the article. We follow the approach of Jackson (2001) 
and provide the basic results necessary for the remainder 
of the paper. The interested reader is referred to Wiede- 
mann (2007a) for a complete and detailed description of 
synchrotron radiation from bending magnets, undulator, 
and wiggler insertion devices (e.g. for the angular distri- 
bution of individual undulator harmonics, for polariza- 
tion or spatial and temporal coherence properties). 



A. Radiation features 

Relativistic electrons can produce bright x-ray beams 
if their motion is appropriately driven. For all the laser- 
based x-ray sources discussed in this article, the radia- 
tion mechanism is the emission from accelerated relativis- 
tic electrons. The features of this relativistically moving 
charge radiation are directly linked to the electron tra- 
jectories. Obtained from the Lienard-Wiechert field, the 
general expression that gives the radiation emitted by 
an electron, in the direction of observation ft, as a func- 
tion of its position, velocity, and acceleration along the 
trajectory is written (Jackson, 2001) 

d 2 I e 2 



du>d£l 167r 3 eoc 



+oo n x 

iuj[t-n.r(t)/c] 



(i-p.ny 



dt 



(1) 



This equation represents the energy radiated within a 
spectral band duj centered on the frequency u> and a solid 
angle dQ centered on the direction of observation ft. Here 
r(t) is the electron position at time t, f3 is the velocity 
of the electron normalized to the speed of light c, and 

P = d(3/dt is the usual acceleration divided by c. We 
stress that this expression assumes an observer placed at 
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a distance far from the electron so that the unit vector 
n is constant along the trajectory. The expression (1) 
for the radiated energy shows an important number of 
generic features: 

1. When f3 — 0, no radiation is emitted by the elec- 
tron. This means that the acceleration is responsi- 
ble for the emission of electromagnetic waves from 
charged particles. 

2. According to the term (1 — (3.ri)~ 2 , the radiated en- 
ergy is maximum when f3 . n — > 1. This condition 
is satisfied when (3 ~ 1 and f3 \\ n. Thus, a rela- 
tivistic electron (/3 ~ 1) will radiate orders of mag- 
nitude higher than a nonrelativistic electron, and 
its radiation will be directed along the direction of 
its velocity. This is simply the consequence of the 
Lorentz transformation: for an electron emitting 
an isotropic radiation in its rest frame, the Lorentz 
transformation implies that the radiation is highly 
collimated in the small cone of typical opening an- 
gle of A8 = 1/7 around the electron velocity vector, 
when observed in the laboratory frame (see Fig. 1). 
In the following, we consider ultrarelativistic elec- 
trons, 7> 1, and all angles which will be defined 
are supposed to be small so that tan# ~ sin 8 ~ 9. 

3. The term (ft — (3) x (3, together with the relation 

/3ii oc -Fy/7 3 and /3± oc F±/"f between applied force 
and acceleration (respectively for a force longitudi- 
nal or transverse with respect to the velocity 
indicate that applying a transverse force F _L f3 is 
more efficient than a longitudinal force. The term 
also shows that the radiated energy increases with 

the square of the acceleration (3. More precisely, 
P oc Fu and P oc 7 2 -Fj, where P is the radiated 
power. Thus, it is much more efficient to use a 
transverse force in order to obtain high radiated 
energy. 

4. The phase term e M*- n r ( t )/ c ] can be locally ap- 
proximated by e 1 '^ 1- ' 3 )*. The integration over time 
will give a nonzero result only when the integrand, 
excluding the exponential, varies approximately at 
the same frequency as the phase term which oscil- 
lates at oj v ~ oj(1 — (3). Given that the velocity 
P of the electron varies at the frequency w e -, the 
condition lo v ~ cj e - is required to have a nonzero 
result. The electron will radiate at the higher fre- 
quency u; = ui e - /{I — j3) ~ 27 2 w e -. Thus, the usual 
Dopplcr upshift is directly extracted from this gen- 
eral formula. This indicates the possibility to pro- 
duce x-ray beams {uj x ~ 10 18 s" 1 ) by wiggling a 
relativistic directional electron beam at a frequency 
far below the x-ray range: w e - ~ ujx / i^l 2 )- 



This analysis underlines the directions for the produc- 
tion of x rays from relativistic electrons: the goal for 
x-ray generation from relativistic electrons is to force a 
relativistic electron beam to oscillate transversally. This 
transverse motion will be responsible for the radiation. 
This is the principle of synchrotron facilities, where a pe- 
riodic static magnetic field, created by a succession of 
magnets, is used to induce a transverse motion to the 
electrons. The laser-based sources presented here rely 
on this principle. In the next sections, the properties 
of moving charge radiation in two different regimes are 
reviewed. The different laser-based x-ray sources which 
will be discussed can work in both regimes depending on 
the interaction parameters. 



B. Two regimes of radiation: undulator and wiggler 

We consider ultrarelativistic electrons with a velocity 
along the direction e z executing transverse oscillations in 
the e x direction. Two regimes can be distinguished. 

The undulator regime corresponds to the situation 
where an electron radiates in the same direction at all 
times along its motion, as shown in Fig. 1. This occurs 
when the maximal angle of the trajectory ip is smaller 
than the opening angle of the radiation cone A9 = 1/7. 
The wiggler regime differs from the undulator by the 
fact that the different sections of the trajectory radiate 
in different directions. Thus, emissions from the differ- 
ent sections are spatially decoupled. This occurs when 
if) I/7. The fundamental dimensionless parameter 
separating these two regimes is K = "ftp. The radiation 
produced in these two regimes have different qualitative 
and quantitative properties in terms of spectrum, diver- 
gence, and radiated energy and number of emitted pho- 
tons. 



C. Qualitative analysis of the radiation spectrum 

The shape of the radiation spectrum can be determined 
using simple qualitative arguments. In most of the cases 
discussed throughout the article the electron trajectory 
can be approximated by a simple transverse sinusoidal 
oscillation of period X u at a constant velocity j3 and con- 
stant 7. The orbit can be written as 

tp K 
x(z) = x sin(k u z) = — sin(fc„z) = — — sin(fc M z), (2) 
ku 7"'« 

where k u — 2tt/X u is the wave- vector norm, xo is the 
transverse amplitude of motion, and tj; is the maximum 
angle between the electron velocity and the longitudinal 
direction e z . Since the electron energy is constant, an 
increase of the transverse velocity leads to a decrease of 
the longitudinal velocity. This can be explicitly derived 




Wiggler 
K»l 



e" orbit 



FIG. 1 Illustration of the undulator and wiggler limits, at 
the top and the bottom, respectively. The lobes represent the 
direction of the instantaneously emitted radiation, ip is the 
maximum angle between the electron velocity and the prop- 
agation axis e z and AO is the opening angle of the radiation 
cone. When ip <C AO (undulator), the electron always radi- 
ates in the same direction along the trajectory, whereas when 
ip 3> AO (wiggler), the electron radiates toward different di- 
rections in each portion of the trajectory. 



from the assumed trajectory (2), 



P, - P ( 1 

T*-P{i- 



-) 



COS 2 (fc u 2:) 
~ 1 



1 , K\ 



(3) 
(4) 



With the trajectory of the electron periodic, the emit- 
ted radiation is also periodic since each time the electron 
is in the same acceleration state, the radiated amplitude 
is identical. The period of the radiated field can be cal- 
culated to obtain the fundamental frequency of the radi- 
ation spectrum. The radiation emitted in the direction 
n, forming an angle 8 with the e z direction, is consid- 
ered, as represented in Fig. 2. The field amplitude A\ 
radiated in the direction ft by the electron at z = and 
t = propagates at the speed of light c. At z = X u 
and t = \ u /(f3 z c), the electron radiates an amplitude 
A 2 = Ax. The distance separating both amplitudes {A\ 
and A2) corresponds to the spatial period A of the radi- 
ated field and is given by 



A 



A,, 

K 



A„ cos ( 



An 
2f 



(i + ^r+7 2 0. (5) 



The radiation spectrum consists necessarily in the funda- 
mental frequency uj = 2t:c/X and its harmonics. To know 
if harmonics of the fundamental are effectively present in 
the spectrum, it is instructive to look at the electron mo- 
tion in the average rest frame, moving at the velocity /3 Z 
in the e z direction with respect to the laboratory frame. 




FIG. 2 Schematic for the calculation of the spatial period A of 
the radiation emitted toward the direction of observation n, 
forming an angle with the e z direction. At the two positions 
marked by a blue point, the electron radiates the same field 
amplitude. The distance between these two amplitudes at a 
given time corresponds to A. 



If K <C 1, the longitudinal velocity reduction is negli- 
gible: /3 Z — P an d 7z = 1/ yl — Pz — 7- The motion 
contains only the fundamental component. Indeed, the 
motion reduces to a simple dipole in the average rest 
frame. The spectrum consists in a single peak at the 
fundamental frequency uj which depends on the angle of 
observation 0. As K — > 1, radiation also appears at har- 
monics. 

If K 1, the longitudinal velocity reduction is sig- 
nificant: 7 2 = 7/ J I + \K 2 . In the average rest frame, 
the trajectory is a figure-eight motion. It can contain 
many harmonics of the fundamental. In the laboratory 
frame, this can be explained by the fact that an ob- 
server receives short bursts of light of duration r. In- 
deed, the instantaneous radiation is contained within a 
cone of opening angle A6 = 1/7 centered on p and points 
toward an observer positioned in the direction ft during 
a time At (see Fig. 3), corresponding to the variation 
of P by an angle Ad = 1/7. Locally, a portion of the 
trajectory can be approximated by a portion of a cir- 
cle of radius p, such that the direction of the velocity 
P changes by an angle AO when the electron travels a 
distance d e = 27rp x (A0/2ir) = p/j, which corresponds 
to a time At = t e — d e /(pc). During the time At, the 
radiation has covered a distance <i 7 = 2psin(l/27) cor- 
responding to a propagation time of t~ f = dj/c. The 
radiation burst duration r as observed by an observer 
reads 

13p 
24 7 V 

The temporal profile (see the inset of Fig. 3) of the radia- 
tion emitted in the wiggler regime has been qualitatively 
obtained. The Fourier transform of this typical profile 
gives a precise representation of the radiation spectrum. 
With the time profile a succession of bursts of duration 
r, the spectrum will contain harmonics up to the critical 
frequency 



(6) 



1/r ~ 7 3 



(7) 
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quency: 



FIG. 3 In the wiggler limit, the radiation cone points toward 
the observer during a time At, which corresponds to a du- 
ration r for the emitted radiation. This is repeated at each 
period: the observer receives bursts of radiation separated 
by a time A/c. The inset gives the temporal profile of the 
radiation power seen by the observer. 



Note that the spectrum that arises from a complete cal- 
culation of the radiation emitted by a relativistic charged 
particle in instantaneously circular motion (Jackson, 
2001) is in agreement with the above estimation. Such 
calculation yields the synchrotron spectrum, which is 
written in terms of radiated energy per unit frequency 
and per unit time, 



dP _ 
duj 

S(x) = 

P-v = 



p 

— S(oj/uj c ), 

2 4 

e q 



K 5/3 (0d£, 

67reop 2 277reoC7 2 ' 
3 o c 
2 p 



(8) 



where we introduced the exact definition of the critical 
frequency lo c of the synchrotron spectrum which will be 
used throughout the review, P 7 = J dP/duj is the radi- 
ated power, and K 5 / 3 is the modified Bessel function of 
the second kind. 

The expression for the radius of curvature p can be 
obtained for an arbitrary trajectory. Its value can be 
calculated at each point of the trajectory, corresponding 
to a particular direction of observation. In the particular 
case of a sinusoidal trajectory, the radius of curvature 
reads p(z) = po[l + ^ 2 cos 2 (fc„z)] 3 / 2 /| sin(k u z)\ and it 
is minimum when the transverse position is extremum 
(x = ±xq) and its value at this point reads 



Po 



(9) 



lu c = -K-f 2itc/\ u 



-Ku 



{K«i,e=o>- 



(10) 
(11) 



Note that for a nonplanar and nonsinusoidal trajectory, 
the spectrum extends up to a critical frequency deter- 
mined by the minimal radius of curvature of the trajec- 
tory. 

The parameter K can be considered as the number of 
decoupled sections of the trajectory. With each section 
radiating toward a different direction, the radiation is 
spatially decoupled and leads to bursts of duration r in 
each direction and to a broad spectrum with harmonics 
of the fundamental up to tu c . 

An infinite periodical motion has been considered so 
far, leading to harmonics that are spectrally infinitely 
thin. For a finite number of oscillation periods N, Fourier 
transform properties imply that the harmonic of number 
n and of wavelength A„ = X/n has a width given by 



AA» 
X n 



1 

nN' 



(12) 



This corresponds to the harmonic bandwidth in a given 
direction characterized by the angle 9. Since the wave- 
length A„ depends on the angle 8 [see Eq. (5)], the inte- 
gration of the radiation over a small aperture of finite di- 
mension broadens each harmonic. When integrating over 
the total angular distribution, harmonics overlap and the 
radiation spectrum becomes continuous, but keeps the 
same extension (up to lj c ). In the undulator limit, in 
which only the fundamental wavelength is present, the 
bandwidth is highly degraded when integrating over the 
total angular distribution. 



D. Duration and divergence of the radiation 

The pulse duration and the divergence of the radiation 
emitted by a single electron can be deduced from the 
previous analysis. If N is the number of periods of the 
trajectory, the radiation consists in N periods of length 
A and the total duration of the pulse is T r \ N =1 = NX/c. 

For an undulator, K <C 1, the direction of the vector (3 
varies along the trajectory by an angle ip negligible com- 
pared to A8 = 1/7. Hence, the radiation from all sections 
of the trajectory overlap and the typical opening angle of 
the radiation is simply 8 r = I/7 1 . For a wiggler, K 3> 1, 



leading to the following expression for the critical fre- 



1 For K < 1, the root mean square angle of the angular distribution 
of the radiated energy at the fundamental frequency is (0 2 ) 1 / 2 = 
l/7z, which simplifies to I/7 for K 1 (Jackson, 2001). 
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the divergence increases in the direction of the transverse 
oscillation e x , but remains identical in the orthogonal di- 
rection e y . In the direction of the motion e x , the typical 
opening angle of the radiation is 9x r — tp — K/j 2 , which 
is greater than A9 — I/7, while the typical opening an- 
gle of the radiation in the direction e y is By r — 1/7. For 
the general case of a transverse motion occurring in the 
e x and e y directions, the angular profile can take various 
shapes depending on the exact three-dimensional (3D) 
trajectory. 



E. Analytical formulas for the total radiated energy and the 
number of emitted photons 

Analytical calculations provide simple expressions for 
the radiated energy and the number of emitted photons 
per period. Using the expression of the radiated power by 
an electron P(t) = (e 2 /6Tre Q c)j 2 [(dj)/dt) 2 -(d-f/dt) 2 } (fis 
the momentum normalized to mc), the averaged radiated 
power P 7 and the total radiated energy per period J 7 can 
be derived for both the undulator and the wiggler case for 
an arbitrary trajectory. In the case of a planar sinusoidal 
trajectory, the result is 

7re 2 c ~/ 2 K 2 
3e K 
Tre 2 ~f 2 K 2 
3e 



P 7 = 



7 T = 



A„ 



(13) 



(14) 



To obtain an estimation of the number of emitted pho- 
tons _/V 7 , a mean energy of photons must be determined. 
For K <C 1, the spectrum is quasimonochromatic in the 
forward direction and the mean energy of photons, af- 
ter integrating over the angular distribution, is equal to 
Tkvg—o/2. The number of emitted photons reads 

2tt 



(15) 



where a = e 2 / (At: eq he) is the fine structure constant. 
For K 1, the spectrum is synchrotronlike with the crit- 
ical frequency oj c . Using the fact that for a synchrotron 
spectrum (fku) = (8/15y/3)hw c , the following estimation 
is derived: 



iV 7 = 



5\/3tt 



aK. 



(16) 



F. Radiation from an ideal electron bunch 

The radiation from a single electron has been discussed 
so far. We now take into account the fact that there arc 



2 Here the typical opening angle is defined as the maximum deflec- 
tion angle of the electron trajectory ij>, such that the full width 
of the angular distribution of the radiated energy in the direction 
e x is 2K/f. 



N e electrons contained in the bunch, assuming that they 
all have exactly the same energy and the same initial 
momentum (zero emittance) . The radiation from several 
electrons is obtained by summing the contribution of each 
electron before taking the squared norm, 



d 2 I 
dujdVL 



e 



N e 

E 



167r 3 eoc 

+00 



iiAj(t-n.fj(t)/c) 



(n- /3 j) x Pj 
(l-4n)* 



dt 



(17) 



This formula expresses the coherent addition of the ra- 
diation field of each electron. It can be considerably 
simplified by considering that all electrons follow simi- 
lar trajectories linked to each other by a spatiotemporal 
translation (tj,Rj) of a reference trajectory r(t): 

r ] {t) = R J +r(t-t J ), 
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dcudn 



N e 



iu>(tj—n.Rj j c) 



e 

167r 3 eoc 



+00 n x 

iu(t-n.r{t)/c) 



(1 - (3.nf 



dt 



(18) 



The radiated energy per unit frequency and unit solid 
angle from an electron bunch is equal to the radiated 
energy from a single electron following the trajectory r(t) 
multiplied by the coherence factor 



c(w) 



gMitj—ri.Rj/c) 

3 = 1 



(19) 



The value of c(w) depends on the electron distribution in 
the (tj,Rj) space. For a uniform distribution, c{uS) = 0, 
whereas for a random distribution, c(uS) — N e on aver- 
age. If the distribution is microbunched at the wave- 
length Ah = 27rc/wb, the summation is coherent for the 
frequency uib and its harmonics, c(nwj) = N 2 for n £ N*. 

In large accelerators or in laser-plasma accelerators, 
electrons are randomly distributed inside the bunch at 
the x-ray wavelength scale, and the radiation is incoher- 
ently summed, 



c{u) = N e , 
dujdil\N<! dudfl I jv.=i 



(20) 
(21) 



The spectrum shape and the radiation divergence 6 r re- 
main unchanged for an electron bunch. The temporal 



profile of the radiation is given by the convolution be- 
tween the electron bunch temporal profile and the ra- 
diation profile from a single electron. Since the typical 
electron bunch length is in the micron range, whereas 
the radiation length from a single electron l r \ N =1 = NX 
is in the nanometer range for x rays, the duration of the 
radiation from an electron bunch is in most cases approx- 
imately equal to the bunch duration, r r \ N = rj. 

In addition, the electron experiences the radiation from 
other electrons in the case of a bunch, which can modify 
its motion and its energy. This interaction of the bunch 
with its own radiation can lead to a microbunching of the 
electron distribution within the bunch at the fundamen- 
tal wavelength of the radiation and its harmonics. This is 
the free-electron laser (FEL) process which produces co- 
herent radiation. Here the coherence factor c(nojb) is N 2 
instead of N e , indicating that the FEL radiates orders of 
magnitude higher than conventional synchrotrons. How- 
ever, the FEL effect requires stringent conditions on the 
electron beam quality and an important number N of 
oscillations. At the present status of laser-plasma accel- 
erators, realizing a FEL represents a technological chal- 
lenge. In the following sections, the interaction between 
the electron bunch and its radiation will not be taken into 
account because in these schemes the conditions required 
for the FEL are not fulfilled: the electron distribution re- 
mains random along the propagation and the radiation 
is incoherent. In Sec. VII, the underlying physics of 
the free-electron laser is presented in more detail and the 
possible realization of such high-brightness coherent ra- 
diation using laser-plasma accelerators is discussed. 



G. Radiation reaction 

The radiation reaction (RR) corresponds to the ef- 
fect of the electromagnetic field scattered by an electron 
on itself, the so-called self-interaction (Hartemann, 2002; 
Jackson, 2001; Landau and Lifshitz, 1994). RR effects 
can modify the electron trajectory and its energy and 
it is therefore important to define the range of param- 
eters for which these effects come into play. The main 
steps, followed by Dirac to derive the relativistically co- 
variant form of the self-force associated with RR, are 
to solve for the self-quadripotential scattered by the 
electron in terms of Green's functions, and then to calcu- 
late the associated electromagnetic force on the electron, 
= -e{d^A s v - d„A^)u v , where u v = dx v /dr is the 
electron quadrivelocity with r the electron proper time. 
This leads to the Dirac-Lorentz equation of motion for a 
pointlikc electron, 



dr 



= ~eFn V u v + t 



dr 2 



( dpv 


dp v \ 


\ dr 


dr J 



where the first term is the Lorentz force with the ex- 
ternal electromagnetic field tensor, the last two terms 
that define the RR self-force correspond, respectively, 
to the Schott term and the radiation damping term, 
= mtip is the electron quadrimomcntum and tq = 
2r e /(3c) = e 2 /(67re mc 3 ) = 6.26 x 10" 24 s, with r e the 
classical electron radius. The Schott term r^d 2 p^ = 
—d T G^ accounts for the change of energy momentum 
G M of the external field, while the radiation damping 
term ~T p ll (d T p^d T p' / )/(m 2 c 2 ) — —d T H^ accounts for 
the change of energy-momentum of the scattered elec- 
tromagnetic wave (Hartemann, 2002). 

Several difficulties appear in the point electron model 
described by the Dirac-Lorentz equation (22). First, 
it admits unphysical runaway solutions with exponen- 
tially increasing acceleration, which can be eliminated 
by requiring the Dirac-Rohrlich asymptotic condition 
linir-^oo d T Pfj, — 0. Second, physical solutions presents 
acausal preacceleration, i.e., that electron momentum 
changes before an external force is suddenly applied, on 
a time scale tq. Equation (22) can be approximated 
by evaluating the RR self-force with the solution of the 
zeroth-order equation d T p^ — —eF^u 1 ' (Landau and 
Lifshitz, 1994). This yields the Landau-Lifshitz equa- 
tion which neither admits runaway solutions nor presents 
acausal preacceleration behavior. 

For the case of an electron undulating according to Eq. 
(2), it is important to define the range of parameters for 
which radiation reaction comes into play and has to be 
included in the description of the electron motion and its 
radiation. For relativistic electrons, the dominant term in 
the radiation reaction comes from the energy momentum 
transferred to the scattered electromagnetic wave [while 
for rest electrons, the radiation reaction describes the di- 
rect exchange of energy momentum between the external 
field and the scattered wave (Hartemann, 2002)]. The 
rate of energy loss v 1 for the electron can be estimated 
from mc 2 d'y/dt = — P 7 , with P 7 the average power ra- 
diated by the electron given by Eq. (13). It leads to 

l{t) = 7o/(l + V-yt) with 



T &-2 

y 7 = -^loK 



/2ttc\ 2 



(23) 



(22) 



and 70 is initial gamma factor of the electron (Esarey, 
2000; Huang and Ruth, 1998; Koga et al, 2005; Michel 
et at, 2006; Telnov, 1997). Therefore, the radiation re- 
action can be neglected when the interaction duration 
or equivalently the number of oscillations satisfy respec- 
tively, r < v- 1 and N < N RR = X u / '(2tt 2 CT -f K 2 ) . 
With conventional undulators (see Sec. V), RR will al- 
ways be negligible; for ~ 10 GeV electrons, X u ~ 1 cm 
and K ~ 1, the limiting number of period A?rr is on the 
order of 1.3 x 10 7 . For current and short-term laser-based 
betatron experiments (see Sec. IV), RR is negligible but 
in the long term or for electron beam driven plasma ac- 
celerators [e.g., parameters of the Facility for Advanced 
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Accelerator Experimental Tests (FACET) (Hogan et ai, 
2010)], - 10 GeV electrons, A„ ~ 1 cm, and K ~ 100 
lead to Nrr ~ 1.3 x 10 3 . For Thomson backscattering 
(see Sec. VI), GeV electrons colliding with a laser pulse of 
strength parameter K ~ 10 and wavelength A = 0.8 /xm 
(A u = A/2) leads to -/Vrr ~ 60. The radiation reaction 
can be neglected in Thomson backscattering for sub-GeV 
electron beams and laser pulses of strength parameter on 
the order of unity, for which N-rr > 6 x 10 3 . 

In the realm of quantum electrodynamics (QED), the 
radiation reaction corresponds to the recoil experienced 
by an electron due to consecutive incoherent photon 
emissions (Di Piazza et ai, 2010). Quantum effects be- 
come important when the electron energy loss associ- 
ated with the emission of a photon is on the order of 
the electron energy. Signatures of quantum effects can 
be observed before entering this quantum regime. In- 
deed, quantum fluctuations, which imply that different 
electrons emit a different number of photons (that carry 
different energies) and hence lose a different amount of 
energy, can lead to an observable increase of the elec- 
tron beam energy spread (Esarey, 2000). Experimen- 
tally, a study of the quantum regime is accessible in the 
framework of Compton scattering (Sec. VI). QED ef- 
fects, such as nonlinear Compton scattering (Bula et ai, 
1996) and the production of electron-positron pairs from 
light (Burke et ai, 1997), were observed in the SLAC 
E-144 experiment, where 46.6 GeV electron beams col- 
lided with relativistic laser pulses with intensities of 10 18 
W.cm~ 2 . 



H. Real electron bunch: longitudinal and transverse 
emittance 

In the previous section, ideal bunches with electrons at 
the same energy and same momentum have been consid- 
ered in order to simply discuss the effect of summation 
of the radiation from each electron. However, in realis- 
tic bunches, electrons have slightly different energies and 
momenta. More precisely, inside the bunch, electrons 
that are at the same location can have different energies 
and momenta. This implies, for example, that the bunch 
cannot be focused and compressed on an infinitely small 
point. This limitation is fundamental and inherent to the 
bunch, it does not depend on practical realization. The 
parameter which accounts for that is called the emittance 
and is related to the volume occupied by the electrons 
in the 6D phase space (x,y, z,p x ,p y ,p z ) at a given time 
(Humphries, 1990). The 6D phase volume is constant 
in time if only smooth external forces are applied and if 
collisions are neglected (the phase volume conservation is 
a consequence of the collisionless Boltzmann equation). 
The emittance reflects the quality of the electron beam 
because it quantitatively indicates if electrons have the 
same coordinates, direction, and energy. 



For a relativistic electron beam traveling in the e z di- 
rection, an emittance is defined for each dimension: the 
z one is called longitudinal and two others are trans- 
verse (x and y). For a uniform distribution with sharp 
boundary, the normalized emittance e a jy is defined as 
the area occupied by electrons in the (a,p a /mc) space 
divided by ir (a — x,y,z). But because realistic elec- 
tron beams have diffuse boundaries, the normalized root- 
mean-square (rms) emittance is used and defined as 

e aN = V(Aa 2 ) (ApJ) - (AaAp^/mc, (24) 

with a — x,y, z and where Aa = a — (a), Ap a = p a — (p a ) . 
For transverse dimensions, it is convenient to use the 
unnormalized emittance e a , which is related to the area 
occupied by electrons in the (a, a') trace space (a = x, y), 
where a' ~ Pa/Pz is the transverse angle with respect to 
the propagation axis z, 

e a = v^Aa 2 ) (Aa' 2 ) - (AaAa') 2 , a = x,y. 

.(25) 

It is generally expressed in 7r.mm.mrad. For cylindrically 
symmetric beams, the emittance e r (defined in the same 
way by putting a — r) can be used. The normalized 
omittance is related to the unnormalized one by ejy = 
7/3e and has the advantage of being conserved during 
acceleration (in systems which preserve the emittance). 
In a focal plane, where there is no correlation between 
position and angle, the emittance is simply the product of 
the rms transverse size o~ a by the rms angular dispersion 
cr a >- 

e a = a a a a >, a = x,y,r. (26) 

For incoherent radiation of electrons oscillating in un- 
dulator or wiggler devices, the emittance and energy 
spread have the following effects. The electron beam 
is focused in the device, with a transverse size a and 
a divergence 9 satisfying e = <j9. Because of the angu- 
lar spread, the radiation angular distributions from each 
single electron are slightly shifted from one another, lead- 
ing to a redshifted broadening of harmonic bandwidths 
[an electron with direction 9 with respect to the prop- 
agation axis contributes on axis with higher wavelength 
A = A e=0 (l +7 2 # 2 ), see Eq. (5)]. The energy spread 
effect is straightforward: electrons with different ener- 
gies radiate at slightly different wavelengths, leading to a 
broadening of harmonic bandwidths. These effects re- 
sult in a modified bandwidth given by (AA n /A„) 2 = 
(l/(nN)) 2 + (2A 7 / 7 ) 2 + ( 7 2 e 2 /cr 2 ) 2 . Hence, the band- 
width of an harmonic at a given direction comes from 
three different effects: the finite number of periods, the 
energy spread, and the angular spread (which depends 
on the emittance). 

The transverse emittance is essential for transport con- 
siderations and applications such as the free-electron 
laser. Concerning the FEL application, required condi- 
tions on the transverse emittance and the energy spread 
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will be given in Sec. VII. A smaller transverse emittance 
permits one to transport or focus the electron beam on 
a smaller focal spot size. 



III. ELECTRON ACCELERATION IN PLASMA 



plasma wave. After a brief description of the character- 
istics of the acceleration mechanism, recent experimental 
progress is presented. We refer the interested reader to 
the recent article of Esarey et al. (2009) for a complete 
review of laser-plasma electron accelerators. 



The possibility to accelerate electrons in laser- 
produced plasmas was originally proposed by Tajima and 
Dawson (1979). They suggested to use the intense elec- 
tric field of a relativistic plasma wave, created in the wake 
of an intense laser pulse, to accelerate electrons to rela- 
tivistic energies. The main advantage of plasmas relies 
on their ability to sustain an accelerating gradient much 
larger (on the order of 100 GeV/m) than a conventional 
radio frequency accelerating module (on the order of 10 
MeV/m). This means that electrons could be accelerated 
up to 1 GeV in millimeter- or centimeter-scale plasmas 
(Clayton et al, 2010; Froula et al, 2009; Hafz et al, 2008; 
Kncip et al, 2009; Leemans et al, 2006) while a few tens 
of meters would be necessary to reach the same energy 
in conventional accelerators. 

This acceleration method has experienced a remark- 
able development over the past decades, mainly thanks 
to the advent of high-intensity lasers and to a better un- 
derstanding of the physical mechanisms driving the accel- 
eration. Different plasma accelerator schemes have been 
developed over the years, 3 leading to electron bunches 
with ever-increasing quality. The most efficient to date 
is the so-called bubble, blowout or cavitated wakefield 
regime (Lu et al, 2006a,b, 2007; Pukhov et al, 2004a; 
Pukhov and Meyer-ter Vehn, 2002). In that regime, de- 
pending on the chosen parameters, electron bunches can 
now be produced with tunable energy in the hundreds of 
MeV range (Faure et al, 2006), low divergence (mrad), a 
relatively high charge (~ 100 pC), and a bunch duration 
of less than 10 fs (Davoine et al, 2008; Faure et al, 2004; 
Geddes et al, 2004; Glinec et al, 2007; Lundh et al, 
2011; Mangles et al, 2004, 2006; Thomas et al, 2007; 
van Tilborg et al, 2006; Tsung et al, 2006, 2004). Be- 
cause the x-ray sources which will be reviewed are based 
on laser-plasma accelerators, this section is dedicated to 
a short description of wakefield acceleration in the cavi- 
tated regime. In particular two important physical mech- 
anisms are introduced: the pondcromotive force and the 



3 For an overview of the historical development of the field, see 
Amiranoff et al. (1998, 1995); Bingham et al. (2004); Clayton 
et al. (1993); Coverdale et al. (1995); Esarey et al. (2009, 1996b); 
Everett et al. (1994); Gahn et al. (1999); Gordon et al. (1998); 
Joshi (2007); Joshi et al. (1984); Kitagawa et al. (1992); Leemans 
et al. (2002); Malka et al. (2002); Mangles et al. (2005); Mod- 
cna et al. (1995); Moore et al. (1997); Najmudin et al. (2003); 
Nakajima et al. (1995); Patcl (2007); Pukhov (2003); Santala 
et al. (2001); Tajima and Dawson (1979); Ting et al. (1997); 
Umstadter et al. (1996a); and Wagner et al. (1997). 



A. Ponderomotive force and plasma waves 

The ponderomotive force is a force associated with the 
intensity gradients in the laser pulse, that pushes both 
electrons and ions out of the high-intensity regions. Ions, 
being much heavier than electrons, still remain for short 
interaction times whereas electrons are cast away. This 
leads, in an underdense plasma, to the formation of a 
relativistic plasma wave whose fields can accelerate elec- 
trons. Here a short description of the ponderomotive 
force and of the excitation of a plasma wave (Kruer, 1988) 
is given. 

Because of the mass of the plasma ions, they can be 
considered motionless for short interaction times. Con- 
sidering a fluid description for the plasma electrons, the 
equation of motion for a fluid element submitted to the 
electromagnetic force reads 4 



dp da -* . , , 
dt = dt +C ^~^ 



(27) 



where p = p/mc is the normalized momentum of an elec- 
tron fluid element, 7 is the relativistic factor of an elec- 
tron fluid element, and cj) — eV/mc 2 and a = eA/mc 
are, respectively, the normalized scalar potential and the 
normalized vector potential of the electromagnetic fields. 
a describes the high-frequency laser pulse, cV</> is the 
Coulomb force associated with the charge distribution, 
and — CV7 is the relativistic ponderomotive force which 
expels electrons away from the laser pulse. In the absence 
of Coulomb and ponderomotive forces, the equation sim- 
plifies to p = a, corresponding to the fast electron oscilla- 
tion in the laser pulse. Depending on the amplitude of the 
laser pulse normalized vector potential do, plasma elec- 
trons oscillate with relativistic velocities \v\ ~ c (do > 1) 
or with velocities much smaller than c (ao <C 1). An 
inhomogenous laser intensity distribution leads to an in- 
homogenous 7 distribution and to a ponderomotive force 
that pushes plasma electrons from the high 7 region (cor- 
responding to high do) to the low 7 region (low ao). This 
slow drift motion of the plasma electrons leads to a charge 
density distribution p responsible for a Coulomb force 
cV0 (by virtue of the Poisson equation A</> = p/eo). 

On the other hand, a small charge density perturba- 
tion in a plasma oscillates at a characteristic frequency, 



4 Equation (27) assumes that V X (p — a) is initially zero, which 
is the case in practice because both p and a are zero before the 



passage of the laser pulse in the plasma. 
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the plasma frequency u> p = \J n e e 2 /meo, where n e is 
the electron density of the plasma. For example, if we 
translate an electron slice (continuously, without cross- 
ing between different electrons) from its initial position, 
the slice will oscillate around its initial position at the 
plasma frequency because of the restoring Coulomb force 
from the plasma ions. These plasma oscillations, called 
plasma waves, are excited by the ponderomotive force of 
the laser pulse, which creates the initial charge density 
perturbation. The phase velocity of the plasma wave v^, 
excited in the wake of the laser pulse, is approximately 
equal to the laser group velocity. 

For a low-intensity laser pulse, with a normalized vec- 
tor potential amplitude a < 1, the excited plasma wave 
is linear and has a sinusoidal shape at the plasma fre- 
quency, 4>{r,i) = tf)o(r) sm(u> p t). For a higher intensity, 
ao > 1, the plasma wave becomes highly nonlinear and 
can involve transverse currents and a quasistatic mag- 
netic field (Gorbunov et at, 1996). 

Relativistic plasma waves can be produced in various 
regimes depending on the laser and plasma parameters. 
For appropriately chosen parameters, electrons can be 
trapped at a proper phase of the plasma wave and ex- 
perience its field over a distance sufficiently long to be 
accelerated up to relativistic energies. 

In practical units, the plasma wavelength X p = 2ttc/lo p 
and the laser strength parameter ao are given by 

\p\pm] = 3.34 x 10 10 /V«e[cm- 3 ], (28) 
a = 0.855^ I[10 ls W/cm 2 ]A|[/im], (29) 

where / is the laser intensity and Al is the laser wave- 
length. 



B. The cavitated wakefield or bubble regime 

To date, the most efficient mechanism to accelerate 
electrons in a plasma wave is called the bubble, blowout 
or cavitated wakefield regime. In this regime the wake 
consists of an ion cavity having a spherical shape (Lu 
et al, 2006a,b, 2007; Pukhov and Meyer-ter Vchn, 2002). 

This regime is reached when the waist wo of the focused 
laser pulse becomes matched to the plasma (k p w — 
2^/oq, with k p = uj p /c) and if the pulse duration is of the 
order of half a plasma wavelength (cr ~ A p /2). In addi- 
tion, the laser intensity must be sufficiently high (ao > 2) 
to expel most of the electrons out of the focal spot. If 
these conditions are met, an ion cavity is formed in the 
wake of the laser pulse, as represented in Fig. 4. Elec- 
trons can be trapped at the back of the cavity and accel- 
erated by the high electric field (the space-charge force) 
until they reach the middle of the cavity where they start 
to decelerate. The distance over which electrons must 



F z 
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FIG. 4 Principle of the electron acceleration in the bubble 
regime. The laser pulse expels all plasma electrons out of 
the focal spot, leaving in its wake an ion cavity of radius rv 
The longitudinal force F z accelerates the self-injected elec- 
tron bunch in the first half of the bubble, while the bunch 
decelerates in the second half of the bubble. 

propagate before they reach that point is called the de- 
phasing length Ld and is much larger than the bubble 
radius r^, as electrons travel almost at the same speed 
as the wave phase velocity. If the process is turned off 
at that time, by setting the acceleration length close to 
the dephasing length, electrons exit the plasma with the 
maximum energy gain. The maximum electric field and 
the radius of the cavity are, respectively (Lu et al. , 2007) , 

E m = muj p cy/a^/e, (30) 
n = w a = (2/k p )y/aQ. (31) 

For the typical experimental condition accessible with 
present lasers, ao = 4 with a 30 fs laser, the maximum 
electric field is E m ~ 600 GeV/m and the radius of the 
cavity is ri, ~ 7 /im for an electron density of n e = 1 x 
10 19 cm~ 3 . 



C. Experimental production of relativistic electron bunches 

Several techniques have been explored to produce 
plasma waves appropriate to accelerate electrons. A 
typical experiment for laser-plasma acceleration consists 
in focusing a laser pulse into a gas jet, the interac- 
tion parameters being essentially the laser intensity, fo- 
cal spot size and duration, the propagation length, and 
the plasma density. Depending on the choice of these pa- 
rameters the features of the produced electron bunch can 
be very different. Prior to 2004, experiments used rela- 
tively high plasma densities (> 10 19 cm 3 ) and the pro- 
duced electron bunches were characterized by broadband 
spectra, extending up to about 100 MeV. These spectra 
were either nearly Maxwellian in the direct laser accel- 
eration (DLA) regime (Gahn et al., 1999; Kncip et al., 
2008; Mangles et al, 2005; Pukhov, 2003; Pukhov et al, 
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1999; Tsakiris et al, 2000) or non-Maxwellian in the self- 
modulated laser wakefield accelerator (SM-LWFA) and 
forced laser wakefield (FLWF) regimes (Coverdale et al, 
1995; Malka et al, 2002; Modena et al, 1995; Najmudin 
et al, 2003; Santala et al, 2001; Ting et al, 1997; Um- 
stadter et al, 1996a; Wagner et al, 1997). In 2004, ma- 
jor advances were made by setting the laser pulse dura- 
tion close to the plasma period, increasing the interaction 
length and matching the dephasing to the propagation 
length. Three groups reported simultaneously on the 
production of monoenergetic electrons in the 100 MeV 
range (Faure et al, 2004; Gcddes et al, 2004; Mangles 
et al, 2004), collimated within a few milliradians and 
with charge on the order of 100 pC. However, despite 
these remarkable progresses, the high nonlinearity of the 
mechanism resulted in electron bunches which were nei- 
ther stable nor tunable in energy; in addition electron 
energies remained below the GeV level. 

To overcome these limitations, improved schemes have 
recently been developed. An external and controlled elec- 
tron injection into the wakefield was proposed (Esarey 
et al, 1997; Fubiani et al, 2004; Umstadter et al, 1996b) 
and recently demonstrated (Faure et al, 2007, 2006; 
Malka et al, 2009). It consists of colliding the main laser 
pulse generating the plasma wave with a second laser 
pulse, creating a beat wave whose ponderomotive force 
can preaccelerate and locally inject background electrons 
in the wakefield (Davoinc et al, 2008). The energy of 
the electron bunch can be tuned by varying the collision 
position and therefore the acceleration length. Exper- 
iments demonstrated that stable electron bunches can 
be produced with energy continuously tunable from a 
few tens of MeV to above 200 MeV (Faure et al, 2007, 
2006; Malka et al, 2009) and an energy spread on the 
order of 1% (Rechatin et al, 2009b). The few femtosec- 
ond duration and the few kA current of these electron 
bunches were experimentally demonstrated by Lundh 
et al. (2011). A cold optical injection providing nar- 
row energy spread for GeV electrons has been proposed 
(Davoine et al, 2009). 

Several other possibilities to control electron injection 
have been investigated: the use of a plasma density down- 
ramp (Brantov et al, 2008; Bulanov et al, 1998; Faure 
et al, 2010; Gcddes et al, 2008; Hemker et al, 2002; Suk 
et al, 2001), ionization-induced injection (Clayton et al, 
2010; McGuffey et al, 2010; Pak et al, 2010; Pollock 
et al, 2011) or magnetically controlled injection (Vicira 
et al, 2011). 

Laser-plasma accelerators have also recently reached the 
GeV level, using either external laser guiding (Leemans 
et al, 2006; Nakamura et al, 2007) or higher laser power 
(Clayton et al, 2010; Froula et al, 2009; Hafz et al, 2008; 
Kneip et al, 2009). In the first case, the experiment re- 
lied on channeling a few tens of TW-class laser pulses 
in a gas-filled capillary discharge waveguide in order to 
increase the propagation distance, and so the accelera- 



tor length, to the centimeter scale (Leemans et al, 2006; 
Nakamura et al, 2007; Rowlands-Rees et al, 2008). In 
the experiment of Leemans et al. (2006), electrons have 
been accelerated up to a GeV. 

In order to further improve the quality of electrons 
from laser-plasma accelerators, several routes are pro- 
posed from the use of a PW-class laser to multistaged 
acceleration schemes (Gordicnko and Pukhov, 2005; Lif- 
schitz et al, 2005; Lu et al, 2007; Malka et al, 2006; 
Martins et al, 2010). These foreseen developments are 
of major importance for the production of a free-electron 
laser based on a laser-plasma accelerator. 



IV. PLASMA ACCELERATOR AND PLASMA 
UNDULATOR: BETATRON RADIATION 

In the bubble acceleration regime, the plasma cavity 
can act as a wiggler in addition to being an accelera- 
tor, reproducing on a millimeter scale the principle of 
a synchrotron to produce x rays (Kiselev et al, 2004; 
Rousse et al, 2004). In this section, we will show a laser- 
produced ion cavity drives the electron orbits in such a 
way that a short pulse of collimated x-ray radiation is 
emitted. Figure 5 represents the principle of the mecha- 
nism. As discussed, the bubble regime is reached when an 
intense femtosecond laser pulse, propagating in an under- 
dense plasma, evacuates plasma electrons from the high- 
intensity regions and leaves an ion cavity in its wake. In 
addition to the longitudinal force, responsible for the ac- 
celeration discussed above, the spherical shape of the ion 
cavity results in a transverse electric field producing a 
restoring force directed toward the laser pulse propaga- 
tion axis. Therefore, electrons trapped and accelerated in 
the cavity are also transversally wiggled. The conditions 
for an efficient production of accelerated charged particle 
radiation, discussed in Sec. II, are therefore met. A colli- 
mated beam of x-ray radiation is emitted by the electron 
bunch. This radiation, which can be directly compared 
to a synchrotron emission in the wiggler regime, is called 
betatron radiation. 

The betatron radiation from laser-produced plasmas 
was simultaneously proposed and demonstrated in 2004 
(Kiselev et al, 2004; Rousse et al, 2004). This repre- 
sented a major step forward in the field of plasma x- 
ray sources since it was the first method allowing one to 
produce bright collimated x-ray (keV) beams from laser- 
plasma interactions. Since then, this radiation has been 
measured and widely characterized in interaction regimes 
from multiterawatt (Albert et al, 2008; Shah et al, 2006; 
Ta Phuoc et al, 2005a, 2008a, 2006, 2007) to petawatt 
lasers (Kneip et al, 2008). According to the theory, fem- 
tosecond pulses of x rays up to a few tens of keV could 
be produced within tens of milliradian beams. In the 
following sections, the properties of the betatron mech- 
anism are reviewed. Following the approach of Sec. II, 
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of a test electron in the cavity is then: 



Betatron 
orbit 



Electron 



FIG. 5 Schematic of the betatron mechanism. When an elec- 
tron is injected in the ion cavity, it is submitted not only 
to the accelerating force, but also to a restoring transverse 
force, resulting in its wiggling around the propagation axis. 
Because of this motion, the electron radiates x rays whose 
typical divergence is 8 = if/7. 



the electron orbit is first calculated using an ideal model 
and then the features of the emitted radiation are de- 
rived. Simulations based on this model and a particle 
in cell (PIC) simulation are presented. Finally, after a 
summary of the experimental results, we conclude with 
the short term developments foreseen. 



A. Electron orbit in an ion cavity 

An idealized model of a wakefield in the bubble regime 
based on the phenomenological description developed by 
Lu et al. (2006a, b, 2007) is assumed. The phase ve- 
locity of the wake is expressed as = v g — w c tch — 



1 



- ujp/ivl is the laser 



c(l — 3Wp/2cj£), where v g — 

group velocity in the plasma, w c t c h — cWp/ w £ is the etch- 
ing velocity due to local pump depletion (Decker and 
Mori, 1 994, 1995; Decker et al, 1996; Lu et al, 2007), 
Lu p = n e e 2 /meo is the plasma frequency, n e the elec- 
tron density of the plasma, and ljl the central frequency 
of the laser field. The ion cavity is assumed to be a sphere 
of radius rt = (2/k p )^/ao {k p = ui p /c) and a cylindrical 
coordinate system (r,9,z), where r is the distance from 
the laser pulse propagation axis is used. The comoving 
variable is defined as £ = z — v^t so that (£ = 0, r = 0), 
(C = rb,r = 0), and (C = — r&,r — 0), respectively, cor- 
respond to the center, the front, and the back of the 
cavity. The electromagnetic fields of the wake are an 
axial electric field E z , a radial electric field E r , and an 
azimuthal magnetic field Bg. They are given by E z /Eq = 
k P C/2, E r /E Q = k p r/A, and Bg/Eo = -k p r/(Ac), where 
E a = moj p c/e is the cold nonrelativistic wavebreaking 
field (Kostyukov et al, 2004; Lu et al, 2006a,b; Pukhov 
et al, 2004a; Xie et al, 2007). The equation of motion 



d ^=-e(E + vxB) 
at 



F« +F, 



(32) 

where p is the momentum of the electron. The last ex- 
pression assumes p± <C p z . The electron is initially in- 
jected at the back of the cavity with space-time coor- 
dinates (ti = Q,Xi,yi,Zi) such that x\ + yf + zf = r\ 
and Ci < 0, and with a energy-momentum quadrivector 
( r yimc 2 ,pi) such that v z > v$. 

The term F\\ is responsible for the electron acceleration 
in the longitudinal direction e z . As the electron becomes 
relativistic, its velocity becomes greater than v^; the term 
— C — v <j)t ~ z decreases and the accelerating force is re- 
duced. The length for the electron to reach the middle 
(C = 0) of the cavity and to become decelerated (C > 0) 
corresponds to the dephasing length L^. The term F± is 
a linear restoring force that drives the transverse oscilla- 
tions of the electron across the cavity axis at the betatron 
frequency ujp ~ u) p / \f2q. 

To derive the analytical expressions of the electron or- 
bit, we start the integration of the equation of motion 
from an initial state where 7^ = 1/yl — (v Z i/c) 2 70 
in order to ensure that the slippage between the test elec- 
tron and the wakefield is constant in time. This approx- 
imation is realistic because an electron quickly attains 
lz 3> l<j> and the transitory period in which the approx- 
imation fails is small compared to the dephasing time 
Ld/c. Hence d(/dt — v z ~ ~ c — v^, ~ 0/(27?) if 
7i > 1, In agreement with the assumption used in the 
last expression of Eq. (32), we perform a perturbative 
treatment in the variable p±/p z -C 1. The variables with 
a hat are normalized by the choice m = c = e = uj p = 1 
and p = v/c is the velocity normalized to the speed of 
light c. Equation (32) projected on each axis reads 



di 

di 

d(M 
di 



x 
2' 

V 
2' 

2 



t 



(33) 
(34) 

(35) 



At zero order (pj_ = 0), Eq. 
grated to 7(f) ~ p z (i) = 7d(l- 

r = (t- 



(35) can be directly inte- 
-r 2 ) with 7^ = ^+7^2/2, 

27?^ the dephasing 



■id)/y£d + 87^ li and i d 

time {id = k p Ld). This parabolic profile of 7(7") can be 
inserted into Eqs. (33)-(34) to obtain the first-order so- 
lution for x(t) and y(r). Note that motion in the e x and 
e v directions are decoupled at this order of calculation. 
Each equation takes the form of a Legendre differential 
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equation, 



d 
dr 


a-r 2 ; 


dx' 
Tr\ 


d 


\l r 2 


f] 


dr 


dri 



+ v(v + l)x = 0, 



(36) 
(37) 



where v{y + 1) = 47?. The solution space is generated 
by the Legendre functions of the first and second kinds 
P v {t) and Q v (t), which have the following asymptotic 
expressions for v ^> 1 (Abramowitz and Stegun, 1972): 

PAt) = 1/ ^(1-^)1/4 cos ^ + V2) cos- 1 r - tt/4], 

(38) 

[(i/+l/2) cos _1 T + 7r/4]. 

(39) 



Q*(r) 



2i/(l -t 2 )V4 



Hence, the solutions for x (r) and z/(r) can be written as 



£(t) 



(1 -T 2 )!/4 
(1 - T 2 )V4 



3(27^ cos 1 r + <p iB ), (40) 
3(27^ cos -1 r + <p y ), (41) 



where (A x , tp x ) and (^4 y , (/?y) are constants which have to 
be determined by initial conditions (xi,p X j) and (yi,p y i). 
They can be put into the simpler form 

x(t) = Xf} {i) cos f / u>f)(P)dP + <p x ) , (42) 



y(t) = yp(t) cos / Ljp{t')dt' +<f/ v ], (43) 



Xf)(t) 

hit) 



^/7(^) 1/4 , 

4/7(^) 1/4 - 



(44) 
(45) 
(46) 



The transverse motion consists of sinusoidal oscilla- 
tions in each direction with time-dependent amplitude 
Xp{t),yp(t), and frequency u)p(t). They depend on t only 
through 7. For an arbitrary 7(f) profile, Eqs. (42)-(46) 
can be derived from the adiabatic approximation and the 
WKB method [valid if (l/uip)du)p/dt <C 1], using the con- 
servation of the action variables (Kostyukov et al., 2004; 
Thomas, 2010). Thus, Eqs. (44)-(46) in terms of 7 are 
more general than the frame of the above calculation [i.e., 
the case of a parabolic j(t) profile]. The betatron ampli- 
tudes xp,yp decrease as 7- 1 / 4 during acceleration, while 
the betatron frequency uip decreases as 7- 1 / 2 . Accord- 
ing to Eqs. (40)-(41), the number of betatron oscillations 
between t = and the dephasing time t = is approx- 
imately 7^/2. Depending on the values of ip' x and ip' , 
the motion in each direction e x and e y can be in phase 



(motion confined in a plane) or not (helical motion) . The 
latter case corresponds to an initial state with nonzero 
angular momentum L z — xpy — yPx • 

Finally, the lon- 
gitudinal motion z(t) can be obtained by integrating f} x 
which is deduced from 7 2 = 1/(1 — fi 2 x ~ /3 2 — /3 2 ). 

This trajectory is noticeably different from the one 
studied in Sec. II. First, the betatron amplitude and 
frequency are time dependent. Nevertheless, in the wig- 
gler regime, the resulting radiated energy dl/du can 
be seen as an integral over time of the radiated power 
dP/doj(ujp, xp,yp) for each instantaneous amplitude and 
frequency if the characteristic time scale of acceleration is 
much longer than the betatron period. Second, x(t) and 
y(t) are sinusoidal functions of t and not z. Even if the 
detail of the electron dynamics is slightly different from 
the case of Sec. II, the betatron trajectory has the same 
properties. As in Sec. II, the trajectory is a figure-eight 
motion in the electron local average rest frame. The radi- 
ation features derived in Sec. II from the parameters K, 
A M , and 7 remain therefore valid. Assuming dt ~ dz/c in 
Eqs. (42)-(43) and identifying the terms with the ones 
of Eq. (2), the local electron period X u (t) and the local 
strength parameter K(t) read in physical units 



\ u {t) = V2^{tjx P , 

K(t) = r (t)k p ^(t)/2, 



(47) 
(48) 



in which we assume a motion confined in a plane with 
a betatron amplitude rp(t) = yjxp(t) 2 + yp{t) 2 . We re- 
cover the results obtained by Esarey et al. (2002) and 
Kostyukov et al. (2003) in the case of an ion channel 
without longitudinal acceleration, but in which the pa- 
rameters are time dependent through the variation of 7. 

In the articles of Esarey et al. (2002) and Kostyukov 
et al. (2003), the constant of motion ^ z = -Jl + p 2 z ap- 
pears in Eqs. (47) and (48) instead of 7. The relative dif- 
ference between 7 and j z is of the order of p\ /p 2 and can 
be neglected at the order where the trajectory has been 
calculated. In this section, the definition 7 Z = yT + p 2 
has not been used because it is different from the usual 
one •fz = 1 / \J 1 — f3 z (which is small compared to 7 for 
wigglers, as discussed in Sec. II. C), commonly used in 
the synchrotron and FEL communities. 
In practical units, Eqs. (47) and (48) read 



X u [fim] = 4.72 x 10 10 vVMcm-' 3 ], (49) 
K = 1.33 x 10~ 10 vVUcm- 3 ]^ {pan] . (50) 



In the equation of motion (32), the electromagnetic 
fields correspond to an idealized cavitated regime. It as- 
sumes three conditions: (i) the matching conditions in 
terms of focal spot size, laser pulse duration, and den- 
sity are perfectly met; (u) the cavity is free of electrons; 
and (m) the trapped electrons are only submitted to the 
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wakefields (there is no interaction with the electromag- 
netic fields of the laser pulse). In addition, in this ap- 
proach, initial conditions have to be arbitrarily chosen 
since the injection mechanism is not taken into account. 
If the three above conditions are not met, the following 
effects occur. 

First, the interaction regime is slightly different when 
the plasma wavelength is smaller than the focal spot size 
and the pulse duration (slightly off the matching condi- 
tions). In this regime, called the forced laser wakefield 
regime (Malka et al, 2002), the shape of the wake is not 
spherical anymore, the laser pulse does not necessarily ex- 
pel all electrons from the focal spot. The wakefields are 
therefore different from the case of a cavity free of elec- 
trons. Nevertheless, the bubble model remains a good 
approximation to roughly determine the electron trajec- 
tories and the betatron radiation features. 

Second, the cavity cannot be considered free of elec- 
trons when injection occurs. Trapped electrons can mod- 
ify the fields of the cavity by the beam loading mechanism 
(Rechatin et al, 2009a; Tzoufras et al, 2008, 2009): they 
modify the motion of background electrons and hence the 
resulting wakefields. As a result the acceleration is re- 
duced and, for massive electron injection, the term Fy in 
(32) can be largely overestimated. 

Third, electrons that have been accelerated can catch 
up with the back of the laser pulse. The interaction of 
an electron with the laser pulse can increase its oscilla- 
tion amplitude rp in the direction of the laser polariza- 
tion, enhancing the betatron radiation. The interaction 
is complex and strongly dependent on the phase velocity 
of the laser pulse. It corresponds to the mechanism of 
direct laser acceleration (DLA), in which some electrons 
can be in betatron resonance, if the induced transverse 
motion of the laser pulse fields is in the same direction 
as the betatron motion at all times (Gahn et al, 1999; 
Kneip et at, 2008; Mangles et al, 2005; Pukhov, 2003; 
Pukhov et al, 1999; Tsakiris et al, 2000). In the forced 
laser wakefield regime, electrons can catch up with the 
laser pulse more easily, resulting in a brighter betatron 
radiation. 

Three-dimensional PIC simulations can be used to ac- 
curately describe al these effects, and to compute the 
electron orbits and the emitted radiation. Simple nu- 
merical simulation, integrating the equation of motion, 
can be used as well to provide the additional effects by 
modeling the laser pulse and by extracting effective ac- 
celeration and transverse fields from a PIC simulation. 
Nemeth et al (2008) followed this approach and showed 
that the betatron motion can be driven by the laser pulse 
in the polarization direction if electrons interact with it. 
Such a model, without the laser pulse fields, has also been 
considered by Wu et al. (2009) who, for a given config- 
uration with ultrahigh intensity (a = 20), reduced the 
acceleration field by a factor of 0.35 and the transverse 
fields by a factor of 0.9 in order to fit the electron tra- 



jectories with the ones extracted directly from the PIC 
simulation. This highlights how the wakefields are over- 
estimated when Eq. (32) is considered without any cor- 
rection. 

B. Radiation properties 

The betatron motion derived in the above simple 
model leads to the emission of synchrotron radiation re- 
ferred to as betatron radiation. It can be calculated 
within the general formalism of the radiation from a mov- 
ing charge. The ion cavity acts as an undulator or a 
wiggler with a period X u {t) and strength parameter K(t) 
which depends on the electron initial conditions upon in- 
jection into the cavity. We first consider the radiation 
properties for constant A u , K, and 7 and then consider 
the radiation produced by the electron with acceleration. 

1. Without acceleration 

The spectrum of the emitted radiation depends on the 
amplitude of K. For a small amplitude of the betatron 
oscillation K -C 1, the radiation is emitted at the fun- 
damental photon energy Hui with a narrow bandwidth in 
the forward direction (8 — 0). As K 1 harmonics of 
the fundamental start to appear in the spectrum, and for 
K ^> 1 the spectrum contains many closely spaced har- 
monics and extends up to a critical energy fko c . These 
quantities are given by 

hu = (2j 2 hc/\ u )/(l + K 2 /2) fovK<l, 

hu c = ^Kj 2 hc/X u for K » 1, (51) 

which gives in practical units 

huj[eV] = 5.25 x lO^V' V«e[cm- 3 ] forAT<l, 

hu c [eV] = 5.24 x 10" 21 7 2 n e [cm" 3 ]r^[//m] for K » I. 

(52) 

According to Sec. II, the radiation is collimated within 
a cone of typical opening angle 6 r = I/7 in the undulator 
case. For a wiggler, the radiation is collimated within a 
typical opening angle K/j in the electron motion plane 
i^xi^z) an d I/7 in the orthogonal plane (e y ,e z ). 

The x-ray pulse duration equals the electron bunch du- 
ration. Depending on the parameters and the electron in- 
jection process, it can be extremely short, a few fs (Lundh 
et al, 2011), corresponding to a bunch length equal to a 
very small fraction of the bubble radius, or using massive 
self-injection at high density, the bunch length can be on 
the order of the bubble radius and r r ~ r^J c. 

The number of emitted photons can be estimated using 
the expressions of Sec. II. In practical units, the number 
of photons emitted per period and per electron (at the 
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mean photon energy (hui) 
is 



0.3huj c for the wiggler case) 



1.53 x 10~ 2 K 2 forK<l. 



N~ = 3.31 x 10~ 2 K for K » 1. 



(53) 



From the above expressions, an estimation of the radia- 
tion properties for constant X u , K, and 7 can be obtained 
for a typical parameter regime. A 100 MeV electron 
(7 ~ 200) undergoing betatron oscillations in a plasma 
of density n e = 2 x 10 19 cm~ 3 is considered. The spatial 
period of the motion is A„ ~ 150 /im. For K <C 1, this 
electron will emit betatron radiation at a fundamental 
energy fvjj ~ 650 eV. In the same parameter regime but 
with K ~ 10, typical of our experimental conditions, the 
critical energy of the radiation is ftw c ~ 5 keV. Along the 
interaction length L acc , the electron accelerates and ra- 
diates but its main contribution in terms of energy to the 
betatron radiation comes from the part of the trajectory 
where its energy is maximal. If the electron is around 
its maximal energy during ~ 3 betatron periods, the to- 
tal number of photons emitted per electron at the mean 
energy hcu ~ 1.5 keV is ~ 1. Considering that the 
number of electrons trapped into the ion cavity is on the 
order of 10 8-9 , the number of x-ray photons expected is 
in the range 10 8-9 as well. Finally, in a typical parameter 
regime where K ~ 10 and 7 ~ 200, the betatron emis- 
sion is collimated within a cone of typical solid angle of 
50 mrad x 5 mrad for a single electron. 



2. With acceleration 

For an electron accelerating, the situation is more com- 
plex. In what follows, we consider the wiggler limit and 
we derive an approximation for the radiation spectrum 
taking into account the acceleration. This will be useful 
for discussing the scalings of betatron properties which 
will be presented in Sec. IV. E. The acceleration has a 
characteristic parabolic profile 7(7-) = 7^(1 — t 2 ), with 

r = (t- i d )/ J¥ d + 87^ ji and i d = -2^C,i the dephasing 

time (td — k p Ld), as shown in Sec. IV. A. Thus, from 
Eqs. (13), (45)-(48), and (51), we deduce 



^(t)=u, m (1-t 2 )- 1/2 , 
K{T) = K d {l-T 2 ) 1 '\ 
Huj c (t) = fiuc, d (l - r 2 ) 7 / 4 , 
P 7 ( T ) = P 7id (l_ T 2 )3/ 2 , 



(54) 
(55) 
(56) 
(57) 



where cJp,d, K d , hu! c ^ d , and P~ ( ,d are the values of the 
parameters at the dephasing time t = 0, and P 7 (r) 
is the radiated power averaged over one oscillation pe- 
riod [whose expression is given in Eq. (13) for a pla- 
nar trajectory, and is increased by a factor of 2 for a 



helical trajectory]. We consider that the spectrum (in- 
tegrated over angles) radiated per unit time is given by 
dP/du = (Pj/u] c )S(u)/u) c ), i.e., a synchrotron spectrum 
[see Eq. (8)] in which P 7 is replaced by P™, which is ex- 
act for a helical trajectory (for which the radius of cur- 
vature p does not depend on the phase of the oscillation) 
but only approximate for a sinusoidal trajectory (because 
p varies during an oscillation and uj c corresponds to the 
minimal value of p in the oscillation). The radiation spec- 
trum then reads 



dl 

did 



~' d dP 
du 



(t)dt 



dr 



^5( W / Wc (r)) 



P~l, d td c / 



U c ,d 



S (cj/cj CjC j), 



where the function S' is defined as 
'° dr 



S'(x) = 



(1-T 2 ) 1 /' 



S[(l 



2)-7/4 a 



(58) 



(59) 



Equation (58) highlights the fact that the radiation prop- 
erties are encoded into the values of the parameters 0Jf}, d , 
Kd, 7d, fi^c.d and td- The shape of the spectrum is 
not described by the usual universal function S of the 
synchrotron spectrum anymore, but by the function S' 
defined above. This latter function takes into account 
the photons emitted at low energies during the accelera- 
tion. Indeed, Eq. (53) shows that the number of emit- 
ted photons in one period has a weak dependence on 7 
(only through K ex 7 1 / 4 ), whereas the radiated energy 
is strongly dependent on 7 because photons are emit- 
ted at much higher energies for higher 7. As a result, 
the angular distribution of the radiated energy is dom- 
inated by the electron oscillations around maximal val- 
ues of the parameters, i.e., around the dephasing time. 
Hence, to obtain rough estimations of the radiation prop- 
erties, the formulas given above for constant parameters 
can be used, in which the values at the dephasing time 
are inserted. The information on the spectrum will be 
valid at high photon energies, but will underestimate the 
number of low-energy photons. 

Finally, it is interesting to compare the scalings of 
dl/dco and A^ 7 for both cases of constant parameters and 
acceleration, 



dl_ 
dio 
dl 
du 



7— constant 



oc NpKS{u/u c 



acceleration 



oc wp t dt d KdS'(w/w Cf d) 



which leads to, respectively, 



N. 



1 oc NrK, 

1 1 7— constant ^ 



N 7 \ (xufj yd tdKd, 

' acceleration ' ' 



(60) 
(61) 

(62) 
(63) 
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FIG. 6 Electron trajectory on the left and its gamma factor 
as a function of the longitudinal coordinate z on the right. 
The set of initial conditions is a = 2 /mi, yi = 0, Zi = C» = 
—Tb + 0.3 /im = —6.42 fim, p xi = 0, p yi = 0, p zi = 25 mc 



and the laser-plasma parameters are n e 
a = 4. 



and 



for the dependence of the total number of emitted pho- 
tons, where Np is the number of betatron oscillations in 
the 7 = const case. The value uJp^td is also characteristic 
of the number of betatron oscillations for the acceleration 
case. Thus, we see that the scalings are identical for both 
situations (parametrized either by the constant parame- 
ters or by the parameters at the dephasing time). The 
only difference resides in the exact shape of the spectrum, 
given by the function S' instead of S. 



C. Numerical results 

In this section, the case of an electron experiencing 
the longitudinal acceleration force is treated numerically. 
The complete equation of the electron motion (32) is first 
integrated and then the features of the emitted radia- 
tion are calculated using the general formula of moving 
charge radiation (1). Finally, a particle in cell simulation 
of the betatron radiation taking into account all possi- 
ble effects, including the injection process, is presented. 
This latter simulation is based on the extraction of the 
trapped electron trajectories from the PIC code and the 
calculation of the corresponding radiation via Eq. (1) in 
postprocessing. 



1. Test-particle simulation 

An interaction regime accessible with tens of TW-class 
lasers and parameters typical of current experiments is 
considered. The left panel of Fig. 6 represents the orbit 
of a test electron in the ion cavity obtained by integrat- 
ing Eq. (32). The electron density is n e = 1 x 10 19 
cm -3 , the propagation length is set at 1.2 mm, and the 
laser strength parameter is do = 4 (only used to calcu- 
late the bubble radius r& = 6.72 /Ltm). The test electron 
enters at the back of the cavity with Xi — 2 /im, = 0, 
Zi = Q = -r b + 0.3^m = -6.42 /xm and p xi = 0, p yi = 0, 
p Z i — 25 mc. The electron drifts along the longitu- 
dinal direction and oscillates in the transverse direction 
(Ta Phuoc et ai, 2005a). The betatron amplitude rp 
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FIG. 7 X-ray angular distribution calculated from the trajec- 
tory of Fig. 6. The color scale represents the radiated energy 
per unit solid angle. 
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FIG. 8 Transverse trajectory (on the left) and angular profile 
of the corresponding emitted radiation (on the right, the color 
scale representing the radiated energy per unit solid angle) for 
two different cases. A: The trajectory is three dimensional and 
helical (elliptical in the transverse plane). B: The trajectory 
is planar. Both electrons are accelerated up to 7 ~ 380 and 
the density is 2 x 10 19 cm -3 . 



decreases as the electron gains energy (rp oc 7 -1 ^ 4 ) and 
presents a minimum at the maximum electron energy. 
The right panel of Fig. 6 shows the 7 factor as a function 
of the longitudinal coordinate z. The electron is rapidly 
accelerated and acquires an energy up to ~ 240 MeV. It is 
decelerated after the dephasing length has been reached. 
As discussed in Sec. IV. A, the acceleration field may be 
largely overestimated and several effects not taken into 
account here can affect the acceleration process. Thus, 
the predicted maximum electron energy can be largely 
overestimated in this test-particle simulation. The shape 
of the trajectory is determined by the initial conditions 
(xi,yi,Zi) and (p X i,Pyi,Pzi) and can take various forms: 
planar as in the present simulation, circular, or helical 
(Ta Phuoc et al, 2008a, 2006). 

The radiation has been calculated using the general 
formula (1) in which the orbit calculated above is in- 
cluded. Figure 7 represents the spatial distribution of 
the radiated energy. The transverse profile of the emit- 
ted radiation has the same symmetry as the transverse 
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FIG. 9 Betatron radiation spectrum, integrated over angles 
and corresponding to the trajectory of Fig. 6, in number of 
photons per 0.1% bandwidth and per electron. 



orbit. In our particular case, the main divergence is in the 
(sx,£z) plane and the radiation beam is confined within 
a typical angle 1/7 in the perpendicular direction. The 
typical opening angles in each direction are 8xr ~ 18 
mrad and 9y r ~ 2 mrad. The spatial distribution is a 
direct signature of the transverse electron orbits in the 
cavity (Ta Phuoc et ai, 2008a, 2006). As an example, 
the spatial distributions of the radiation produced by an 
electron undergoing two types of transverse orbits are 
represented in Fig. 8. For minimum radius of curvature 
along the helical orbit (case A), the transverse electron 
momentum is along e y , so the radiation is emitted up 
and down and the radiated energy is maximal. For max- 
imum radius of curvature, the transverse momentum is 
along e x , so the radiation is emitted left and right and the 
radiated energy is minimal. Therefore, the x-ray beam 
profile is mapping out p x /p z and p y /p z . Also note that 
the increase of the radiated energy with the electron en- 
ergy appears in the radiation angular profile. Measure- 
ment of the spatial distribution in the x-ray range can 
therefore provide information on the transverse electron 
momentum during their oscillation in the wakefield cav- 
ity, provided that the energy and p z are known. It can 
also provide some insight on the initial conditions of elec- 
tron injection that determines the type of trajectory. 

Figure 9 presents the spectrum of the radiation in- 
tegrated over the spatial distribution of Fig. 7. In- 
stead of the radiated energy primarily given by Eq. (1), 
Fig. 9 shows the number of emitted photons per 0.1% 
bandwidth and per electron. The spectrum extends 
well above the keV range with a significant number of 
x-ray photons per electron. An estimate of the total 
number of photons can be obtained on the basis of the 
number of electrons that are currently accelerated in a 
laser-plasma accelerator. Assuming 10 9 electrons gives 
~ 10 5 photons/0.1% bandwidth (BW) between 1 and 10 
keV. The total photon number is on the order of 10 9 . 



2. Particle In Cell simulation 

We now look at a 3D PIC simulation in a typical exper- 
imental parameters range. We are using the CALDER- 
Circ code (Lifschitz et ai, 2009), which, for the cost of 
a few 2D simulations, can provide the fully 3D trajecto- 
ries of particles in the plasma. Indeed, this model uses 
Fourier decomposition of the electromagnetic fields in the 
azimuthal direction. The axially symmetric (or cylin- 
drical) assumption consists in computing only the single 
mode m = 0. In order to be able to account for planar 
fields, such as the linearly polarized laser, we assume the 
so-called quasicylindrical geometry and compute only one 
additional mode, m = 1. The linearly polarized laser's 
wakefield in a ultrahigh-intensity regime is well described 
by those first two modes since fields never stray far from 
the axial symmetry. The macroparticles, each represent- 
ing an assembly of electrons, evolve in a full 3D space 
in which the 3D fields have been reconstructed from the 
two known Fourier modes. Monitoring macroparticle tra- 
jectories gives us the same information as a test-particle 
code but with realistic wakefield, energies, and injection 
properties in space and time. The injected particles in 
the code are detected by their energy. As soon as a par- 
ticle reaches 45 MeV, it is considered as being injected 
and we keep track of its trajectory for the rest of the 
simulation time. These trajectories thus begin after the 
injection time of each single particle and finish at simula- 
tion end, after the particles come out of the plasma. Ap- 
proximately 10 5 macroparticles are tracked. Their three 
position coordinates and their 7 factor are recorded ev- 
ery A^traj = 100 dt comput ~ 5 fs, where dt comput is the 
simulation time step. Their contribution to the radia- 
tion is computed similarly to the test-particle code [i.e., 
using Eq. (1)]. We note, however, that to reduce the 
computing time, it is possible to use the synchrotron 
radiation formulas instead of Eq. (1), since betatron 
radiation occurs in the wiggler regime. For each elec- 
tron, the angle-integrated spectrum can be calculated by 
dl/duj = J dP/duj{t)dt, where dP/d 10 is the synchrotron 
spectrum given in Eq. (8) and depends on the radius of 
curvature which can be evaluated numerically along the 
arbitrary electron trajectory. The angle-dependent spec- 
trum d 2 1 1 dcudfl can also be obtained in the wiggler limit 
by using the "saddle point" method (Kostyukov et ai, 
2003) . The saddle points are defined as the points in the 
electron trajectories where the velocity vector v is paral- 
lel to a given observation direction n. The total radiation 
spectrum is then obtained as the sum of the synchrotron- 
like bursts of all saddle points. 

A 30 fs (FWHM) laser pulse is focused at the begin- 
ning of the density plateau of a 4.2 mm wide gas jet, to 
a peak normalized electric field of do = 1.2. It is linearly 
polarized in the e x direction and the focal spot size is 18 
/j,m (FWHM). The gas jet density profile consists in a 
density plateau of 3 mm at n e = 1.5 x 10 19 cm -3 and the 
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FIG. 10 Electronic density in the z — y plane averaged over 
5 |im in the x direction. Densities are given in cm" 3 , z is the 
distance from the beginning of the first density ramp. Differ- 
ent steps of the simulation are shown. First the bubble forms, 
then injection starts, a whole bunch of accelerated electrons 
is created inside the bubble and finally the bunch comes out 
of the bubble and is scattered. The dynamics of the bunch is 
such that the trailing electrons have a much higher oscillation 
amplitude than leading electrons. 

edges are described by ramps varying linearly from to 
n e = 1.5 x 10 19 cm -3 over 600 fim. The simulation de- 
scribes all the successive steps of the bubble regime (see 
Fig. 10): penetration of the laser pulse in the plasma, 
formation of the bubble, injection of electrons at the back 
of the bubble and acceleration, wiggling of the electron 
bunch in the bubble, laser energy depletion and scatter- 
ing of the bunch in the plasma when it comes out of the 
focusing transverse electric field of the bubble. 

Figure 11 shows the electron beam spectrum before 
scattering. Most of the electrons are around 100 MeV 
but the tail extends up to 350 MeV. The total charge of 
the bunch is 640 pC, i.e., approximately 4x 10 9 electrons. 

Figure 12 focuses on a single particle. It illustrates 
the fact that a particle in the bunch, first accelerated 
and then dephased, follows a sinusoidal orbit around the 
propagation axis and radiates the most when it has both 
a high 7 and a strong perpendicular acceleration (when 
the curvature radius is small). 

In a more global point of view, Fig. 13 shows typical 
trajectories of electrons during their acceleration inside 
the bubble in the polarization plane. Each line repre- 
sents the trajectory of one macroparticle in the code. 
The color of a line changes as the macroparticle is ac- 
celerated to higher energy (represented here by its rela- 
tivistic 7 factor). The 500 trajectories plotted in Fig. 13 
are computed from macroparticles injected around the 
same instant and thus in the same conditions. As the 
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FIG. 11 Normalized electron beam spectrum before scatter- 
ing at 2 — 1800 /im. 
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FIG. 12 Trajectory, energy and radiated power of a single 
electron before scattering. The tube shows the particle's tra- 
jectory. Its radius is proportional to the particle's energy and 
the color measures the instantaneous radiated power of the 
particle. The orbit of this particle is in the polarization plan. 

behavior of a particle strongly depends on its injection 
conditions, it is not surprising to see that, in majority, 
the represented particles follow very close orbits during 
the acceleration phase. Moreover, the oscillation is mod- 
ulated and slightly amplified by the interaction with the 
laser pulse. This interaction was neglected in the test- 
particle simulation. 



20 




FIG. 13 Projection in the z~x plane of trajectories and ener- 
gies of 500 particles from injection until scattering. The laser 
is polarized in the e x direction and modulates the particle 



trajectories. Note the different scales in x and z axes. FIG. 15 Angular profile of the emitted radiation. 
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FIG. 14 Trajectories and radiated power of the same 500 
particles as Fig. 13. The laser is polarized in the e x direction. 




Figure 14 shows trajectories and radiated power of the 
same particles. It can be directly compared to Fig. 13. 
As expected, most of the radiation is produced inside the 
bubble at high electron energies and far from the axis, 
where the transverse acceleration is the highest. 

As shown in Sec. II, the radiation is emitted in the 
direction of the electron velocity. Figure 15 shows the 
angular distribution of the radiated energy. It is very 
well centered on the propagation axis and has an angular 
width of 30 mrad (FWHM). 

The frequency distribution of the radiated energy in- 
tegrated over all angles also yields interesting informa- 
tion. Figure 16 represents the calculated x-ray spectrum 
emitted by all injected macroparticles. The spectral dis- 
tribution is given by the number of emitted photons per 
0.1% bandwidth. The total number of emitted photons 
is N 7 ~ 4.5 x 10 9 and we have 10 6 photons/0. 1%BW 
at 1 kcV and 3 x 10 5 photons/0. 1%BW at 10 keV. This 
simulation reproduces the global behavior observed in the 
experiments (see Sec. IV. D) but tends to produce more 
photons with a spectrum extending to higher energies. 
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FIG. 16 Spectrum of the x-ray radiation emitted by all in- 
jected macroparticles. 



D. Experimental results 

The production of betatron radiation requires a rela- 
tivistic electron bunch and an ion cavity or channel. The 
laser-plasma interaction is not the only mechanism that 
can create this ion cavity. Another approach was inves- 
tigated first experimentally: an electron beam propagat- 
ing in a plasma can expel background electrons via the 
repulsing electric space-charge force and create an ion 
cavity if its own density is larger than the plasma den- 
sity. The cavity being created by the front edge of the 
bunch, its main body propagates in the self-produced 
cavity and experiences its fields. This scheme of beta- 
tron x-ray emission was demonstrated by Wang et al. 
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FIG. 17 Betatron x-ray beam (round spot), produced by the 
28.5 GeV SLAC electron beam propagating through a plasma 
of density ~ 2 x 10 13 cm -3 , and recorded by a fluorescent 
screen imaged onto a visible CCD camera. From Wang et al. 
(2002). 



(2002). The experiment, carried out at the SLAC Na- 
tional Accelerator Laboratory, was based on the use of 
the 28.5 GeV SLAC electron beam, containing 1.8 x 10 10 
electrons per bunch and focused near the entrance of a 
preformed lithium plasma of length 1.4 m. For a plasma 
density of 1.7 x 10 14 cm" 3 , the strength parameter of 
betatron oscillations was K = 16.8, and the number of 
emitted photons and the peak spectral brightness were 
estimated to be respectively 6 x 10 5 photons/0. 1%BW 
and 7 x 10 18 photons/ (s mrad 2 mm 2 0.1%BW), at 14.2 
keV. The betatron x-ray beam (shown in Fig. 17 for a 
smaller plasma density, ~ 2 x 10 13 cm" 3 ) was found to 
be collimated in the divergence angle of ~ (1 — 3) x 10~ 4 
rad. The expected quadratic density dependence of the 



spontaneously emitted betatron x-ray radiation [P 7 oc 
7 2 -ff 2 /A 2 cx n 2 at constant 7 and 773; see Eqs. (13) and 
(49)-(50)] was also observed, with a detection sensible in 
the spectral range 5-30 keV. These features were found 
to be in good agreement with theoretical predictions. 




FIG. 18 Experimental setup for the generation and observa- 
tion of betatron radiation. A 50 TW / 30 fs class laser is 
focused onto the front edge of a gas jet. Electron spectra are 
obtained by deflecting electrons with a magnet and using a 
LANEX phosphor screen imaged onto a CCD. Betatron ra- 
diation is recorded on axis using a charged coupled device 
sensitive up to 10 keV (CCD X). A filter (Be, Al, etc.) is 
placed before the CCD to block infrared light. 
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FIG. 19 Left: A raw image of the angular distribution of the 
betatron radiation beyond 3 keV recorded by the CCD in a 
single shot. Right, divergence (full width at half maximum) 
of radiation as a function of density. The dotted line shows 
the dependence of 9 on n e for constant rp and 7 according to 
the equations of Sees. IV. A and Sec. IV. B. 



The first experiment based on a fully laser-plasma ap- 
proach discussed above was performed by Rousse et al. 
(2004) using a 50 TW / 30 fs laser focused onto a 3 mm 
helium gas jet at an electron density of ~ 10 19 cm" 3 . 
A typical betatron experimental setup is shown in Fig. 
18. An intense femtosecond laser pulse is focused onto a 
supersonic helium gas jet. Electrons accelerated during 
the interaction are deflected toward a phosphor screen 
using a permanent magnet with a field approaching 1 
T. Because the magnet bends the electron trajectories in 
only one direction, the image on the screen reveals the 
electron spectrum and the beam divergence. The x-ray 
radiation is observed using an x-ray CCD placed on axis. 



With the camera sensitive in a broad spectral range, all 
radiation from the laser and the plasma below the keV 
range is blocked using filters (beryllium and aluminum). 

The left part of Fig. 19 shows the x-ray beam profile 
observed on the camera (at energy > 3 keV as a 500 /im 
beryllium filter was placed in front of the camera). The 
right part of the figure represents the evolution of the 
beam divergence (full width at half maximum) as a func- 
tion of the electron density. The radiation is collimated 
on axis with a divergence in the range [10 : 100] mrad. 
The divergence increases with the electron density due 
to the reductions of both the betatron period and the 
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FIG. 20 X-ray signal (proportional to radiated energy beyond 
1 keV per unit solid angle) as a function of plasma electron 
density. Each point corresponds to an average value over ten 
shots. The dotted line corresponds to the results obtained 
using 3D PIC simulations. From Rousse et al. (2004). 
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FIG. 21 Raw electron spectra for three different densities and 
fixed laser parameters. Horizontal axis, electron energy; ver- 
tical axis, exit angle; and color scale, number of counts. This 
latter gives an indication of the beam charge. No electron is 
observed below the trapping threshold (n e = 7.5x 10 18 cm -3 ), 
a quasimonoenergetic electron beam is produced just above 
(n e — 1.25 x 10 19 cm -3 ), and a broadband spectrum is 
recorded for a higher density (n e = 1.5 x 10 19 cm -3 ). The 
threshold density is not the same as in Fig. 20 because the 
laser parameters are different. 



electron energy as observed in the experiment. These 
experimental observations are in good agreement with 
the numerical predictions presented above (Albert et al, 
2008; Ta Phuoc et al, 2005a, 2008a, 2006) and with PIC 
simulations (Kiselev et al, 2004; Pukhov et al, 2004b; 
Rousse et al, 2004; Ta Phuoc et al, 2005a). 

The betatron radiation has a specific dependence with 
plasma density (for given laser parameters) (Rousse 
et al, 2004; Ta Phuoc et al, 2005a). Figure 20 shows 
the maximum of the radiated energy per unit solid angle 
beyond 1 keV (in arbitrary units), as a function of the 
electron density, obtained in the experiment and in PIC 
simulations. The observed x-ray emission is peaked at 
n e ~ 1.1 x 10 19 cm -3 . Below 5 x 10 18 cm -3 , the x-ray sig- 
nal vanishes simply because no electron is trapped (the 
laser intensity is too low for such densities in order to en- 
sure electron trapping). This is confirmed in the experi- 
ment for which no electron was detected by the spectrom- 
eter. Just above this threshold, electrons are trapped in 
the bubble regime and monoenergetic electron bunches 
as well as betatron radiation of relatively low intensity 
can be obtained. When the density is higher, the plasma 
wavelength is smaller and the focal spot size and pulse 
duration are not exactly matched. It first corresponds to 
the FLWF regime (Malka et al, 2002; Najmudin et al, 
2003), and then for even higher density to SM-LWFA 
(Coverdale et al, 1995; Modena et al, 1995; Najmudin 
et al, 2003; Santala et al, 2001; Ting et al, 1997; Um- 
stadter et al, 1996a; Wagner et al, 1997). This results 
in acceleration of electron bunches of poor quality com- 
pared to the bubble regime (broad spectrum, high diver- 
gence, and large shot-to-shot fluctuations). The coun- 
terpart of this poor quality electron beam is that the 



number of electrons and the amplitude of motion rp are 
higher and the betatron period is smaller, such that the 
betatron radiation is brighter even if the mean electron 
energy is smaller. An additional effect which can occur 
is the interaction of the electron bunch with the laser 
pulse. From n e ~ 1.1 x 10 19 to 2.7 x 10 19 cm -3 , the 
x-ray signal drops down and a plateau is reached. The 
PIC numerical simulations (Rousse et al, 2004; Ta Phuoc 
et al, 2005a) shown in Fig. 20 reproduce this experi- 
mental behavior: a sharp increase of the x-ray intensity 
followed by a smoother decrease of the signal. Three 
typical electron spectra which can be obtained simulta- 
neously with the measurement of betatron radiation are 
represented in Fig. 21. The first shows a broadband spec- 
trum obtained at relatively high density, characteristic of 
the forced laser wakefield regime, and the second corre- 
sponds to a monoenergetic electron bunch (signature of 
the bubble regime) obtained for a density just above the 
trapping threshold. Below this threshold, no electron is 
observed (third spectrum). For some broadband spectra, 
a correlation of the electron output angle with the elec- 
tron energy can be observed. It is attributed to off-axis 
injection of the electron beam due to asymmetric laser 
pulse intensity profile or tilted energy front. For these 
shots, direct measurement of the betatron motion can be 
obtained from the oscillation observed in the spectrum 
(Glinec et al, 2008). 

The spectrum of the betatron radiation was first mea- 
sured using a set of filters (Rousse et al, 2004; Ta Phuoc 
et al, 2005a), for an electron density n e = 1 X 10 19 cm -3 . 
The top panel of Fig. 22 presents the number of x- 
ray photons obtained experimentally, integrated over the 
beam divergence and over the spectral bandwidths deter- 
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FIG. 22 Top: Total number of photons obtained experimen- 
tally within the spectral bandwidths determined by filters: 25 
/jmBe (1<£< 10 keV), 25 fim Be + 40 fim Al (4 < E < 10 
keV), and 25 /im Be + 25 )im Cu (6 < E < 10 keV). From 
Rousse et al. (2004). Bottom: Averaged spectrum (in red 
line) of the radiation obtained by photon counting, with a fit 
(in green line) by the function S{uj/uj c ) with fno c = 5.6 ± 1 
keV. From Fourmaux et al. (2011b). 



mined by (1) 25 /imBe (1< E < 10 keV), (2) 25 Be 
+ 40 nm Al (4 < E < 10 keV), and (3) 25 Be + 25 
^m Cu (6 < E < 10 keV). The bottom part of the figure 
represents a measurement of the spectrum performed by 
photon counting, at higher resolution (350 eV), in the 
range 8 to 21 keV. In this experiment, betatron radia- 
tion was produced with a 80 TW laser system (2.5 J and 
30 fs) at the Advanced Laser Light Source (Fourmaux 
et al., 2011b), for an electron density n e — 5.4 x 10 18 
cm~ 3 . The experimental spectrum (averaged over 10 
shots) was well described by a synchrotron function at 
hu) c = 5.6 ± 1 keV in the range 8 to 21 keV, demonstrat- 
ing the synchrotron-type spectrum of betatron radiation 
(Fourmaux et al., 2011b). Depending on the experiment 
and available laser energy in the focal spot, the total 
number of photons range from 10 8 (Rousse et al, 2004) 
to 10 9 with 2.2 x 10 8 photons/(sr 0.1%BW) at 10 keV 
(Fourmaux et al., 2011a). The latter result also showed 
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FIG. 23 Synchrotron spectrum fitted from transmission val- 
ues through six different filters. The fit uses the synchrotron 
function (uj/lo c ) 2 Ky 3 (tu/uj c ) instead of S(w/ui c ), and the fit- 
ted critical energy is 29 ± 13 keV. From Kneip et al. (2010). 



a critical energy above 10 keV (Fourmaux et al., 2011a). 

Betatron radiation has also been observed in a recent 
experiment performed on the Michigan Hercules laser 
system at the University of Michigan (Kneip et al., 2010). 
The results demonstrated high photon energies and small 
x-ray beam divergence. The betatron radiation had a 
broadband spectrum, characterized using a set of filters. 
For an electron density n e = 1 x 10 19 cm~ 3 , the trans- 
mission values through the six filters were fitted with the 
synchrotron function (Lo/u! c ) 2 Ky 3 (u>/uj c ) [which differs 
from S(oj/uj c )] with a critical energy Huj c ~ 29 ± 13 keV 
(see Fig. 23). The divergence angle was found to be 
~ 10 mrad (FWHM), and the electron spectrum, simul- 
taneously recorded, indicated electron acceleration up to 
400 MeV. The peak spectral brightness was estimated 
to be ~ 10 22 photons/(s mrad 2 mm 2 0.1%BW), which 
is comparable to currently existing third-generation con- 
ventional light sources. 

The duration of short-pulse betatron x rays was first 
estimated using time resolved x-ray diffraction. In this 
experiment (Ta Phuoc et al., 2007), an ultrafast phase 
transition (nonthermal melting of InSb) was used as a 
Bragg switch to sample the x-ray pulse duration. The 
dynamics of this specific transition is well known. The 
InSb is initially crystalline and diffracts the x-ray radia- 
tion according to the Bragg law. However, when irradi- 
ated by a femtosecond laser pulse (at a fluence in the 100 
mJ/cm 2 range), the surface of the InSb is melted and dis- 
orders appears in a few tens of femtoseconds. As a conse- 
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quence, x-ray radiation can no longer be diffracted. The 
duration of this phase transition has been determined 
in several experiments based on the use of femtosecond 
synchrotron radiation (Lindcnberg et al., 2005) and fem- 
tosecond K a radiation (Roussc et al., 2001a). Here, as- 
suming the dynamics of the phase transition is known, 
the duration of the betatron x-ray pulse was estimated 
to be less than 1 ps with a best fit below 100 fs (Ta Phuoc 
et al., 2007). Recently, the electron bunch duration was 
measured (Lundh et al., 2011) using spectral characteri- 
zation of coherent transition radiation (CTR), in a laser- 
plasma accelerator working with external colliding pulse 
injection (Faure et at, 2006). The CTR result demon- 
strates a root-mean-square duration of 1.5 fs for the elec- 
tron bunch. Betatron radiation was shown to be strongly 
correlated with the properties of these electron bunches 
(Corde et al., 2011a), showing that the radiation is effec- 
tively produced by these ultrashort electron bunches and 
that it inherits its duration of a few femtoseconds. 

The source size has been determined with two differ- 
ent methods. The first uses the Fresnel edge diffraction 
of the x-ray beam and is based on the spatial coherence 
of the radiation (Shah et al., 2006). Because electrons 
emit x rays incoherently (see Sec. II. F) in the betatron 
mechanism, the source size determines the degree of spa- 
tial coherence which results in interference fringes in the 
Fresnel diffraction experiment. This latter experiment 
corresponds to the measurement of the intensity profile 
of the shadow made by a knife edge (here the edge of a 
cleaved crystal) placed in the x-ray beam onto a detector 
at a distance sufficient to ensure an appropriate resolu- 
tion. Using this technique, it was demonstrated that the 
source size was less than 8 /im (FWHM). A knife-edge 
technique was also used in betatron experiments with 
the Michigan Hercules laser system at the University of 
Michigan (Kncip et al., 2010), and with the Advanced 
Laser Light Source laser system (Fourmaux et al., 2011a), 
showing an x-ray source size of ~ 1 — 2 /im (FWHM). The 
second method relics on the information contained in the 
spatial distribution of the radiation. Indeed, because the 
betatron radiation is emitted in the direction of the elec- 
tron velocity, its profile represents a direct signature of 
the electron orbits in the cavity. In addition, the diver- 
gence and the structure of the distribution depends on 
the amplitude of the electrons orbits. This method, dis- 
cussed in the previous section, is presented in references 
Ta Phuoc et al. (2008a, 2006). It is shown that the spa- 
tial distribution of the radiation implies a limited choice 
of electron orbits. Therefore, by determining the electron 
orbits, the source size was deduced and the result was in 
the range of ~ 1 - 2 /xm (FWHM). 

Betatron x-ray radiation has interesting properties for 
application experiments: high peak spectral brightness, 
ultrashort duration, very small source size, as well as 
femtosecond time-scale synchronization in pump-probe 
experiments. The potential of betatron radiation was 
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FIG. 24 Bee imaged with the betatron x-ray beam with an 
edge line out indicated by the rectangular area: (a) one laser 
shot; (b) 13 laser shots, (c) Line out of the images taken at the 
indicated rectangular area. From Fourmaux et al. (2011a). 

demonstrated with the example of single shot phase con- 
trast imaging of biological samples (Fourmaux et al., 
2011a; Kneip et al., 2011), where the experiment takes 
advantage of the very small source size and high spatial 
coherence of the x-ray beam. For an x-ray source size of 
^1 — 2 /im, the coherence length at 1 m from the source is 
~ 10 /xm, which allows to perform phase contrast imaging 
with a compact setup. Figure 24 shows a phase contrast 
of a bee in single shot [Fig. 24(a)] as well as with accu- 
mulation over 13 shots [Fig. 24(b)] and the contrast is 
found to be 68 % in single shot for the indicated line out 
(Fourmaux et al., 2011a). Betatron x-ray radiation is also 
a powerful tool to study laser-plasma accelerator physics 
and for noninvasive measurement of these properties. A 
novel method has shown the possibility to measure the 
longitudinal profile of the x-ray emission region, giving 
insight into the history of the laser-plasma interaction 
(Corde et al., 2011b). Betatron radiation from optically 
injected, tunable, and monoenergetic electron bunches 
have demonstrated the strong correlation between elec- 
tron and x-ray properties (Corde et al., 2011a), showing 
that betatron radiation can be used to perform noninva- 
sive measurements of some of the electron beam param- 
eters, such as the normalized transverse emittance of the 
electron beam (Kneip et al., 2012). 

Betatron radiation was also observed in the petawatt 
regime using a laser delivering 300 J within 600 fs pulses 
(Kneip et al., 2008). In this parameter regime, the laser 
pulse is much longer than the plasma period and the 
physical mechanisms are different. For a long laser pulse 
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(tl ^ 1/wp), electrons can be accelerated and wiggled 
in the SM-LWFA and/or DLA regimes, depending on the 
laser strength parameter clq. Electrons experience both 
the transverse focusing force of the channel and the elec- 
tromagnetic field of the laser. As a consequence, the 
transverse amplitude of the electron orbits can be signif- 
icantly increased due to betatron resonance. The results 
showed the possibility to produce betatron radiation with 
a spectrum fitted by the synchrotron function S(uj/uj c ) 
with a critical energy up to hw c — 14.5 keV [note that 
in Kneip et al. (2008), their convention for Huj c leads to 
a value of 29 keV] and with a divergence of ~ 1 rad. In 
this experiment, the total number of photons was found 
to be more than 5 x 10 8 photons/ (mrad 2 0.1%BW) in 
the spectral range [7 keV; 12 keV] and a typical electron 
charge of 2 nC was observed with electron energies in the 
50 MeV range. From these properties, the amplitude of 
motion and the strength parameter K were deduced to 
be ~ 30 /im and ~ 130 respectively, in agreement with 
PIC simulations. 



E. Scalings and perspectives 

The betatron source offers several routes of develop- 
ment in terms of number of photons, spectral range and 
divergence. According to Sec. IV. B, these quantities 
scale in the wiggler limit as (Rousse et at, 2007) 

A 7 oc N e N K oc N e Np^fT e rp, (64) 
fiw c a 7 L5 ^if oc 7 2 n e r,3, (65) 
B r oc K/sy oc y^rp/y/j, (66) 

where Np is the number of betatron periods and N e is 
the number of electrons. As discussed in Sec. IV. B, these 
formulas are valid for time-constant parameters, but can 
be used with the values of the parameters at the de- 
phasing time for the realistic case of an electron being 
accelerated in the ion cavity. In addition to these formu- 
las, the scalings of the acceleration process are required. 
Considering the phenomenological model of the bubble 
regime described by Lu et al. (2006a, b, 2007), higher elec- 
tron energy gain requires lower plasma density. Its scal- 
ing is given by 7 ~ A£/mc 2 — (2w| / /3Wp)ao oc ao/n e . 
The matched spot size condition reads k p wo = 2y/a~o 
and the bubble radius should equal the laser pulse waist 
7*6 = wo oc a^/ 2 n e 1 ^ 2 . The number of trapped electrons, 
obtained from an energy balance between the electron 
beam and the bubble electromagnetic fields, scales as the 
ionic charge of the bubble: N e oc rln e oc a^ 2 n e 1 . The 
dephasing length scales as Ld oc a^rie and the beta- 
tron period as A u oc a\ n~ x , such that the number of be- 

— 1/2 

tatron oscillations obeys Np oc n e . Different electron 
injection mechanisms can lead to different scalings for rp, 
but betatron radiation is optimized when rp is maximum, 



and this latter is limited by rb- We therefore assume that 
the transverse size of the electron bunch roughly scales as 

the bubble radius for optimized conditions for betatron 

,. 1/2 -1/2 

radiation: rp oc rj, oc a n e 

The scalings of Eqs. (64), (65) and (66) can then be 

reformulated as functions of ao and n e as 

A 7 oc a 5 /2 n~ 3/2 , (67) 
hu c oc a 5 Q /2 n- 3/2 , (68) 
B r oc n\ /2 , (69) 

while the needed laser energy .Eiaser oc a^w^T obeys 

£aascr <x a 7 /2 n-^ 2 . (70) 

This expression uses the optimal pulse duration cr = 
2wq/3 corresponding to the situation where the dephas- 
ing length Ld matches the pump depletion length L p d- 
These scalings take into account neither the possible in- 
teraction between the electron beam and the laser pulse 
nor beam loading effects, and describe the betatron per- 
spectives in the purely bubble regime. However, they 
highlight the route toward harder x rays and higher radi- 
ated energies: increasing the laser pulse strength param- 
eter and/or decreasing the density. Equation (70) indi- 
cates that it is more efficient in terms of laser energy to 
decrease the plasma density and not to increase the laser 
strength parameter. The efficiency r\x of the conversion 
of laser energy into x rays, defined as the ratio of the radi- 
ated energy by the laser energy, scales as <Zo rig 3 . The 
betatron mechanism will become more and more efficient 
when going to higher and higher photon energies. 

Petawatt class lasers and capillaries both allow one to 
decrease the plasma density. With petawatt lasers pro- 
ducing sub-100 fs pulses, the higher laser energy can be 
focused into a larger focal spot matched to the lower 
plasma density according to k p wo = For self- 

injection to occur, ao ~ 4 — 5 is needed, while the 
self-guided propagation requires a > {ric/n,,) 1 ^ , where 
n c = meoui^/e 2 is the critical density (this condition ex- 
presses that energy loss at the front of the laser pulse due 
to pump depletion is higher than the loss due to diffrac- 
tion) (Lu et al. , 2007) . Capillaries allow one to guide the 
laser pulse over distances much larger than the Rayleigh 
length such that no condition apply on ao- This tech- 
nique has already been used to accelerate electrons up 
to 1 GeV in a few centimeters capillary (Lecmans et al., 
2006). Combined with external injection, a strength pa- 
rameter as small as ao = 2 can be used, which permits 
one again to focus the laser energy into a larger focal spot 
and to decrease the plasma density. A drawback of these 
natural scalings is the increase of the x-ray source size 
and pulse duration as rp,r r oc a^ 2 n e 1 ^ 2 . As a conse- 
quence, the peak spectral brightness (number of photons 
divided by the pulse duration, by the solid angle, and 
by the source size in 0.1% spectral bandwidth) scales as 
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ao/n e . Nevertheless, an external injection could provide 
small duration and therefore better time resolution and 
brightness for application experiments. 

Using 15 J of laser energy focused into a gas jet of den- 
sity n e — 1.1 x 10 18 cm -3 with a waist wq = 21 /xm and 
compressed to a duration r = 48 fs (300TW on target), 
the laser strength parameter is ao = 4.4, ensuring self- 
guiding and self-injection, and trapped electrons reach 
the maximum energy of 2.4 GeV after 23 mm of plasma 
(Rechatin, 2009). Assuming a reasonable value for the 
betatron amplitude rp = 3 fim, we expect K ~ 30 and 
x-ray photons emitted up to Huj c ~ 400 keV, collimated 
within a typical opening angle of 6 mrad, and containing 
jV 7 ~ 1 photon per electron and per betatron period. On 
the other hand, a configuration with 3 J of laser energy, 
n e = 5.1 x 10 17 cm" 3 , w = 21 fxm, r = 47 fs (60TW on 
target) corresponding to a Q = 2 and assisted by exter- 
nal guiding (capillaries) and external injection [colliding 
pulse injection (Faure et al., 2006) or density gradient 
injection (Geddes et at, 2008)] provides the same out- 
put electron energy of 2.4 GeV after 52 mm of plasma. 
With the same value of betatron amplitude rp = 3/Ltm, x 
rays are emitted with the following properties: K ~ 20, 
huj c ~ 200 keV, 6 r ~ 4 mrad, and iV 7 ^0.7 photon 
per electron and per betatron period. Compared to the 
self-guided self-injected situation, the number of trapped 
electrons N e is 2 times smaller, the number of betatron 
periods Np is 1.4 times higher (see their respective scal- 
ings above), and the total number of emitted photons 
is 2 times smaller, but within a better divergence (6 r is 
1.4 times smaller) and for a needed laser power 5 times 
smaller. For such configurations, the number of betatron 
photons could reasonably reach the level of 10 10 photons. 



Another route toward brighter betatron radiation and 
higher photon energies makes use of density-tailored plas- 
mas (Layer et al., 2007) as recently suggested (Ta Phuoc 
et al., 2008b). This method relies on the control of 
the electron orbits by varying the density and therefore 
the forces acting on the electrons along the propagation. 
For appropriately chosen density modulations, numeri- 
cal simulations show that the amplitude of the betatron 
oscillations rp can be significantly increased. The basic 
idea is to freeze the wakefield during a short duration 
in order to let the electron beam defocus and acquire 
larger betatron amplitude. An additional effect is that, 
for large betatron amplitude, the mean longitudinal ve- 
locity is lowered and electrons dephase slower and can 
reach higher energies. In an optimistic configuration de- 
scribed by Ta Phuoc et al. (2008b), the critical energy 
of the radiation can be increased by a factor of « 9, the 
number of photons by a factor of ~ 5, but with an in- 
crease of the divergence by a factor sa 4. 
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FIG. 25 Schematic of the undulator source. A laser-plasma 
accelerated electron beam is injected in a conventional undu- 
lator. Electrons oscillate transversally and emit x rays due to 
this motion. 



V. PLASMA ACCELERATOR AND CONVENTIONAL 
UNDULATOR: SYNCHROTRON RADIATION 

Synchrotron radiation can be produced from laser- 
plasma accelerated electrons propagating and oscillating 
in a conventional undulator or wiggler. A conventional 
undulator or wiggler is a periodical structure of magnets 
generating a periodic static magnetic field (Clarke, 2004). 
Figure 25 represents the principle of the scheme. Elec- 
trons are first accelerated in a plasma accelerator (gas jet, 
capillary or steady-state-flow gas cell) and then trans- 
ported into an undulator or wiggler. The size of such 
a device is on the meter scale, with a typical magnetic 
period in the centimeter range. Moreover, the values of 
X u , technologically limited, imply that GeV electrons are 
necessary to produce radiation in the x-ray range. These 
represent the disadvantages of this scheme, which how- 
ever is of crucial interest because it represents a possible 
route toward the production of a compact free-electron 
laser based on a laser-plasma accelerator (Griiner et al., 
2007; Malka et al, 2008; Nakajima, 2008) (see Sec. VII). 
In the following sections, the electron orbits in undulators 
and the main properties of synchrotron radiation pro- 
duced from laser-plasma electron bunches are described 
in a presently realizable case. After a summary of the ex- 
perimental results, we conclude with the short term de- 
velopments foreseen. For detailed properties of undulator 
radiation, such as the angular distribution of individual 
undulator harmonics, the different undulator types (pla- 
nar or helical) , the polarization properties and the spatial 
or temporal coherence properties, see Chao and Tigncr 
(2009); Clarke (2004); and Wiedemann (2007a). 



A. Electron motion 

We consider a laser wakefield accelerator producing a 
monoenergetic electron bunch of energy £ — jimc 2 with 
velocity directed in the e z axis. The simplest model to 
describe the electron dynamics in the undulator is to ap- 
proximate the static magnetic field near the undulator 
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axis by 



B = B sm(k u z)e y , 



(71) 



where k u = 2tt/X u with A u the magnetic period of the 
undulator. The equation of motion for a test electron in 
this idealized model is given by 



dp 



= —ev x B. 



(72) 



The normalized Hamiltonian describing the electron dy- 
namics is given by 



(73) 



where P is the normalized canonical electron momentum 
and a = —ao cos(k u z)e x is the normalized vector poten- 
tial verifying B = (mc/e)V X a. H does not depend on x, 

y, and t. The transverse canonical momentum P± and 
the normalized energy 7 are therefore conserved. Inte- 
grating these constants of motion with the initial condi- 
tions P± = and 7 = 7i leads to, in physical units, 



7 = etc = 7i, 
K 

x(z) ~ — sin(k u z), 

7«u 



K = a 



eB 
k u mc 



(74) 
(75) 

(76) 



The general description developed in Sec. II is recovered 
and the practical expression of K in terms of experimen- 
tal parameters and fundamental constants has been ob- 
tained. In particular, K does not depend on 7 and the 
radiation regime (undulator or wiggler) is the same for 
all electron energies. The strength parameter K is deter- 
mined by the amplitude of the magnetic field Bq and the 
period A u and is, in practical units, given by 



K = 0.934A„[cm] J B [T] 



(77) 



B. Radiation properties 



The radiation properties are the same as in Sec. II. 
The spectrum of the emitted radiation depends on K. 
For a small amplitude of oscillation K -C 1, radiation 
is emitted at the fundamental photon energy hw with a 
narrow bandwidth in the forward direction (8 = 0). As 
K — > 1 harmonics of the fundamental start to appear 
in the spectrum, and for if > 1 the spectrum contains 
many closely spaced harmonics and extends up to a crit- 
ical energy Huj c . These quantities are given by 



Hlo = (2 7 2 /ic/A„)/(l + K 2 /2) forif<l, 
3 

hjj c = -K~/ 2 hc/\ u for K > 1. 



(78) 



In practical units, this gives 

hu[eV] = 2.48 x KTV/^Ucm] f or k < 1, 
twj c [eV] = 1.74 x l(T 4 7 2 £o[T] for K > 1. (79) 

The radiation is collimated within a cone of typical 
opening angle 8 r = 1/7 in the undulator case. For a wig- 
gler, the radiation is collimated within a typical opening 
angle if/7 in the electron motion plane (e x , e z ) and I/7 
in the orthogonal plane [e v , e z ). 

The radiation temporal profile is simply given by the 
convolution between the electron bunch temporal pro- 
file and the radiation profile from a single electron. For x 
rays, the radiation length from a single electron l r \ N =1 = 
iVA is in the nanometer range, much smaller than the typ- 
ical electron bunch length which is in the micron range, 



and therefore r. 



Tb- However, for UV and XUV ra- 



diation, the slippage between the radiation and the elec- 
tron bunch is not necessarily negligible, and the radiation 
duration can be stretched by the single electron radia- 
tion duration. The electron bunch duration depends on 
the acceleration mechanism and on its transport from its 
source to the undulator. Presently, the shortest electron 
bunches are obtained using a laser-driven plasma-based 
accelerator with an external optical injection mechanism 
(Davoine et ai, 2008; Faure et al, 2007, 2006; Malka 
et ai, 2009; Rechatin et at, 2009b). In this scheme, elec- 
tron bunches with rms (root-mean-square) duration as 
small as 1.5 fs can be produced (Lundh et ai, 2011). 
However, the necessary transport of the electron bunch 
in the undulator can deteriorate the duration (and the 
transverse emittance), because different electrons in the 
bunch travel different distances. For example, consider- 
ing usual transport devices (quadrupoles), placed at ~ 1 
m from the laser-plasma accelerator, an initially 1-^m- 
long bunch has its length increased to ~ 10 /im during 
its transport to an undulator placed several meters after 
the laser-plasma accelerator. An ultracompact transport 
system, placed as closed as possible to the plasma source, 
is required to avoid such a deterioration. 

The number of emitted photons follows the expressions 
given in Sec. II. In practical units, the number of photons 
emitted per period and per electron (at the mean photon 
energy (huj) — 0.3fiiv c for the wiggler limit) is given by 



N-y = 1.53 x IQ- 2 K 2 for K < 1, 
A 7 = 3.31 x 1Q- 2 K foriT>l. 



(80) 



From these expressions the features of the produced ra- 
diation can be obtained for a realizable case: an electron 
bunch at an energy of 1 GeV containing 100 pC and an 
undulator of 100 periods with A„ = 1 cm and K = 0.5. 
The emitted radiation has a narrow band spectrum cen- 
tered at Ex ~ 1 keV in the forward direction (6 = 0), is 
collimated within a typical opening angle of ~ 500 /zrad, 
and contains ~2x 10 s photons. 
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FIG. 26 Electron trajectories. The orbit is plotted as a red 
solid line for the undulator case (K = 0.2) and green dashed 
line for the wiggler case (if = 2). At the entrance, the electron 
has 7 = 1000 and p = \/"f 2 — le z in both situations. The 
undulator or wiggler is ten periods long. 



C. Numerical results 

In this section, the electron motion and the emitted 
radiation are calculated numerically. Both the cases of 
an undulator and a wiggler are considered, with strength 
parameter of K = 0.2 and 2, respectively. The magnetic 
structure has a period X u = 1 cm and the test electron 
propagates in the e z direction with 7 = 1000 (~ 500 
MeV). Figure 26 shows the electron motion along N = 10 
oscillation periods (the number of periods of an undula- 
tor or wiggler is usually on the order of 100). The motion, 
with 7 constant, consists of a transverse oscillation at a 
period X u combined with a longitudinal drift. The am- 
plitude of the transverse motion is given by X u K/2tt^/ 
(here ~ 1.6 x K[fxm\). For K = 0.2 and 2 the transverse 
excursions are 0.32 and 3.2 fjxa, respectively. 

The spatial distributions of the radiation produced by 
the electron undergoing these trajectories are presented 
on the left of Fig. 27. For K = 0.2 (top), the radia- 
tion has a divergence of typical opening angle 1/7 ~ 1 
mrad. For K = 2 (bottom), the radiation has a diver- 
gence of typical opening angle K/j ~ 2 mrad in the plane 
of the electron motion and 1/7 ~ 1 mrad in the orthog- 
onal direction. The right part of Fig. 27 represents the 
spectra of the emitted radiation for each situation (un- 
dulator on the top and wiggler on the bottom). For the 
undulator case, the spectrum is plotted for three different 
solid angles of integration (0.3 x 0.3, 0.8 x 0.8, and 2x2 
mrad 2 ), whereas for the wiggler case it is integrated over 
the overall spatial distribution. For K = 0.2 the spec- 
trum is nearly monochromatic at the energy hut ~ 200 
cV in the forward direction, with a width of the emission 
line given by AA = X/N (where N = 10 is the number of 
oscillation periods), and it is stretched after integration 
over the emission angles due to the angular dependence 
of the radiated wavelength A ~ (A u /27 2 )(1 + ^ + 7 2 6> 2 ) 
[see Sec. II, Eq. (5)]. For K = 2 the spectrum becomes 
broadband with a critical energy of around 300 eV. An 
estimate of the total photon number can be obtained on 
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FIG. 27 On the left is represented the angular distributions 
of the radiation (radiated energy per unit solid angle) and 
the corresponding spectra are shown on the right (in num- 
ber of photons per 0.1% bandwidth, per electron, and per 
shot). The undulator case, K = 0.2, is on the top whereas 
the wiggler case, K — 2, is on the bottom. In the first sit- 
uation, three different solid angles of integration have been 
considered to highlight the stretching of the spectrum. In the 
wiggler situation, the integration is performed over the total 
angular distribution. 

the basis of the number of electrons that are currently ac- 
celerated in a laser-plasma accelerator. Assuming 1 x 10 8 
electrons gives ~ 10 5 photons at 200 eV within a spectral 
bandwidth of 10% in the undulator case (selection of the 
0.3 x 0.3 mrad 2 solid angle). In the wiggler case there 
are ~ 5 x 10 3 photons/0. 1%BW at 600 eV and a total 
photon number in the range 10 7-8 . 



D. Experimental results 

Synchrotron facilities have been developed worldwide 
and have been providing x-ray radiation to users for sev- 
eral decades. These installations are very robust and the 
x-ray beams have high quality. They are based on con- 
ventional accelerator technology. 

The first demonstration of the production of syn- 
chrotron radiation from laser-accelerated electrons was 
performed by Schlenvoigt et al. (2008a, b). In this ex- 
periment, laser-produced electron bunches between 55 
and 75 MeV were injected, without transport, into a 
1-m-long undulator having a period X u — 2 cm and 
a strength parameter K = 0.6. They obtained syn- 
chrotron radiation in the visible and infrared part of the 
spectrum (the wavelength was in the range 700 — 1000 
nm) and estimated a peak spectral brightness of 6.5 x 
10 16 photons/(s mrad 2 mm 2 0.1%BW) and a total num- 
ber of photons of 2.8 x 10 5 . The radiation wavelength 
was observed to scale with the electron bunch energy as 
expected. 

More recently, synchrotron radiation from laser-plasma 
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accelerated electrons was produced up to ~ 130 eV 
(Fuchs et ai, 2009). This experiment represents an im- 
portant advance compared to the first demonstration. 
First, a stable electron beam was produced by the in- 
teraction of a laser at ~ 2 x 10 18 W.cm -2 with a 15-mm- 
long hydrogen-filled steady-state-flow gas cell (Osterhoff 
et at, 2008). The electron bunch was in the 150 - 220 
MeV energy range and contained 30 pC. Second, a pair 
of miniature permanent-magnet quadrupole lenses was 
implemented in order to make the electron beam colli- 
mated at a selected energy and to transport it through 
a 30-cm-long undulator having a period X u = 5 mm and 
a strength parameter K = 0.55 (Eichner et ai, 2007). 
Thanks to the high energy of the electron and the short 
period of the undulator, the radiation was produced down 
to the XUV range with a tunable wavelength. Note 
that this tunability resides in the possibility to select, 
with the magnetic lenses, a specific electron energy in 
the 150 - 220 MeV range (electrons with different en- 
ergies are not collimated and will not produce a col- 
limated radiation beam). The radiation was analyzed 
using a spectrometer based on a transmission grating. 
The fundamental wavelength (A = 18 nm) as well as 
the second harmonic were measured. In addition, the 
parabolic dependence of the radiation wavelength on the 
angular direction A = (A„/27 2 )(l + K 2 /2 + 7 2 # 2 ) was 
observed with good agreement. This soft-x-ray undu- 
lator source delivers ~ 10 5 photons per shot integrated 
over the fundamental within a detection cone of ±0.7 
mrad, corresponding to a relative bandwidth of 30% 
(FWHM). The peak spectral brightness was estimated 
to be - 1.3 x 10 17 photons/ (s mrad 2 mm 2 0.1%BW). 



E. Perspectives 

The experiments performed open perspectives for the 
production of a tunable synchrotron-type x-ray source. 
However, further progress is necessary to produce tun- 
able and reproducible monoenergetic electron bunches in 
the GeV range. The foreseen possibilities to achieve these 
features rely on the use of higher power lasers (petawatt) 
and/or capillaries. This gain in laser power compared to 
actual systems would allow one to decrease the plasma 
density, increase the acceleration length, and as a conse- 
quence increase the electron energy. The use of an exter- 
nal injection will be required as well to easily control the 
electron energy (optical injection (Faure et ai, 2006) or 
plasma density gradient injection (Geddcs et ai, 2008), 
as presented in Sec. III). Finally, a compact electron 
transport system from the plasma source to the undula- 
tor has to be implemented in order to avoid deterioration 
of the electron beam parameters and to reduce the radi- 
ation source size. 

In the short term, this approach will allow university- 
scale laboratories to access x-ray sources with angstrom 



wavelengths and sub-10-fs pulse durations for four- 
dimensional imaging with atomic resolution. Besides this 
future development, the main and more challenging per- 
spective relies on the possibility to produce a free-electron 
laser based on a laser-plasma accelerator combined with 
an undulator. This radiation source has a brightness or- 
ders of magnitude higher than a synchrotron, thanks to 
the coherent addition of radiation from each electron. 
Electron bunches produced with laser-plasma accelera- 
tors are interesting candidates for a FEL source since 
such bunches are naturally femtosecond and have high 
current. In Sec. VII, we discuss the principle of the free- 
electron laser and its possible realization using electrons 
from laser-plasma accelerators which could have revolu- 
tionary impacts in many fields of science, technology, and 
medicine. 



VI. ELECTROMAGNETIC WAVE UNDULATOR: 
NONLINEAR THOMSON SCATTERING AND 
THOMSON BACKSCATTERING 

In this section, the production of femtosecond x-ray 
and gamma-ray beams from electrons oscillating in an 
electromagnetic wave is reviewed. The involved radiative 
mechanism is referred to as Compton scattering. In the 
framework of quantum electrodynamics, it corresponds 
to the absorption of one (linear) or several (nonlinear) 
photons by an electron and to the emission of a single 
photon. When, in the rest frame, the electron experiences 
a negligible recoil (which happens when Hui <C mc 2 in the 
rest frame), the mechanism is referred to as Thomson 
scattering, the low-energy limit of Compton scattering. 
In the following, we consider Thomson scattering and 
use a classical description, assuming that quantum and 
radiation reaction effects are negligible (see Sec. II. G). 

At the laser-plasma interaction, two schemes of fem- 
tosecond x-ray or gamma-ray sources based on Thomson 
scattering have been proposed and demonstrated. Based 
on the same principle, the two methods differ by the ini- 
tial energy of the electron and the intensity of the laser 
pulse scattering off the electrons. 

In the first scheme, electrons are initially at rest and 
laser wakefield acceleration is not invoked. Here, Thom- 
son scattering occurs in a highly nonlinear regime. In- 
deed, for Bo > 1, electrons have a highly non linear 
motion and the emitted radiation consists of high-order 
harmonics. The harmonic spectrum can extend up to 
the x-ray range. This scheme offers simplicity and the 
possibility to produce a large x-ray photon flux because 
the number of electrons participating in the emission is 
on the order of the number of electrons in the focal vol- 
ume. However, reaching the keV energy range requires 
very high laser strength parameters, typically a,Q > 10. 

The second scheme relies on Thomson scattering an 
intense laser pulse off a counterpropagating relativistic 
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electron beam. Here the x-ray range is reached even 
with modest-energy electrons thanks to two successive 
Doppler shifts. In the average rest frame of the electrons, 
moving with a relativistic factor 7 with respect to the 
laboratory frame, the frequency of the incident electro- 
magnetic wave (the laser pulse) becomes u • = 270^ . The 
light is scattered by the electrons at the same frequency 
u)' r = lu[ (for a low intensity electromagnetic wave corre- 
sponding to the linear Thomson scattering and the undu- 
lator regime) or at harmonics u>' r — nw- (for high intensity 
corresponding to nonlinear Thomson scattering and the 
wiggler regime). When observed in the laboratory frame, 
this reflected wave is Doppler shifted and its frequency 
becomes oj r = 47 u>i (linear Thomson backscattering) or 
its harmonics (nonlinear Thomson backscattering) . This 
scheme offers the possibility to produce x rays in the keV 
range even with modest energy electrons (tens of McV) 
or gamma rays if 100 MeV range electrons are used. 

In the following sections the two schemes are discussed. 
For both, electron orbits, radiation properties, numerical, 
and experimental results are presented. 

A. Nonlinear Thomson scattering 

In this section, the production of femtosecond x-ray 
radiation from electrons initially at rest is discussed. In 
this scheme, Thomson scattering occurs in a strongly 
nonlinear regime and electrons are directly accelerated 
and wiggled within an intense laser field. This scheme 
offers the possibility to produce a large photon flux be- 
cause the number of electrons participating in emission 
can be much larger than the number of electrons trapped 
and accelerated in a laser-plasma accelerator. It is of 
the order of the number of electrons in the focal vol- 
ume. However, producing x-ray radiation requires a laser 
strength parameter a Q much greater than unity. Indeed, 
in that regime, the term — ev x B of the Lorentz force, 
negligible at low intensity, becomes comparable to the 
— eE component and the motion is a nonlinear function 
of the driving field in addition to becoming relativistic. 
The electron motion is no longer harmonic and the radi- 
ation emitted consists of high-order harmonics forming a 
broadband spectrum that extends up to the x-ray range. 
The spectrum shifts to higher energies as ag increases. 
This radiative mechanism is commonly called nonlinear 
or relativistic Thomson scattering. The principle of this 
source is displayed in Fig. 28. It simply consists of fo- 
cusing an intense laser on a target (here an underdense 
plasma is considered). 

In the following sections, the properties of nonlinear 
Thomson scattering are reviewed. The electron orbit is 



5 For an overview of the subject see Bardsley et al. (1989); 
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FIG. 28 Schematic of the nonlinear Thomson scattering 
source. A free electron is submitted to a relativistic laser 
pulse (ao S> 1). It is simultaneously accelerated and wiggled, 
and emits XUV radiation due to this motion. This illustrates 
the motion of an electron in the case of a linearly polarized 
laser pulse. 

derived in the case of a circularly polarized laser pulse, 
because it leads to a situation where the concept of wig- 
gling, developed in Sec. II, can be applied. For a linearly 
polarized laser pulse, analytical results can be easily ob- 
tained for the electron orbit, but not for the radiation 
features. In that latter case, the electron movement does 
not correspond to a transverse wiggling around a rel- 
ativistic drift motion along the propagation axis. The 
electron performs violent accelerations from 7 = 1 to 
7 = 1 + a.g/2 and decelerations to 7 = 1 for each half pe- 
riod of the movement. The emitted radiation comes from 
the longitudinal acceleration, and not from the transverse 
one anymore. Hence, in the Sees. VI. A. 1 and VI. A. 2 an- 
alytical expressions are derived for the circular polariza- 
tion but the linear polarization will be discussed in the 
numerical Sec. VI. A. 3. After a brief summary of the ex- 
perimental results obtained, we finally conclude with the 
short-term developments foreseen. 

1. Electron orbit in an intense laser pulse 

We consider a test electron initially at rest (7* = 1) 
submitted to an intense laser pulse modeled by a circu- 
larly polarized plane electromagnetic wave propagating 



Brown and Kibble (1964); Castillo-Herrera and Johnston (1993); 
Esarey et al. (1993a,c); Gibbon (1997); He et al. (2002); Lau 
et al. (2003); Salamin and Faisal (1996); Sarachik and Schappcrt 
(1970); Shen et al. (1997); Tian et al. (2006); Ueshima et al. 
(1999); and Vachaspati (1962). 
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along the e z axis, with a wave vector fej = 2it/\l e z and 
a frequency uii = 2ttc/\l, and with a normalized vector 
potential given by 



—= cos(u>it — kiz)e x H — p sm(u>it — kiz)e y 
\/2 V2 



The equation of motion of the test electron is 
dp 



dt 



= -e(E + vx B) 



(81) 



(82) 



All quantities with a hat are normalized by the following 
choice of units: m = c= e = a>j = l. The Hamiltonian 
describing the test electron dynamics in the electromag- 
netic wave is written as 



l+p : 



l + {P + a) 2 . (83) 



This system of an electron in a electromagnetic plane 
wave is integrable and action-angle variables can be found 
to solve the electron motion. Here, for simplicity, we 
analyze the symmetries of the system to directly yield the 
conserved quantities according to Noether theorem. The 

Hamiltonian H depends on the canonical momentum P 
and on the potential vector a(tp) through the variable tp = 
t—Z. Hence % is independent of x and y, which implies 
that the transverse canonical momentum is a constant of 
motion. The first constant of motion reads 



P 



p± — a = 0. 



(84) 



In addition, H, depends on t and z only through tp = 
t — z. Thus dH/dt = —dH/dz, which leads to the second 
constant of motion C, 



j-p z = C. 



(85) 



For an electron initially at rest, C = 1, and the test 
electron orbit is obtained by integrating the constants of 
motion, 



x(tp) 
2/0) 
z(tp) 
7 



«o 



sm(tp), 



-;-/~cos(<p), 



(86) 
(87) 



= 1 



(89) 



The electron motion in the laser field consists of a drift 
in the longitudinal direction combined with a transverse 
oscillation. The electron orbit is mainly longitudinal and 
relativistic for ao > 1 while it is essentially transverse and 
nonrelativistic for ao < 1. For a < 1, the test electron 
motion is a simple nonrelativistic dipole motion which is 
not relevant for the production of x rays since radiation 



is emitted at the same wavelength as the incident field. 
For ao > 1, the longitudinal motion becomes relativistic 
and nonlinear effects and Doppler shift occur. In the 
following, the interesting nonlinear case ao > 1, which 
results in the generation of smaller wavelength radiation 
than the incident laser pulse, is studied. 
The period of the motion is given by 



\ — \ 



(90) 



The electron orbit is similar to the standard orbits dis- 
cussed in Sec. II since the normalized energy 7 is con- 
stant. However, the trajectory takes place in three di- 
mensions and is helical. The qualitative regime and the 
angular distribution can be derived from the parameters 
Kx = li>x and Ky = ^yipY which are given by 



K 



x 



(1 + og/4) ~ a /V2, (91) 



where ipx and tpy are the maximum deflection angles of 
the orbits in the e x and e y directions. 



2. Radiation properties 

For ao > 1, the relativistic electron motion described 
above leads to the emission of nonlinear Thomson scat- 
tering radiation whose features can be calculated within 
the frame of the general formalism of the radiation from 
a moving charge. As seen in Sec. II, the spectrum of the 
produced radiation depends on the parameter K. Ac- 
cording to the above expressions of Kx and Ky , the mo- 
tion is driven in the wiggler regime since the nonlinear 
case ao > 1 is considered. The spectrum is broadband 
and extends up to the critical energy 



3 3 c 

huj c = x7 ft-, 
2 p 



(92) 



where p = {\ L /2w){a /V2 + V2al/16) ~ (A l /2tt) x 
\/2ao/16 is the instantaneous radius of curvature of the 
electron orbit obtained from the trajectory. Basically, 
the critical energy grows as a^ for a > 1. In practical 
units, it reads for ao 3> 1 



fiw c [eV] = 0.3 



\ L |)im] 



(93) 



Being an emission from a relativistically moving elec- 
tron, nonlinear Thomson scattering radiation is emitted 
in the direction of the electron velocity and has the same 
symmetry as the electron orbit. For circular polariza- 
tion, the radiation is emitted symmetrically around the 
laser propagation axis at an angle 9 ~ 2v / 2/ao, within a 
typical angular width A9 = I/7 ~ 4/a§. 

The number of photons can be estimated by integrat- 
ing, over one oscillation period, the expression of the ra- 
diated power P(t) = {e 2 /6Tre c)j 2 [(dj)/dt) 2 - (dj/dt) 2 } 
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FIG. 29 Top: The electron orbit in an intense linearly po- 
larized laser pulse with ao = 5 (red solid line) and ao = 10 
(green dashed line) . Bottom: The 7 factor as a function of the 
longitudinal position for each case. The electron is initially 
at rest and is submitted only to the linearly polarized laser 
pulse whose temporal profile is Gaussian with 20 fs duration 
(FWHM). 




FIG. 30 Spatial distribution of the nonlinear Thomson scat- 
tering, corresponding to trajectories of Fig. 29. The radiated 
energy per unit surface in the plane situated at 1 m from 
the source and perpendicular to the propagation axis is rep- 
resented. On the top is the case ao = 5 and on the bottom 
the case ao = 10. 



in which the trajectory calculated above is inserted. In 
practical units, the number of photons emitted at the 
mean energy (fku) = 0.3Hui c per one electron undergoing 
one oscillation period is, for ao 3> 1, 

N y = 4.68 x KT 2 a . (94) 

From the above expressions an estimation of the ra- 
diation properties can be obtained for a typical parame- 
ter regime. The laser pulse is focused along a Rayleigh 
length z r = 500 /xm with a — 10, and has a dura- 
tion of 20 fs FWHM and a wavelength X L = 800 nm. 
The period of the electron motion is A u ~ 20 /jm and 
the electron executes about ~ 8 oscillations within the 
laser pulse (FWHM) . The critical energy of the radiation 
is fujj c ~ 400 eV. The emission is emitted at an angle 
~ 280 mrad with respect to the laser propagation axis. 
For a plasma density of 1 x 10 18 cm~ 3 , and if we con- 
sider that the number of electrons participating to the 
emission is the number of electrons in the focal volume 
~ 10 x 10 x 1000 ^m 3 (which is equal to 10 11 electrons), 
the number of x-ray photons per shot is ~ 4 x 10 11 . 



3. Numerical results 

The features of the radiation have been described so 
far for a free electron submitted to a circularly polarized 
plane wave for which the concept of wiggling applies and 
analytical solutions are straightforward. In this section, 
the case of a linearly polarized laser field is studied us- 
ing numerical simulations. Electron orbits and radiation 
features are presented. 

We consider a laser pulse with a Gaussian temporal 



profile, whose normalized vector potential reads 

a = ao exp[— (2 ln2)(z — ct) 2 j (c 2 t 2 )] cos(u;.;i — kiz)e x . 

(95) 

The laser pulse duration (FWHM) is r = 20 fs and the 
polarization is linear along the e x direction. The test 
electron is initially at rest. Figure 29 displays the elec- 
tron orbits for ao = 5 and 10 and the evolution of the 
relativistic factor of the electron 7 along its motion. The 
trajectory consists in successive straight lines with vio- 
lent longitudinal acceleration or deceleration, the velocity 
vanishing at each corner. The electron acquires energy in 
the tens of MeV range while oscillating in the laser pulse, 
and is again at rest once the laser has passed. There is 
no net energy gain. 

The radiation has been calculated using the general 
formula (1) in which the orbits calculated above are in- 
serted. In Fig. 30, the spatial distribution of the radiated 
energy is represented. Because the polarization is linear, 
the radiation consists in two lobes corresponding to the 
two directions pointed by the straight lines of the orbit. 
These two lobes form, with respect to the laser propa- 
gation axis, an angle Ox ~ 450 mrad for ao = 5 and 
9x ~ 230 mrad for ao = 10. As expected, the radiation 
is more collimated for stronger ao- Figure 31 displays the 
spectra of the radiation integrated over the spatial dis- 
tribution (the photon number per 0.1% bandwidth per 
electron and per shot). An estimate of the total photon 
number can be obtained assuming that all electrons in 
the focal volume participate to the emission. Assuming 
a focal volume of dimension 10 x 10 x 1000 /an 3 and a den- 
sity of 10 18 cm -3 gives 10 11 electrons which radiate ~ 10 7 
photon/0. 1%BW at 50 eV for a = 5 and - 1.5 x 10 7 
photon/0. 1%BW at 400 eV for a Q = 10. Comparing 
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FIG. 31 Nonlinear Thomson scattering spectrum, integrated 
over angles and corresponding to the trajectories of Fig. 29. 
The number of photons per 0.1% bandwidth for one electron is 
displayed for the two cases: ao = 5 in red solid line and ao = 
10 in green dashed line. The spectrum near the maximum 
angle of emission is represented in the inset for the case ao = 
10. 



these numerical results with the analytical expressions 
obtained for the circular polarization, it appears that the 
spectral properties depend on the observation angle (see 
the inset of Fig. 31). The shape of the spectrum at the 
maximum angle of emission is not synchrotronlike, be- 
cause it does not correspond to a wiggling motion. The 
radiation for linear polarization is due to the longitudinal 
acceleration from 7 = lto7 = l + aQ/2 and deceleration 
to 7 = 1 occurring for each straight line section of the 
motion. The time profile of the radiated field in a single 
straight line section consists of a double peak structure: 
the first corresponds to the moment when the acceler- 
ation is maximal, and the second to the moment when 
the deceleration is maximal. This complex behavior is 
at the origin of the spectral modulation observed at the 
maximum angle of emission in the inset of Fig. 31 (Lee 
et al, 2003). 

It is important to note that the most simple case of 
nonlinear Thomson scattering has been considered. For 
a more accurate description of the mechanism, several 
effects must be taken into account. In the tightly fo- 
cused case, the exact laser field can strongly modify the 
electron orbit. In particular, the electron is expelled 
from the high-intensity region. This results in a radial 
drift and limits the number of electron oscillations in the 
high-intensity regions. Therefore, the radiated energy is 
smaller than in the ideal case and the divergence of the 
x-ray beam is higher and broadened. The presence of 
the ions from the plasma, or more generally, the collec- 
tive fields of the plasma, can as well have a significant 
impact on the electron orbit. They can reduce the longi- 
tudinal drift or accelerate the electrons. In the first case, 
the longitudinal drift being reduced, the divergence is in- 
creased. Oppositely, if the electron gains energy from the 
plasma, the radiation becomes more collimated. Finally, 
the laser pulse propagates in the plasma with a group 




FIG. 32 Experimental setup for the nonlinear Thomson scat- 
tering experiment. Radiation is produced from the interaction 
of a laser pulse with ao ~ 5 and a helium gas jet. Radiation is 
collected using a grazing incidence spherical mirror and spec- 
trally analyzed using a transmission grating. Filters are used 
to block the infrared light and to select a specific spectral 
range. To scan the angular distribution of the radiation, the 
whole system is mounted on a rotation stage centered on the 
gas jet. 

velocity v g smaller than c, which can strongly modify 
the electron orbit if the electron longitudinal velocity v z 
becomes comparable or greater than v g . 

4. Experimental results 

While nonlinear Thomson scattering was anticipated 
in the 1960s (Brown and Kibble, 1964; Sarachik and 
Schappert, 1970; Vachaspati, 1962), the first experimen- 
tal demonstration of nonlinear Thomson scattering was 
performed in 1998 only (Chen et al., 1998). Using a 4 
TW laser pulse with a duration of 400 fs, Al = 1 /im, 
and ao ~ 2, focused into a helium gas jet with a density of 
a few 10 19 cm -3 , nonlinear Thomson scattering radiation 
was observed for the first time. The measured radiation 
followed the specific features of nonlinear Thomson scat- 
tering and other radiative processes were ruled out. In 
particular, the spatial distributions measured were con- 
sistent with the nonlinear Thomson scattering properties 
obtained from the picture of the single free electron sub- 
mitted to a plane wave. In this experiment, radiation 
was detected at wavelengths up to the third harmonic of 
the laser (few eV range). The efficiency was estimated 
to be of the order of a few 10~ 4 photons per electron 
per pulse. In latter experiments (Bancrjcc et al, 2002, 
2003), the nonlinear Thomson scattering radiation was 
observed up to the 30th harmonic with a spatial distri- 
bution collimated in the direction of the laser propaga- 
tion axis (and maximum on axis, in contradiction with 
the single free-electron model). In 2003, thanks to the 
advent of high-intensity lasers, nonlinear Thomson scat- 
tering was demonstrated in the XUV range (Ta Phuoc 
et al., 2005b, 2003a,b). This experiment was based on 
the interaction of a 50 TW laser (with ao ~ 5, 30 fs, 800 
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the spectral range expected for nonlinear Thomson scat- 
tering. A total of 10 10 photons per shot were measured 
(integrated over the spectrum and the spatial distribu- 
tion). 
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FIG. 33 On the left is represented the angular dependence of 
the radiated energy recorded experimentally by rotating the 
whole detection system. The spectrum of the radiation on 
axis is given in logarithmic scales on the right of the figure 
and the inset shows the total number of photons integrated 
over several spectral bands determined by the choice of filters 
(the horizontal extension of the lines indicates the spectral 
band). 



nm) and a helium gas (density in the range 10 18 — 10 19 
cm -3 ). The experimental setup is presented in Fig. 32. 
The laser was focused within a 6 /jm focal spot (FWHM) 
onto the front edge of a supersonic helium gas jet (3 mm 
in diameter). Radiation was collected using a grazing 
incidence spherical mirror and spectrally analyzed using 
a transmission grating. The whole system was mounted 
on a rotation stage centered on the gas jet in order to 
measure the spatial distribution of the radiation. Non- 
linear Thomson scattering radiation was observed in the 
few tens of eV range up to 1 keV using a set of filters (Al, 
Zr, Ni, and Be). As at previous experiments, radiative 
mechanisms others than nonlinear Thomson scattering 
were ruled out by the variations of the XUV signal as 
a function of experimental parameters. The signal was 
found to increase linearly with the plasma density, the ra- 
diation was anisotropic and the energy range fitted what 
is expected for nonlinear Thomson scattering. 

However, the results differ from the expected emission 
with the model of free electrons in a plane electromag- 
netic wave. The radiation spatial distribution, presented 
in Fig. 33 (left), was collimatcd within a large angle 
(within a cone of half angle ~ 30 deg) and was maxi- 
mum on axis instead of consisting in two lobes centered 
at +23 deg and —23 deg. This result can be explained 
considering the multitude of electron orbits, the focusing 
of the laser and the effect of plasma fields. This was sup- 
ported by the fact that electrons at relativistic energies 
were observed during the experiment. In addition, laser 
filamentation was observed (Faure et al., 2002; Thomas 
et al., 2009, 2007) indicating a complex propagation of 
the laser. 

The spectrum, measured on axis in a limited spectral 
range (30 - 120 eV), is presented in Fig. 33 (right). The 
measurement was performed using a transmission grat- 
ing. In the inset is presented the spectrum measured with 
filters (Al, Zr, Ni, and Be) over a larger bandwidth and 
the number of photons. The radiation detected was in 



5. Perspectives 



Nonlinear Thomson scattering from the laser-plasma 
interaction remains a complex research topic and only 
a few experiments have been performed. Important 
numerical and experimental work has to be done to un- 
derstand the mechanism. This research will be motivated 
by the fact that this source of radiation can become 
a powerful x-ray source. Indeed, the development of 
petawatt-class lasers will open novel perspectives for 
nonlinear Thomson scattering. Assuming that do = 20 
can be reached, the produced radiation is expected to 
have the following features. The source is extended from 
the XUV to the x-ray part of the spectrum with photon 
energies up to ~ 10 keV. In addition, the number of 
photons will be significantly increased: ~ 1 photon per 
electron, per period, and per shot, and the radiation 
will be collimated within a half angle of ~ 150 mrad. 
Such a source has the main advantage to be very simple 
to realize. Further developments would rely on the 
control of the preacceleration of the electron inside the 
plasma in order reduce the divergence of the x-ray beam, 
enhancing the brightness. Flattop or annular transverse 
laser profile should be also interesting to maintain 
electrons in the high-intensity regions and increase 
the number of x-ray photons. Finally, several groups 
reported on the possibility to produce attosecond pulses 
via nonlinear Thomson scattering (Lan et al., 2006; Lee 
et al., 2003). In fact, the nonlinear Thomson scattering 
from a single electron is naturally a train of attosecond 
pulses (Lee et al, 2003). However, when summing the 
radiation from all electrons, the property disappears 
unless electrons radiate coherently (instead of incoherent 
radiation, as presented throughout this section). If this 
latter condition is satisfied, it becomes possible to gener- 
ate single intense attosecond pulse (Lan et al., 2006). To 
fulfill the coherence condition, the use of ultrathin solid 
targets has been suggested (Lee et al., 2005). A scheme 
named "lasetron" for generating coherent zeptosecond 
pulses has also been proposed by Kaplan and Shkolnikov 
(2002), and is based on the wiggling of electrons in a 
subwavelength-size solid particle or thin wire submitted 
to two counterpropagating circularly polarized petawatt 
laser pulses. 
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FIG. 34 Schematic of the ail-optically driven Thomson 
backscattering source. A laser-plasma accelerated electron 
beam is injected in a counterpropagating laser pulse. Elec- 
trons oscillate transversally and emit x rays or 7 rays due to 
this motion. 



B. Thomson backscattering 

X-ray or 7-ray radiation can be produced by scattering 
an electromagnetic wave off a counterpropagating rela- 
tivistic electron bunch. This scheme was first proposed in 
1963 (Arutyunian and Tumanian, 1963; Milburn, 1963) 
and reinvestigated in the 1990s thanks to the advance 
in laser technology and high-brightness electron accelera- 
tors, suggesting the possibility to realize a source of high- 
brightness pulsed x rays or 7 rays (Esarey et a/., 1996a, 
1993a,b; Gao, 2004; Hartemann, 2002; Hartcmann et ai, 
2001; Hartemann and Kerman, 1996; Hartcmann et ai, 
1996; Kim et al, 1994; Krafft et a/., 2005; Li et ai, 2002; 
Ride et ai, 1995; Sprangle et ai, 1992; Yang et at, 1999). 
More recently, when the research on electron acceler- 
ation from laser-plasma interaction became mature, it 
was proposed to use these laser-plasma accelerated elec- 
tron bunches to develop an all-optically driven scheme 
(Catravas et al, 2001; Hafz et ai, 2003; Hartemann et ai, 
2007; Lccmans et ai, 2005; Tomassini et ai, 2005). The 
principle of this scheme is represented in Fig. 34. Two 
laser pulses are required: the first drives the plasma ac- 
celerator, as discussed in Sec. Ill, and the second scat- 
ters off the accelerated electrons. In order to derive the 
typical properties and analytical features of the Thomson 
backscattered radiation, the electron dynamics in a coun- 
terpropagating laser pulse in the simplest model of plane 
electromagnetic wave is presented. Then, numerical and 
experimental results are discussed. Foreseen perspectives 
of this method are described in the conclusion. 



1. Electron orbit in a counterpropagating laser pulse 

We consider a laser-plasma accelerator producing a 
monoenergetic electron bunch of energy £ — Citric 2 with 
a velocity directed along e z (toward the positive val- 
ues). The counterpropagating laser pulse is modeled by 
a linearly polarized plane electromagnetic wave of wave 
vector ki = —2-k/Xl e z and frequency u>i = 2ttc/\l, 



and with a normalized vector potential given by a — 
aocos(u!it + kiz)e x . The equation of motion of the test 
electron and the Hamiltonian T~L describing the test elec- 
tron dynamics in the electromagnetic wave are given by 
Eqs. (82) and (83). We use the same normalized units, 
with the choice m = c= e = u>i = l. As in Sec. VI.A.l, 
the system is integrable and action-angle variables can be 
found to solve the electron motion. Here we analyze the 
symmetries of the system to directly yield the conserved 
quantities according to Noether theorem. The Hamilto- 
nian % depends on the canonical momentum P and on 
the potential vector a(tp) through the variable ip = t + z. 
Hence T-L is independent of x and y, which implies that 
the transverse canonical momentum is a constant of mo- 
tion. In our problem, it is assumed that the electron has 
a velocity directed along e z exclusively before entering 
the laser pulse. Therefore the first constant of motion 
reads 



P± = p± - a = 0. 



(96) 



In addition, % depends on i and z, only through (p = 
t + z. Thus dW/dt = dH/dz, which leads to the second 
constant of motion C, 

7 + Pz = C. (97) 

The constant C is obtained using the initial conditions 
7 = 7j, p z = \J ^il — 1 . The integration of the constants 
of motion gives the electron orbit: 



Hv) = ^sin(v3), 




l + al/2 



liv) = 77 + 



2C 2 j 

C 1 + <2g COS 2 (<p) 



8C- 



2C 



(98) 
(99) 

sin(2(^), (100) 
(101) 



where C = 7, + VtF^I = 2 7t - V(27i) + °(l/ 7l 2 )• 

The trajectory obtained cannot exactly take the form 
x(z) — K / ( r yk u ) sm(k u z) corresponding to the standard 
sinusoidal trajectory studied in Sec. II. Here 7 is not 
constant along the orbit and Eq. (97) implies that 7 and 
p z vary oppositely. The oscillation motion z takes, how- 
ever, a form similar to the sinusoidal trajectory, leading 
to the same characteristic figure-eight motion in the aver- 
age rest frame of the electron. The description of Sec. II 
can be recovered with the approximation ip ~ 2 z ~ 2 t. 
The counterpropagating laser pulse can be seen as an 
undulator with a spatial period given by 



A„ - A L /2, 
and a strength parameter K given by 



K = a = 0.855 d I[10 18 W/crn ]Af>m]- 



(102) 



(103) 
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2. Radiation properties 

The qualitative features of the radiation remain the 
same as in Sec. II and can be parametrized by K, A u , 
and 7 ~ ji (for ciq <C 7,). The spectrum of the emitted 
radiation depends on the amplitude of the parameter K. 
For K <C 1, the electromagnetic wave acts as an undula- 
tor. In the average rest frame of the electron, Thomson 
scattering occurs in the linear regime. In the laboratory 
frame, the radiation is emitted at the Doppler shifted fun- 
damental frequency corresponding to a photon energy hio 
in the forward direction (9 = 0). As K — > 1, harmon- 
ics of the fundamental start to appear in the spectrum. 
For if > 1, the electromagnetic wave acts as a wig- 
gler, and the spectrum contains many harmonics closely 
spaced and extends up to a critical energy fuj c (a nonlin- 
ear Thomson scattering occurs in the average rest frame 
of the electron). These energies are given by 

hw = {A 1 2 hc/\ L )/{1 + K 2 /2) for K < 1, 
huj c = 3Kj 2 hc/X L for K > 1. (104) 

This gives in practical units 

huj[eV] = 4.96 7 2 /A L [Mm] for K < 1, 



8xl0- 4 - 
4xl0- 4 - 




fiw c [eV] = 3.18 7 2 y7[10 18 W/crn] forif>l. (105) 

The radiation is collimated within a cone of typical 
opening angle 8 r = 1/7 in the undulator case. For a wig- 
gler, the radiation is collimated within a typical opening 
angle K/j in the electron motion plane (e x ,e z ) and I/7 
in the orthogonal plane (e y , e z ). 

As described in Sec. II, the duration of the emitted 
radiation is approximately equal to the duration of the 
electron bunch. The duration of the electron bunch being 
femtoseconds at the position of interaction (Lundh et a/., 
2011), the duration of the x-ray pulse is femtoseconds as 
well. 

The number of emitted photons follows the expressions 
of Sec. II. In practical units, the number of photons 
emitted per period and per electron (at the mean photon 
energy E = 0.3fiw c for the wiggler limit) is given by 



1.53 x 10~ 2 K 2 forK<l. 



AL = 3.31 x IQ- 2 K for K > 1. 



(106) 



From these expressions the features of the emitted ra- 
diation can be estimated. For a 100 MeV electron bunch 
(7 ~ 200) containing a charge of 100 pC and a counter- 
propagating laser pulse with ao = 0.5 and A^ = 800 nm, 
the electron oscillation period is \ u = 400 nm and the 
emitted radiation has an energy hui ~ 200 keV. The ra- 
diation is collimated within a cone of typical solid angle 
of 5 mrad x 5 mrad and 2 x 10 6 photons are emitted per 
oscillation period. Considering a laser pulse containing 
around ten optical cycles, the number of emitted photons 
is - 10 7 . 
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FIG. 35 Electron trajectories. The orbit is plotted as the 
red solid line for the undulator case (ao = 0.2) and as the 
green dashed line for the wiggler case (ao = 2). Before the 
interaction, the electron has 7 = 200 and p = \/^f 2 — le z 
in both situations. The counterpropagating laser pulse has a 
Gaussian temporal profile with FWHM duration of 20 fs and 
a central wavelength of 800 nm. 



3. Numerical results 

The electron equation of motion can be numerically 
integrated. The case of an electron initially propagat- 
ing in the e z direction (toward the positive values) with 
7 = 200 is investigated. The counterpropagating laser 
pulse is modeled by a linearly polarized plane wave with a 
Gaussian temporal envelope with = 20 fs at full width 
at half maximum (FWHM), at the wavelength 0.8 /j,m. 
Two laser strength parameters are considered: ao = 0.2 
for the undulator case and ao = 2 for the wiggler case. 

Figure 35 shows the orbits of a test electron traveling 
in the counterpropagating laser pulse for each situation. 
The transverse motion consists of an oscillation at a pe- 
riod X u = 0.4 /nm and its amplitude follows the envelope 
of the laser pulse and increases as ao increases. The max- 
imum amplitudes are 6.4 x 10 -5 and 6.4 x 10 -4 /jm for 
a = 0.2 and 2, respectively. 

The spatial distributions of the radiation produced by 
the electron along these trajectories are presented in Fig. 
36 (left). For ao = 0.2 (top), the radiation has a di- 
vergence of typical opening angle 1/7 ~ 5 mrad. For 
ao = 2 (bottom) , the radiation has a divergence of typ- 
ical opening angle ao/7 ~ 10 mrad in the plane of the 
electron motion and 1/7 ~ 5 mrad in the orthogonal di- 
rection. Figure 36 (right) represents the spectra of the 
emitted radiation, integrated over three selected solid an- 
gles for the undulator case a = 0.2 (1.5 x 1.5, 4x4, and 
10 x 10 mrad 2 ), and over the overall angular distribution 
for the wiggler case ao = 2. For ao = 0.2, the spectrum 
is nearly monochromatic at the energy hw ~ 200 keV in 
the forward direction, with a width of the emission line 
given by A A = X/N (where N ~ 7 — 8 is the number 
of oscillation periods in the FWHM laser pulse duration 
of 20 fs), and it is broadened after integration over the 
emission angles due to the angular dependence of the ra- 
diated wavelength A ~ (A u /27 2 )(1 + K 2 /2 + 7 2 6> 2 ) [see 
Sec. II, Eq. (5)]. For a = 2, the spectrum becomes 
broadband with a critical energy of about 300 keV. An 
estimate of the total photon number can be obtained on 
the basis of the number of electrons that are currently ac- 



37 




ao = 0.2, U ;S 



9- 

OB " 

-So 
J? a. — 2x|0- 6 



Wakefield cavity 
Relativistic e- bunch 



200 400 600 



o e - 8><io- ! 

40 ib 2 o 

. Q 

H ft U 

20Jr" 4x|0- ! 



ao = 2 




6x (mrad) 



400 800 1200 
Energy (keV) 



FIG. 36 On the left are represented the angular distributions 
of the radiation (radiated energy per unit solid angle) and 
the corresponding spectra are shown on the right (in num- 
ber of photons per 0.1% bandwidth, per electron, and per 
shot). The undulator case, ao = 0.2, is on the top whereas 
the wiggler case, ao = 2, is on the bottom. In the first sit- 
uation, three different solid angles of integration have been 
considered to highlight the broadening of the spectrum. In 
the wiggler situation, the integration is performed over the 
total angular distribution. 



celerated in a laser-plasma accelerator. Assuming 1 x 10 8 
electrons gives ~ 5 x 10 4 photons within a spectral band- 
width of ~ 10-15 % in the undulator case (selection of 
the 1.5 x 1.5 mrad 2 solid angle). In the wiggler case, 
- 2 x 10 3 photons/0. 1%BW at 600 keV and a total pho- 
ton number in the range 10 7 ~ 8 are expected. 



4. Experimental results 

Many experiments of Thomson backscattering have 
been performed using electrons from conventional ac- 
celerators. The first experiment to actually use laser- 
backscattered photons as a beam in a physics measure- 
ment was conducted at SLAC in 1969 (Ballam et al, 
1969). Many other experiments followed this approach 
to produce x rays and 7 rays (Albert et al, 2010; Ander- 
son et al, 2004; Babzicn et al, 2006; Bula et al, 1996; 
Burke et al, 1997; Glotin et al, 1996; Kotaki et al, 2000; 
Leemans et al, 1997; Litvinenko et al, 1997; Pogorel- 
sky et al, 2000; Sakai et al, 2003; Sandorfi et al, 1983; 
Schoenlein et al, 1996; Ting et al, 1995; Yorozu et al, 
2003, 2001), and nowadays several 7-ray facilities are 
based on this method. 

A significant milestone has been achieved with the 
first demonstration of femtosecond x-ray generation by 
Thomson backscattering (using electrons from conven- 
tional accelerators) in 1996, by Schoenlein and coworkers 
(Leemans et al, 1997; Schoenlein et al, 1996). They re- 
ported on the production of femtosecond x-ray pulses by 
crossing, at a right angle, a 50 MeV electron beam with a 
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FIG. 37 Principle of the all-optical Thomson source based on 
the combination of a laser-plasma accelerator and a plasma 
mirror. 



100 fs duration, 0.8 /xm wavelength, terawatt laser. The 
laser pulse contained 60 mJ and was focused on a 90- 
/im-diameter electron bunch (1.3 nC in a 20 ps FWHM 
bunch). In this experiment, ~ 5 x 10 photons were pro- 
duced with a maximum energy of 30 keV. However, due 
to the geometry and the mismatch of the electron bunch 
and laser pulse durations, the number of electrons par- 
ticipating in the production of radiation was limited to 
the portion of the bunch selected by the laser pulse. The 
x-ray pulse duration was theoretically predicted to be 
~ 300 fs due to the extension of the overlapping region 
and the radiation was found to be collimated within ~ 35 
mrad. 

Whereas many experiments have been performed us- 
ing conventional accelerators, only two experiments pro- 
duced Thomson backscattering radiation using an all- 
optically driven scheme (Schwoerer et al, 2006; Ta Phuoc 
et al, 2012). The first, performed in 2006, was based on 
the use of two laser pulses in a counterpropagating geom- 
etry (Schwoerer et al, 2006). The first pulse (85 fs, 800 
nm, and ao = 3) was used to produce an electron bunch 
with a Maxwellian spectrum whose temperature equals 6 
MeV and with a total charge in the picocoulomb range. 
The second laser pulse (85 fs, 800 nm, and a = 0.8) was 
counterpropagating and focused on the electron bunch. 
In this experiment, about 3 x 10 4 x-ray photons in the 400 
eV - 2 keV range were observed within a 80 /xsr detection 
cone at 60 mrad with respect to the laser axis. However, 
this type of experiment remains challenging because of 
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FIG. 38 Spectrum of the Thomson backscattering radiation 
measured using a set of copper and aluminum filters. The in- 
set is the x-ray beam angular profile measured using a phos- 
phor screen. 



x-ray signal transmitted by niters of different thicknesses 
and materials, is presented in Fig. 38. With the electron 
spectrum broadband, the produced radiation spectrum 
was broadband as well. The spectrum extends, as ex- 
pected, up to a few hundreds of keV and the total num- 
ber of photons is about 10 8 . The x-ray spectrum was 
reproduced numerically considering the measured elec- 
tron spectrum and charge, and a 15 fs (FWHM) dura- 
tion and a = 1.2 backreflected laser pulse. Figure 38 
(inset) presents the x-ray beam profile. The radiation 
was found to be collimated within a cone of divergence 
18 mrad FWHM. Finally, the source size, measured using 
the knife-edge technique, was estimated to be less than 
3 /itm (FWHM). The peak brightness is estimated to be 
of the order of 10 21 photons/(s mm 2 mrad 2 0.1%BW) at 
100 keV, which is more than 10000 brighter than Thom- 
son 7-ray sources from conventional accelerators (Albert 
et al, 2010; Schoenlcin et al, 1996). This high bright- 
ness is obtained thanks to the micron source size and the 
femtosecond duration reached in this all-optical scheme. 



the difficulty to overlap in time and space the electron 
bunch and the counterpropagating laser pulse. 

Recently, a much more simple and efficient method, 
based on the use of one laser pulse only, has been demon- 
strated (Ta Phuoc et al, 2012). This method has, for the 
first time, allowed the generation of ail-optically driven 
high-energy femtosecond x-ray beams. The scheme re- 
lies on the marriage of a laser-plasma accelerator and 
a plasma mirror. The principle is illustrated in Fig. 
37. An intense femtosecond laser pulse, focused into a 
millimeter-scale gas jet, drives a wakefield cavity in which 
electrons are trapped and accelerated. Then, the laser 
strikes a foil placed at nearly normal incidence with re- 
spect to the laser and electron beam axis. Ionized by the 
rising edge of the laser pulse, the foil turns into a plasma 
mirror which efficiently reflects the laser pulse (the reflec- 
tivity is up to 70%). Naturally overlapped in time and 
space with the backreflected laser pulse, the relativistic 
electrons oscillate in the laser field and emit a bright fem- 
tosecond x-ray beam by Thomson backscattering. 

The experiment was performed at the Laboratoire 
d'Optique Appliquee. A laser pulse (1 J, 35 fs, 810 nm, 
and ao — 1.5) was focused onto the front edge of a super- 
sonic helium gas jet (2.1 mm plateau and 600 /im gradi- 
ents on both sides) to generate the electron bunch. The 
electrons were accelerated in the forced laser wakefield 
regime (Malka et al., 2002). The electron spectrum was 
broadband with a maximum energy of 100 MeV and a to- 
tal charge of about 120 pC. A foil (either glass 1 mm or 
CH 300 /im) was placed in the region of the exit gradient 
of the gas jet. At the interaction of the backreflected laser 
pulse with the electron bunch, a bright high-energy fem- 
tosecond x-ray beam is generated via Thomson backscat- 
tering. The x-ray spectrum, obtained by measuring the 



5. Perspectives 

Thomson backscattering is a remarkable source for the 
production of energetic radiation, since it allows one to 
convert electrons into x rays or 7 rays in an efficient 
manner. Among the foreseen development, a promising 
route would be to take advantage of the latest progress 
made in laser-plasma acceleration. Monoenergetic elec- 
tron bunches with tunable energy have recently been pro- 
duced with an optical injection scheme (Faure et al., 
2006). The all-optically driven Thomson backscatter- 
ing source, based on these electron bunches, would allow 
one to produce tunable and nearly monochromatic radi- 
ation. This would represent significant progress since to 
date none of the existing laser-based femtosecond x-ray 
sources can produce monochromatic x-ray beams. Thom- 
son backscattering offers as well the possibility to gener- 
ate 7-ray radiation. For example, considering an electron 
energy of 1 GeV and a counterpropagating laser pulse of 
wavelength = 800 nm, 7 rays with an energy as high 
as 24 MeV can be produced. In this photon energy range, 
100 keV and beyond, radiography applications, of inter- 
est for industrial or medical applications, will be possi- 
ble. Taking advantage of the very small source size, high- 
resolution images can be obtained. Such high-energy ra- 
diation sources would also have many applications in nu- 
clear physics. An additional foreseen source development 
relies on the production of a high repetition rate x-ray 
source, delivering bright photon beams in the tens of keV 
range. In that case, electron bunches in the few tens of 
MeV will be produced using compact 100 Hz-kHz laser 
systems. Generating such a compact femtosecond x rays 
source at high repetition rate, accessible to a university- 
scale laboratory, will certainly have unprecedented im- 
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pact for applications in time-resolved x-ray science, high- 
resolution radiography, and phase contrast imaging. 

Finally, a fully optical free-electron laser scheme based 
on Thomson backscattering has also been proposed 
(Petrillo et ai, 2008) and will be discussed in Sec. VII. 
This scheme remains a feasibility study since it requires 
an electron bunch quality which has not been achieved 
yet. However, such a source, if it becomes realizable, 
would revolutionize work in the community of x-ray ra- 
diation. 



VII. COHERENT RADIATION: TOWARD A COMPACT 
X-RAY FREE-ELECTRON LASER 

In this section, the topic of free-electron laser (FEL) 
is introduced, especially in the XUV to x-ray range, and 
we discuss the perspectives to get FEL light from laser- 
accelerated electron bunches. The free-electron laser con- 
cept was first proposed by Madey (1971) who, having 
experience with insertion devices (undulators and wig- 
glers) for light sources, realized that a laserlike amplifier 
could be constructed by combining a high-quality elec- 
tron beam with an undulator or wiggler and an input 
radiation beam. A few years later, the principle was 
demonstrated experimentally by Madey and co-workers 
(Deacon et al., 1977). 

The electron beam inside the undulator acts as an ac- 
tive medium and provides the amplification of the radia- 
tion beam. The undulator is the coupling component, 
which permits energy exchange between the radiation 
and electrons. In practice, injecting an electron beam 
in an undulator will naturally lead to incoherent syn- 
chrotron emission, as discussed in the previous sections. 
This incoherent emission provides the input radiation 
beam and can be amplified by the FEL process when 
certain conditions are fulfilled (Bonifacio et ai, 1994, 
1984; Kondratenko and Saldin, 1980; Murphy and Pel- 
legrini, 1985; Saldin et at, 1998). The mechanism can be 
explained as follows. A single electron experiences the 
radiation from other electrons in the bunch. This inter- 
action between the bunch and its own radiation can lead 
to a microbunching of the electron distribution within 
the bunch at the fundamental wavelength of the radiation 
and its harmonics. This microbunching in return leads to 
coherent emission from the overall electron bunch, whose 
power is orders of magnitude higher than the usual inco- 
herent synchrotron radiation. However, this free-electron 
laser mechanism requires restrictive conditions on the 
electron beam quality and an important number N of 
oscillations to take place. This situation, where the inco- 
herent synchrotron or spontaneous emission is amplified 
in a single pass, is commonly named the self-amplified 
spontaneous emission (SASE) configuration, and brings 
a fluctuating radiation output (shot to shot) with poor 
temporal coherence. In the opposite way, injecting in the 



undulator an input coherent radiation beam, i.e., an ex- 
ternal seeding beam, gives a stable radiation output with 
good temporal coherence. Nevertheless, the latter case 
requires a coherent source at a given radiation wavelength 
with a power higher than the synchrotron one (Lambert 
et ai, 2009), and such coherent sources are not yet avail- 
able in the x-ray domain. Hence, the SASE configuration 
seems to be the simplest route toward the production of 
hard x rays. 

Many projects with a compact free-electron laser in the 
x-ray range of the spectrum (X-FEL) and operating FEL 
facilities are based on the SASE approach. The first facil- 
ity which has provided SASE radiation in the soft x-ray 
domain is the Free-Electron Laser at Hamburg (FLASH) 
which is currently working at wavelengths down to A = 6 
nm (Tiedtke et al, 2009). The Linac Coherent Light 
Source (Arthur et ai, 2002; Emma et al, 2010; Mc- 
Neil, 2009) and the Spring-8 Angstrom Compact Free- 
Electron Laser (Ishikawa et ai, 2012; Pile, 2011) have 
demonstrated FEL saturation at angstrom wavelengths, 
while the European X-ray Free-Electron Laser Altarelli 
et al (2007) is under preparation. On the other hand, 
seeded free-electron lasers are considered to improve the 
beam quality [for example, seeded FLASH (sFLASH) and 
FERMI@Elettra (Bocchetta et al, 2007)], but they gen- 
erally aim to less energetic photons (soft x rays). 

For our purposes, the electron beam is generated us- 
ing a laser-plasma accelerator and the undulator can be a 
conventional meter-scale undulator (see Sec. V), an elec- 
tromagnetic wave undulator (see Sec. VI), or the plasma 
undulator of the blowout regime (see Sec. IV). This latter 
differs significantly from the others because the strength 
parameter of the undulator depends on the amplitude of 
the betatron oscillation and on the electron energy. The 
amplification process has been studied for electrons in a 
ion channel, and has been named ion channel laser (ICL) 
for this case. In addition, the spectral ranges of interest 
are from XUV to x rays. 

In the following, we first present a short description of 
the underlying physics of the FEL process and summa- 
rize the conditions needed to be applied on the electron 
bunch for setting the FEL amplification. The two differ- 
ent modes of operation are discussed: the amplification 
of a seeding radiation beam and the amplification of the 
shot noise, i.e., of the spontaneous synchrotron radiation 
(SASE regime). Finally, the required conditions for the 
FEL amplification are used to investigate its feasibility 
from laser-accelerated electron bunches for three specific 
cases of undulator: the conventional one, the electromag- 
netic wave one, and the ICL one. We focus especially on 
the conventional undulator, since it represents the most 
promising route for near-future developments toward a 
compact x-ray free-electron laser. 
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A. The FEL amplifier 

1. Principle of the free-electron laser process 

To understand the physical mechanism of FEL ampli- 
fication, it is instructive to begin by the interaction of 
a single electron with the radiation beam and the undu- 
lator (Brau, 1990a; Huang and Kim, 2007; Saldin et ai, 
2000). In the following, a conventional undulator is con- 
sidered, i.e., a periodic structure of magnets generating 
a periodic static magnetic field (as shown in Sec. V), in 
order to present the FEL amplification process and the 
required conditions on the electron beam parameters. 

The undulator forces the electron to follow a sinusoidal 
trajectory whose period corresponds to the undulator pe- 
riod X u and whose amplitude is equal to K\ u / (27r7). The 
term K — e\ u Bo/(2i:mc) is the undulator strength pa- 
rameter, Bq is the peak magnetic field generated by the 
undulator, and 7 = £/mc 2 is the relativistic factor of the 
electron with £ the electron energy. The electron prop- 
agates in the z direction and the magnetic field of the 
undulator is oriented in the y direction, such that the 
motion of the electron is confined in the Oxz plane. In 
terms of velocity, the motion of the electron inside the 
undulator is written as 
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cos(fc u z), 
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where k u — 2w/X u , j3 x — j3.e x , and j3 = v/c is the elec- 
tron velocity normalized to the speed of light c. We then 
consider the effect of the radiation beam. This latter 
is modeled by a linearly polarized plane electromagnetic 
wave whose electric field is E = Eosm(kz — uot)e x (one- 
dimensional analysis), where k — 2tt/X is the norm of the 
wave vector, uj = 2ttc/\ is the frequency, A is the radi- 
ation wavelength, and Eq is the electric field amplitude. 
Assuming that the radiation field is small compared to 
the undulator field {K 3> a® — eEo/mcuj), the expres- 
sion of the electron velocity (107) can be used in order 
to obtain the rate of electron energy change using the 
kinetic energy theorem, 
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where ip = kz — Lot + k u z is the ponderomotive phase. An 
energy exchange can occur between the electron and the 
radiation if the sine terms are nonzero after averaging 
over z. The phase ip can remain constant if the increase 
of the undulator phase k u z compensates the decrease of 
the radiation phase kz — tot, so that the first sine term is 
manifestly nonzero after averaging. This latter condition 
is written, for the phase ip = kz — ut + k u z defined with 



the averaged arrival time t = J dz/v Zl as 
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The above resonance condition states that when the elec- 
tron advances one undulator period, the radiation field 
has slipped by one radiation wavelength A. It can be put 
in the following form: 
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We emphasize here that this equation is nothing else but 
the fundamental wavelength of the electromagnetic field 
radiated by the electron moving in the undulator [see Eq. 
(5) in Sec. II]. This reinforces the physical intuition that 
energy exchange can occur between the electron and the 
radiation when the condition (110) is fulfilled. A more 
careful analysis of the sine terms of Eq. (108) [see Brau 
(1990b) and Wiedemann (2007b)] shows that energy ex- 
change can also occur if A n = (A u /2n7 2 )(l + K 2 /2), i.e., 
for harmonics of the electromagnetic field radiated by the 
electron. This exchange at harmonics of the fundamental 
wavelength is possible because of the z oscillatory motion 
of the electron in the undulator. In addition, the analysis 
shows that Bessel functions appear in the multiplicative 
factor, so that near the resonance at the fundamental 
wavelength Eq. (108) becomes 
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where A u = K[J (x) - Ji{x)]/y/2 with x = !f 2 /[4(l + 
K 2 /2)], Jo and J\ are Bessel functions (for K <C 1, A u ~ 
K/y/2 and for K > 1, A u ~ K/2), and the bar over %p 
is now implicitly assumed. In the following, we will not 
consider harmonics and concentrate on the resonance at 
the fundamental wavelength. 

The ponderomotive phase tp plays a major role in the 
FEL dynamics and it is relevant to locate the electron 
using the variable ip instead of z, and to use z as the new 
time variable. If £q is the energy at resonance satisfying 
(110) and A£ = £ — £ is the energy of the electron 
relative to the resonance energy, then for A£ <C £0 Eqs. 
(109) and (111) become 
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where j z = 1/ y 1 — [3 Z 2 . By differentiating Eq. 
the pendulum equation is obtained, 
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Hence, the electric field of the radiation Eg (all other pa- 
rameters being constant) can be interpreted as the am- 
plitude of the effective or ponderomotive potential of the 
pendulum. There are two classes of trajectory for the 
pendulum equation: trapped and untrapped. The curve 
in phase space (ip, £) at the limit of both regimes is the 
separatrix. The higher the electric field is, the higher the 
phase space volume of trapped conditions (volume inside 
the separatrix) and the induced acceleration or decelera- 
tion for the electron. 

This analysis of the interaction of a single electron with 
the radiation beam and the undulator can be used to de- 
scribe qualitatively the collective effect of FEL amplifi- 
cation occurring with a bunch of electrons. The initial 
condition is the electron distribution in the phase space 
(ip, £). As discussed in Sec. II. F, the electron bunch 
coming from a laser-plasma accelerator has a random 
distribution. For simplicity we consider the case of a 
uniform distribution in phase tp (such that there is no 
synchrotron radiation) and electrons at exact resonance 

£ = £ . 

Imposing that the electric field of the radiation Eq is 
constant (the so-called small signal gain regime), there 
are as many electrons which lose energy as electrons 
which gain energy. As a result, there is no energy gain 
for the radiation. In this regime, increasing the radia- 
tion energy requires £ > £q and decreasing the radiation 
energy requires £ < £q. 

However, the FEL amplification takes place in the 
high-gain regime and the process significantly differs. 
When the electrons of the initially uniform beam start 
to move inside the ponderomotive potential, the electron 
distribution is modulated with a period 2tt in ip, which 
corresponds at a fixed time t to a modulation with period 
A in z. This small modulation implies that the electron 
beam radiates partially coherently such that the electric 
field of the radiation Eq increases and the mean energy 
of the electron beam decreases [in Eq. (113), the en- 
ergy variation averaged over all electrons is nonzero be- 
cause (sin 7/;) 7^ 0]. As a consequence, on the one hand, 
the amplitude of the ponderomotive potential increases 
and, on the other hand, electrons lose more energy on 
average than predicted by the single electron model. Be- 
cause the ponderomotive amplitude increases, the volume 
of trapped conditions also increases and electrons which 
initially were on untrapped trajectories can be finally 
trapped in the ponderomotive potential well. Therefore, 
after a stage of modulation of the electron distribution, 
electrons tend to be trapped in the wells of the pon- 
deromotive potential. The amplification process can be 
summarized as the following: 

• The electromagnetic radiation induces a modula- 
tion of the electron beam. At the beginning, the 
modulation is sinusoidal and the regime is linear. 




Position inside the undulator z (arb. units) 



FIG. 39 Schematic of the FEL amplification process. A typ- 
ical profile for the electric field amplitude Eo along the un- 
dulator axis in linear scale (red solid line, left y axis) and 
logarithmic scale (green dashed line, right y axis). 

The higher the electric field amplitude Eq is, the 
faster the amplitude of the sinusoidal modulation 
grows. 

• The electric field amplitude growth is proportional 
to the beam modulation. Hence, it results in an 
exponential amplification, Eq oc exp(z/2Z g ), where 
l g is the gain length in power. 

• The modulation of the electron density becomes on 
the order of the mean electron density, i.e., becomes 
nonlinear. This happens when electrons become 
trapped in the wells of the ponderomotive poten- 
tial. The electron beam is microbunched at the ra- 
diation wavelength A and radiates coherently. The 
largest fraction of output radiation is generated 
during this last stage. This leads to an important 
loss of energy for the electron beam in such a way 
that it becomes off-resonance. This corresponds to 
the saturation: the radiation energy cannot be in- 
creased further in normal operation. 

• In addition, letting the process go on, electrons off- 
resonance fall in accelerating phases and extract 
energy from the electromagnetic radiation. When 
they recover sufficient energy to satisfy the reso- 
nance condition, they can amplify the radiation 
again. Hence, the radiation amplitude starts to os- 
cillate after the saturation point. 

This FEL amplification process is illustrated in Fig. 39, 
where the radiation amplitude is plotted as a function 
of z. This emphasizes the different stages of the process: 
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linear mode of amplification, nonlinear mode, saturation, 
and then oscillations. 

The efficiency of the amplification process, defined as 
the ratio of the output total radiation energy to the 
total energy of the electron bunch, is limited at satu- 
ration by the fact that electrons lose energy and be- 
come off-resonance. As presented next, in the ideal 
fD case the efficiency at saturation equals the Pierce 
parameter p which is small compared to unity in all 
cases. Recalling that the resonance condition reads 
A = (A„/27 2 )(f + K 2 /2), the loss of electron energy (de- 
creasing of 7) can be compensated by varying undulator 
parameters (for example, the undulator gap and so K) 1 
which allows the process to remain in resonance after the 
saturation point and to continue to radiate at the wave- 
length A in a fully coherent manner. This approach is 
called tapering and it permits the FEL efficiency 77 to be 
increased well above the untapered one r\ ~ p and even 
to approach efficiency of unity. 



2. Required conditions on the electron beam parameters 

In our simple illustration of the amplification process, 
any practical effects such as energy spread in the elec- 
tron beam have not been considered and tolerances of 
the FEL process to deviations from the ideal case have 
not been discussed. Here the conditions requested from 
the electron beam parameters in order to allow amplifi- 
cation of an electromagnetic field are summarized. The 
Pierce parameter p (also named efficiency parameter and 
FEL universal parameter) is defined as 



A u ^u 
4A„ 



2/3 



I 

2^ 



1/3 



A u \ u 

27TCT 



2/3 



(115) 



where A p = 2nc/ yjn e e 1 /meo is the plasma wavelength 
(without relativistic correction), / is the electron beam 
current, I a — income 3 / e ~ 17.0kA is the Alfven current, 
a is the rms radius of the transverse distribution of the 
electron beam, and A u — K[J (x) — Ji(x)]/y/2 is defined 
in Sec. VII. A. 1. In practical units, the Pierce parameter 
reads 
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a. Steady-state ID FEL theory 

In the steady-state ID FEL theory, the Pierce parame- 
ter provides the following conditions applied on the elec- 
tron beam and the output values of the main parameters: 



the rms relative energy spread of the electron beam 
has to be smaller than p, 
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• the relative detuning of the mean energy £ m of the 
electron beam from the resonance (£q is the en- 
ergy at resonance corresponding to a given radia- 
tion wavelength A) has to be smaller than p, 
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the power gain length l g [defined by P cx exp(z/l g ), 
where P is the radiation power] is given by 
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• the efficiency of energy conversion from the electron 
beam to the radiation 77 = 7 ra diation /{N e ^mc 2 ) is 
finally given at saturation by 



(121) 



The output radiation does not depend on the energy of 
the input radiation if the length of the undulator is cho- 
sen in such a way that saturation is attained at its exit. 
This is because, whatever the input radiation is, the pro- 
cess ends when electrons have lost sufficient energy to be 
off-resonance, and the amount of energy that they lose to 
become off-resonance is roughly given by p£. The total 
radiation energy /radiation scales as N^ 3 since the Pierce 
parameter p depends on the electron current as J 1 / 3 . 
The undulator length providing the saturation depends 
on the energy of the input radiation, and there is no sim- 
ple expression to estimate it. Nevertheless, its typical 
value is on the order of 20 gain lengths, Z sat ~ 20 l g . 
According to the above formulas, the Pierce parameter 
p has to be as high as possible. Indeed, this leads to 
less stringent conditions on the energy spread and detun- 
ing, smaller gain length, and so smaller undulator length 
and higher efficiency at saturation. The parameter p in- 
creases with the current or density as I 1 / 3 or n~J 3 , with 
the strength parameter as A^/ 3 and decreases with 7 as 
I/7 for constant undulator period A u . However, for con- 
stant radiation wavelength A, p increases with 7 as 7 1 / 3 
(the dependence on current or density is the same, while 
the dependence on K becomes more complex). 

Our next step will be to determine the limits of valid- 
ity of the steady-state ID FEL theory. Working inside 
the limits and the conditions of the ID theory constitutes 
the best configuration for a high FEL performance. For 
out-of-limit configurations, one can use the Ming Xie for- 
mulas (Xie, 2000) (given below) which coarsely account 



43 



for 3D effects, or simulations with codes such as Gen- 
esis (Rciche, 1999) for example, which give precise and 
quantitative answers. 

b. Three-dimensional effects 

Three-dimensional effects and the associated limits of 
the ID theory arise when considering the transverse di- 
mensions. First, do the electron beam and the radiation 
correctly overlap? Indeed, in three dimensions diffrac- 
tion occurs and radiation can escape the electron beam 
and become lost for the amplification process. To avoid 
this, the Rayleigh length of the radiation has to be higher 
than the gain length (to ensure gain guiding). Second, 
the electron beam having a finite emittance, it has an an- 
gular spread: electrons have slightly different directions, 



resulting in different values for 7 Z = 1/y 1 — j3 z and in 
an effective energy spread. Imposing that this effective 
energy spread is smaller than p provides a condition on 
the unnormalized rms transverse beam emittance e (see 
Sec. II. H). The above conditions or limits of the ID the- 
ory can be written as 
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The Ming Xie formulas allow these effects on the FEL 
process to be estimated, by giving the modified gain 
length and the modified efficiency at saturation (Xie, 
2000), 
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where A^, A e , and A 7 are parameters taking into account, 
respectively, diffraction, emittance, and energy spread ef- 
fects, and are given by 
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These formulas show how 3D effects degrade the FEL 
performance, increasing the gain length and decreasing 
the efficiency. However, they do not take into account 
time-dependent effects and space-charge effects. 

c. Time- dependent effects 



Time-dependent effects describe how the short dura- 
tion of electron bunches and variation of electron beam 
parameters inside the bunch influence the FEL amplifi- 
cation process. The previous formulas assume a steady- 
state electron beam, i.e., infinitely long with constant 
density, for which slippage between radiation and elec- 
trons can be neglected. Realistic bunches have finite du- 
ration and the electron density varies along the bunch. 
The cooperation length is defined as the distance slipped 
by the radiation with respect to electrons in one gain 
length, 
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Electrons interact with each other via the radiation 
field. However, this radiation field slowly slips on the 
bunch. The cooperation length indicates the size inside 
the bunch where interaction between electrons can occur 
during a gain length. As a direct consequence of this 
notion, one could consider that the bunch is made of in- 
dependent sub-bunches of length l c . 
The steady-state model is valid if the electron beam pa- 
rameters are approximately constant in a cooperation 
length inside the bunch, and if the electron bunch length 
is much longer that the cooperation length. This latter 
condition reads 
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where r is the duration of the electron bunch. 
For bunch durations smaller than the above limit, the 
FEL process is highly degraded because radiation es- 
capes the bunch in less than a gain length (edge effects) . 
This regime is called weak superradiance (Bonifacio and 
Casagrande, 1985; Bonifacio et al, 1990, 1988, 1989). Su- 
perradiance, as defined by Dicke (1954), is a spontaneous 
emission from a coherently prepared system, whose en- 
ergy scales as the square of the number of electrons. FEL 
superradiance is slightly different, because electrons be- 
gin to bunch when interacting with the spontaneous emis- 
sion and evolve toward a coherently prepared system. 
The radiation energy of the weak superradiance regime 
then scales as N 2 instead of the N^ 3 steady-state depen- 
dence, but is much lower than the steady-state radiation 
energy. In contrast, for long electron bunches, strong 
superradiance can occur if the weak superradiance radia- 
tion emitted by the trailing edge of the bunch propagates 
and slips through the overall bunch and is subsequently 
amplified. In that case, the peak power is much greater 
than the steady-state saturation value. 

The notion of cooperation length also indicates that 
the condition on the energy spread does not have to be 
applied on the overall electron bunch. It is sufficient that 
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each sub-bunch of length l c has an energy spread smaller 
than p. The energy spread calculated within a coopera- 
tion length is called the slice energy spread. The condi- 
tion (118) is relaxed and replaced by 

(f) r <». (133) 

\ C / slice 

Requiring transverse emittance can also be replaced with 
the same conditions applying slice transverse emittance, 
defined in the same way as slice energy spread. 

d. Space-charge effects 

Another effect which influences the FEL performance 
is space charge. When the FEL amplification process 
tends to microbunch the electron beam at the radiation 
wavelength, longitudinal space-charge forces, on the con- 
trary, tend to debunch the beam. The inverse plasma 
frequency w" 1 (time scale of plasma oscillations) can in- 
dicate over which distance debunching occurs, while the 
gain length characterizes the distance over which bunch- 
ing forces are significant. In the electron rest frame, the 
longitudinal distance is 7 times higher than in the labora- 
tory frame, so the density is divided by 7 and the plasma 
frequency by 7 1 / 2 . Hence plasma oscillations take place 
on a time scale 7 1 / 2 /w p in the rest frame, and because of 
time dilatation, they take place on a time scale 7 3//2 /cj p 
in the laboratory frame, corresponding to a covered dis- 
tance of 7 3 / 2 c/w p . Longitudinal space-charge forces are 
then negligible when 

7 3/2 c/w P > l g , (134) 

where ui p = \J n e e 2 /meo is the plasma frequency without 
any relativistic correction. 

Transverse space-charge forces can also influence the FEL 
performance. In the envelope equation for the elec- 
tron beam, the focusing field, emittance term, and the 
space-charge term appeared. To neglect transverse space- 
charge forces, the corresponding term in the envelope 
equation has to be small compared to the other terms, 
i.e., the focusing one and the emittance one (Humphries, 
1990). 

e. Undulator errors and wakefields 

The practical realization of a free-electron laser con- 
sists in many undulator sections and transport devices. 
The magnetic field of undulators has small errors and 
sections can be slightly misaligned, resulting in a degra- 
dation of the FEL performance. In addition, a high- 
current electron beam propagating between the two re- 
sistive walls of small-gap undulators induces a wakefield 
which can modify the beam properties along the propa- 
gation, altering the FEL process. See Huang and Kim 
(2007) for a description of these effects. 

/. Quantum FEL 



The FEL process described above is based on a clas- 
sical treatment, so it is relevant to ensure that quantum 
phenomena are unimportant by determining the limit 
for which the dynamics is effectively described by clas- 
sical mechanics. The FEL dynamics becomes quantum 
if the momentum recoils Ap = —Hk of an electron emit- 
ting a photon of wave vector k is greater than the FEL 
bandwidth p\p\ = pE/c, i.e, if the electron becomes off- 
resonance after emitting a photon. Hence, in a quantum 
FEL, each electron radiates only one photon. In all real- 
izable cases, hk <C pp and quantum effects are completely 
negligible so that the classical FEL theory correctly de- 
scribes the process. 

To conclude, we emphasize that ID formulas are use- 
ful for simple estimations, and Ming Xie formulas (Xie, 
2000) for corrections due to energy spread, diffraction, 
and emittance. For quantitative study of a given config- 
uration and for time-dependent effects and space-charge 
force effects, self-consistent, 3D, and time-dependent 
codes such as Genesis can be used (Reiche, 1999). 

3. Seeding or self-amplified spontaneous emission 
configurations 

In the above description of the FEL amplification pro- 
cess we did not discuss where the input radiation field 
came from. The FEL radiation can be started either by 
an initial seeding field or by a random modulation of the 
electron beam (SASE regime). 

The above description of the FEL amplification process 
directly applies to the seeding case. This case provides 
high-quality output radiation, but requires that a rela- 
tively intense coherent radiation field at the given wave- 
length A can already be generated, which is not the case 
for high-energy photons as x rays. In principle, high- 
order harmonics (HHG) generated in gas allow one to 
seed FEL up to the extreme ultraviolet part of the spec- 
trum at the present state of the art. A recent experi- 
ment has demonstrated the FEL amplification of a seed- 
ing beam at 160 nm by 3 orders of magnitude (Lambert 
et at, 2008). These schemes, and, in particular, HHG 
on solid targets (Dromey et at, 2006; Tarasevitch et al., 
2000; Teubner and Gibbon, 2009), could be extended to 
x-ray photons in the near future and could provide a seed 
for X-FEL. However, the seeding power has to be higher 
than the power of the incoherent synchrotron or sponta- 
neous emission (Lambert et al., 2009), and this latter is 
proportional to the desired photon energy. This means 
that it will be more and more difficult to have a suit- 
able seeding beam for smaller radiation wavelength and 
higher photon energy because the required power will be 
higher. 

The case of an initially randomly modulated electron 
beam (SASE configuration) is different and the theoret- 
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ical description is more complex. No radiation input en- 
ters with the electron beam inside the vmdulator. As dis- 
cussed in Sec. II. F, electrons inside bunches delivered by 
conventional or laser-plasma accelerators are randomly 
distributed in space. Whereas if the electron distribu- 
tion in space is uniform no radiation is emitted, for a 
random one, incoherent radiation is emitted (synchrotron 
radiation, whose power is proportional to the number of 
electrons) and this incoherent radiation provides a seed 
for the FEL amplification process. Moreover, the inco- 
herent radiation is emitted at the wavelength given by 
Eq. (110), and therefore the resonance condition is auto- 
matically satisfied. The main advantage of this method 
is straightforward: it allows one to extend FEL radia- 
tion to high-energy photons because no external coher- 
ent radiation is needed. A major experiment demon- 
strating self- amplified spontaneous emission gain, expo- 
nential growth, and saturation was achieved by Milton 
et al. (2001) and highlighted the possibility to develop 
an operational X-FEL. Nowadays, most of the X-FEL 
projects (A < 1 nm) are based on this approach. How- 
ever, because the input is a fluctuating quantity shot to 
shot, the quality of the output radiation is very low com- 
pared to the seeding case, in which the output inherits 
the quality of the external input. For the SASE case, 
the output power has very large fluctuations and very 
poor temporal coherence. The temporal profile and the 
spectrum are constituted of many spikes; the spikes of 
the temporal profile have typical durations of l c /c while 
the spikes of the spectrum have typical relative width 
of X u /lb, where 4 is the bunch length (the FEL process 
smooths the parameters of the electron beam and the 
radiation over the scale of the cooperation length) . 



B. Free-electron laser from a laser-plasma accelerator 

From the different required conditions on the electron 
beam parameters (see Sec. VII. A. 2), the feasibility of 
a FEL from a laser-plasma accelerator can be investi- 
gated for each specific undulator or wiggler type. The 
above results which are strictly valid for the conventional 
undulator give only good estimations for the other ones 
[for example, see Bacci et al. (2008) for the comparison 
between electromagnetic wave and conventional undula- 
tors] . 

1. With a conventional undulator 

Conventional undulators have typically centimeter pe- 
riod and strength parameter of unity. Hence, according 
to the resonance condition (110), using electron bunches 
at the 100 MeV level, the FEL radiation wavelength is 
in the UV to XUV range, whereas at the 1 GeV level it 
reaches the soft x-ray and hard x-ray ranges. In the fol- 



lowing, we consider two different scenarios, X-FEL 1 and 
2, which are speculative and can be seen as future per- 
spectives in the short term and long term, respectively, of 
laser-plasma accelerators. All input parameters of these 
scenarios and the corresponding output FEL parameters 
are given in Table I. Table I contains the results obtained 
by 3D time-dependent Genesis simulations in the SASE 
configuration. 

The choice of the electron beam parameters is moti- 
vated by extension of the 100 MeV range state-of-the- 
art laser-plasma accelerator performance to higher elec- 
tron energies. The electron beam peak current is main- 
tained constant when going to higher electron energies 
because the requirement of small energy spread imposes 
an optimized beam loading which flattens the accelerat- 
ing electric field along the bunch (Rechatin et al. , 2009a) . 
The optimized peak current scales as the normalized vec- 
tor potential ao of the laser pulse and is independent of 
the electron density of the plasma n e . Because the elec- 
tron accelerator efficiency scales as 1/ao, ao, and there- 
fore the electron beam peak current, should not be in- 
creased when designing a higher energy laser-plasma ac- 
celerator. The absolute energy spread and the normal- 
ized emittance can remain constant in the acceleration 
process if beam loading is optimized and if quasista- 
tionary forces are applied to the electron beam. There- 
fore, in principle, the low absolute energy spread (a few 
MeV) (Rechatin et al., 2009b) and normalized emittance 
eyv = 1 7r.mm.mrad (Brunctti et al, 2010) achievable for 
electrons at the 100 MeV range could also be obtained 
at 1 and 5 GeV. 

Note that the bunch rms transverse size a is exper- 
imentally chosen by the transport configuration. Of 
course, this choice determines the rms angular spread 
6 = e/a. In order to have the Pierce parameter as high 
as possible, a has to be small, but, on the other hand, 
emittance, diffraction, and space-charge effects increase 
while a decreases. In addition, practical considerations 
on the transport realization along the undulator put a 
limit on the smallest transverse size achievable. 

Since the topic of short wavelength free-electron laser 
using conventional undulators has been extensively stud- 
ied by means of conventional accelerators, we discuss here 
the differences which arise when using electron bunches 
produced by laser-plasma accelerators. The first differ- 
ence is straightforward: because of the short duration of 
the electron bunch, currents can be higher which result in 
higher p, smaller gain lengths, smaller undulator lengths, 
higher efficiency at saturation, less stringent conditions 
on the energy spread and detuning, but also in higher 
space-charge effects. 

Second, the short duration of the electron bunch can be- 
come comparable or smaller than the cooperation length 
divided by c, leading to important edge effects which de- 
grade the FEL performance. 
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TABLE I Parameters of different scenarios for realizing free-electron lasers from laser-driven plasma-based accelerators in the 
XUV and x-ray domains with the corresponding FEL quantities obtained from Genesis simulations in SASE configuration. For 
all scenarios, the undulator period X u is equal to 1 cm and the strength parameter K is equal to 1. 





X-FEL 1 


X-FEL 2 


Electron beam input parameters 


Electron energy 


1 GeV 


5 GeV 


Peak current 


10 kA 


10 kA 


Bunch rms duration 


4 fs 


4 fs 


Bunch rms transverse size 


30 /im 


30 /im 


rms energy spread 


0.2% 


0.04% 


Normalized rms transverse emittance 


1 7r.mm.mrad 


1 7r.mm.mrad 


FEL output parameters 


Radiation wavelength A 


1.96 run 


0.0783 nm 


Pierce parameter p 


0.00225 


0.000450 


Undulator length 


12.5 m 


100 m 


Peak power 


25 GW 


15 GW 


Pulse energy 


50 nJ 


25 ,uJ 


rms pulse duration 


2 fs 


2 fs 



Third, the electron bunch transport from the plasma 
source to the undulator is specific: even if the emittance 
is comparable to that of conventional accelerators, at the 
exit of the plasma, the transverse size is very small (a 
few microns) when the divergence is higher (a few mil- 
liradians). The transverse size grows fast: 10 cm after 
the plasma it reaches a few hundreds of microns, and 1 
m after it reaches a few millimeters. For such size, chro- 
matic aberrations of quadrupoles (which are proportional 
to the electron beam energy spread) are large, leading to 
an enormous growth of emittance, and because electrons 
cover different distances, the longitudinal bunch length 
also grows. It is therefore necessary to design an ultra- 
compact transport system which acts near the plasma 
source and does not let the beam transverse size grow. 
Space-charge forces, higher in the case of laser-plasma 
accelerators, also strongly affect the beam transport, in- 
creasing emittance, longitudinal bunch length, and total 
energy spread [they induce a linear energy chirp: the 
head of the bunch moves faster than the tail (Griincr 
et al, 2007, 2009)]. 

We comment now on the results obtained for the sce- 
narios considered in Table I. They are illustrated in 
Fig. 40, which shows the longitudinal profile of radiation 
power as a function of the position inside the undula- 
tor for each scenario. In both configurations, X-FEL 1 
and 2, the small values of the Pierce parameter impose 
stringent conditions on the electron beam energy spread. 
As can be seen in Table I, the energy spread is close to 
p in both cases, and using a higher electron beam en- 
ergy spread strongly reduces the FEL performance. In 
X-FEL 1, the output FEL radiation consists of a few 
SASE spikes, meaning that the bunch length is of the 
order of a few cooperation lengths and that edge effects 
do not degrade the FEL performance. Moreover, diffrac- 
tion, emittance, and longitudinal space-charge effects are 




Position relative to the electron beam center s (/im) 




Position relative to the electron beam center s (jim) 

FIG. 40 Longitudinal profile of radiation power as a function 
of the position inside the undulator obtained from Genesis 
simulations. The corresponding electron beam and undulator 
parameters are indicated in Table I. Top: X-FEL 1. Bottom: 
X-FEL 2. The position relative to the electron beam center 
is defined as s = ct e — ct, where t e and t are, respectively, the 
arrival time of the electron beam center and of the radiation, 
at a given position z in the undulator. 
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found to have negligible impact on FEL output, while 
the main degrading effect comes from the energy spread. 
Strong FEL radiation, with 25 GW peak power and 50 
/iJ pulse energy, is produced at a wavelength of 1.96 nm 
with an undulator length of 12.5 m. 
In X-FEL 2, the radiation wavelength reaches the 
angstrom range, a spectral region of great interest for 
diffraction experiments, more particularly on biological 
samples. In these conditions, the bunch length equals 
many cooperation lengths and the notion of slice energy 
spread is helpful: it is sufficient that only the slice en- 
ergy spread is below 0.04%, not the total energy spread. 
Here, because of the higher 7 and the smaller p, the re- 
quirements on the relative slice energy spread and on the 
undulator length are more stringent. The output FEL ra- 
diation has a 15 GW peak power and 25 pJ pulse energy 
at a wavelength of 0.0783 nm, for an undulator length of 
100 m. In this situation, diffraction, longitudinal space- 
charge and edge effects are found to be negligible while 
both emittance and energy spread have an impact on the 
FEL performance, which may explain the slightly lower 
performance than in the previous scenario. 

In summary, laser-plasma accelerators have promising 
perspectives in the realization of compact free-electron 
lasers, by combining them with conventional undulators. 
The scenarios X-FEL 1 and 2 highlight the possibility 
to produce strong FEL radiation from a laser-plasma ac- 
celerator. The most challenging condition is the rela- 
tive energy spread of the electron beam, which has to be 
decreased to attain the x-ray part of the spectrum. In 
addition, high-quality electron beams at the GeV level 
are required [1 GeV electron beams have already been 
generated (Lccmans et ai, 2006) but not with the same 
quality, i.e., energy spread and stability, as at the 100 
MeV level]. The notion of slice energy spread and slice 
omittance can facilitate reaching the FEL requirements 
for the highest-energy case. 

2. With an electromagnetic wave undulator 

In Sec. VI, the incoherent radiation of an electron 
beam submitted to a counterpropagating electromagnetic 
wave was presented. This wave acts as an undulator or 
wiggler whose strength parameter K equals the dimen- 
sionless vector potential amplitude an of the laser pulse 
and whose period equals Al/2, where Al is the wave- 
length of the counterpropagating laser pulse. Currently, 
only two laser systems can deliver such a pulse with an on 
the order of unity: Ti:Sa and CO2, with wavelengths of 
0.8 and 10 pm, respectively. From the resonance condi- 
tion (110), electron beams at the 10 MeV level radiate in 
the keV (nanometer) range for CO2 lasers and in the 10 
keV (angstrom) range for Ti:Sa lasers. At the 100 MeV 
level, the radiation is in the 100 keV range and 1 MeV 



range for, respectively, CO2 and Ti:Sa lasers, and at the 
GeV level the radiation is in the 10 and 100 MeV range 
for, respectively, CO2 and Ti:Sa lasers. The use of an 
electromagnetic wave undulator in the FEL process was 
proposed in the 1980s (Danly et ai, 1987; Gallardo et ai, 
1988; Gea-Banacloche et ai, 1987; Mima et ai, 1988) and 
reinvestigated recently (Bacci et ai, 2006, 2008; Maroli 
et ai, 2007; Petrillo et ai, 2008). It was also suggested 
to couple an electron beam coming from a laser-plasma 
accelerator and a relativistic laser pulse, in the conditions 
where FEL collective amplification takes place, and this 
system has been called an all-optical free-electron laser 
(AOFEL) (Petrillo et al, 2008). 

According to Pierce parameter scalings (115), for a 
given radiation wavelength A, using a smaller undulator 
period A„ and smaller 7 is unfavorable: the Pierce pa- 
rameter p decreases and the requirements on the electron 
beam parameters are more stringent for the electromag- 
netic undulator than for the conventional undulator (re- 
member that for constant A, p scales as 7 1 / 3 ). The most 
realistic scenarios for the electromagnetic undulator use 
high values of A, i.e. CO2 laser and electron beams at 
the lOMeV level. 

Petrillo et al. (2008) simulated the FEL collective am- 
plification taking place in the Thomson backscattering of 
a relativistic laser pulse by an electron beam, this latter 
generated by a laser-plasma accelerator. First, they use 
a 2.5D particle in cell simulation (3D in the fields and 
2D in the coordinates) to obtain electron beam parame- 
ters. They considered the case of an injection realized by 
a down ramp in the plasma density profile and showed 
that, after selecting a particular slice of the electron beam 
of longitudinal size 0.5 pm, it is possible to obtain a slice 
energy spread of 0.4%, a slice transverse emittance of 
0.3 7r.mm.mrad, and a current of 20 kA at electron ener- 
gies of about 30 MeV. Then, inserting the exact shape of 
the slice in the FEL code Genesis 1.3 [using the equiva- 
lence between conventional and electromagnetic undula- 
tors developed by Bacci et al. (2008)], they obtained for 
the SASE case the radiation output as a function of the 
undulator length: for a CO2 laser pulse with normalized 
vector potential an = 0.8, coherent radiation at wave- 
length A = 1.35 nm with subfemtosccond duration and 
200 MW peak power is generated for a saturation length 
of a few millimeters. These results open perspectives for 
an x-ray AOFEL using millimeter-scale electromagnetic 
wave undulators. However, their simulation does not take 
into account the short propagation from the plasma exit 
to the electromagnetic wave undulator, which could de- 
grade the beam parameters at such low energies and high 
currents. 

To conclude, the AOFEL has the main advantage of 
being compact and not needing a transport system. Be- 
cause of the short undulator period, high photon energies 
can be achieved with small electron energies. Although 
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it allows one to theoretically generate hard x rays and 
7 rays, the FEL amplification for such wavelength de- 
mands extremely stringent conditions on energy spread 
and emittance (the Pierce parameter is smaller for higher 
photon energies). Hence, it is more reasonable to use elec- 
tron energies at the 10 MeV level and CO2 laser pulses 
in order to succeed in operating a free-electron laser. Ex- 
perimentally, good quality electron beams (with percent 
energy spread) have been demonstrated at the 100 MeV 
level (Rcchatin et al, 2009a), but not yet at the 10 MeV 
level. It is generally thought that, because acceleration 
conserves absolute energy spread, smaller relative energy 
spread is easier to obtain at higher electron energies. 
Hence, achieving very small relative energy spread at the 
low-energy 10 MeV level seems to be difficult. From the 
perspective of an AOFEL, it is therefore necessary to ex- 
perimentally investigate the production of high-quality 
electron beams at rather small electron energies, when 
currently the tendency is toward higher and higher elec- 
tron energies. 

3. With a plasma undulator 

As shown in Sec. IV, the ion cavity, generated be- 
hind the laser pulse when the bubble or blowout regime 
is attained, acts as an undulator or wiggler whose pe- 
riod is A u = y/2r/\ p and whose strength parameter is 
K = rpk p \J~f /2. Assuming that we can prepare exper- 
imental conditions in which the electron beam keeps a 
constant energy inside the ion cavity (7 is constant and 
A. lt is too), then FEL amplification can in principle take 
place. For the case of a ion channel, the process has been 
studied (Whittum, 1992; Whittum et al, 1990) and is 
called an ion channel laser (ICL). Additional conditions 
have to be fulfilled: all electrons have to oscillate in the 
same plane (for the radiation to be polarized identically 
for all electrons) and with the same strength parameter 
(to radiate at the same wavelength) (Esarcy et al., 2002; 
Kostyukov et al., 2003). 

The main difference between ICL and FEL is that, in 
ICL, the strength parameter K depends on both the en- 
ergy and the transverse amplitude of motion r@, when 
FEL it is a constant defined uniquely by the undulator 
parameters. In ICL, different values of K lead to an ef- 
fective energy spread in ^ z = 7/yl + K 2 /2. Imposing 
that this energy spread is smaller than the Pierce param- 
eter p leads to the condition Arp/rp < p{2 + K 2 )/K 2 for 
the admissible relative spread in the amplitude of motion 
rp. Hence, the realization of ICL needs a configuration in 
which there is no acceleration to satisfy the condition of 
constant 7, where there is an off-axis injection with am- 
plitude rp and corresponding relative amplitude spread 
roughly below p. To our knowledge, a method to im- 
plement such a configuration is currently not known. In 
addition, conditions analogous to those presented in Sec. 



VII. A. 2 regarding the natural energy spread, emittance, 
and undulator length, also have to be fulfilled. 

VIII. CONCLUSION 

In this paper, mechanisms that can produce short- 
pulse x-ray or gamma-ray radiation using laser- 
accelerated electrons have been described. We have seen 
that all these sources are based on the radiation from a 
relativistically moving charge, and that even if the de- 
tails of the electron orbits are different for each source, 
the main features can be obtained using only five param- 
eters: the relativistic factor 7 of electrons, the number of 
electrons N e , the strength parameter K, the period X u , 
and the number of periods N of the undulator or wiggler. 
In addition, the incoherent radiation provided by these 
sources can eventually be amplified to coherent radia- 
tion via the FEL mechanism for sufficiently high-quality 
electron beams; the use of conventional undulators rep- 
resents the most promising route toward such goal. Here 
a summary of the features of these sources is given. 

Table II gives typical parameters and features as well 
as scaling laws for the incoherent sources based on the 
use of a plasma undulator (betatron radiation), a con- 
ventional undulator or an electromagnetic wave undula- 
tor (Thomson backscattering) . For the produced sources, 
interaction parameters accessible with currently available 
tens of TW-class lasers have been considered. The pho- 
ton spectra obtained with each mechanism show that an 
energy range from eV to MeV can be covered. The main 
difference between the different sources reviewed resides 
in their photon energy. Depending on the spectral range 
required for a desired application, Table II helps to select 
the most suitable type of source. 

In the XUV domain, the conventional undulator cou- 
pled with a laser-plasma accelerator has been discussed. 
The emitted radiation can be monochromatic if the elec- 
tron bunch is monoenergetic. Up to ~ 10 9 photons col- 
limated within a few milliradians could be produced in 
a single shot using 10 9 electrons (160 pC) in the hun- 
dreds of MeV energy range and an undulator of strength 
parameter K = 1 with 100 periods. The radiation 
wavelength can be tuned by varying the electron energy. 
Because the acceleration and the wiggling are dissoci- 
ated, a stable and/or tunable regime of acceleration can 
be used, resulting in a stable and/or tunable radiation 
source. From the perspective of GeV electron beams, 
this scheme will allow one to reach the keV range. How- 
ever, this source is more complex to realize because the 
electron bunch must be transported into the undulator. 
The transport can degrade the electron bunch duration 
(and so the x-ray pulse duration) and the transverse emit- 
tance, leading to a broader bandwidth in the radiation 
spectrum and a smaller brightness. Finally, this scheme 
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TABLE II The first subtable gives typical leatures for the different sources achievable with a 50 TW-class laser. The values 
represent the orders of magnitude of the parameters. The second subtable provides the scalings of the radiation features with 
the relevant parameters K, \ u , and 7 derived in Sec. II for sinusoidal trajectories. The last subtable presents the scalings 
of the relevant parameters with the practical parameters for each schemes derived from simple ideal models which have been 
described in each corresponding section. The practical parameters of each source are indicated in brackets. We refer to each 
section for the definition of the quantities used in this table. 
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Scalings of the radiation features 



Fundamental radiation energy hu> for undulators (K < 1) (2-y hc/\ u )/{l + K /2) 

Critical radiation energy huj c for wigglers (K ^> 1) ^K^ 2 hc/\ u 

Typical opening angle of the radiation 6 r for undulators (K < 1) 1/7 

Typical opening angle of the radiation Q r for wigglers [K > 1) K/y 
Number of emitted photons per electron and per period AT 7 , 

for undulators (K < 1) 1.53 x lO -2 ^ 2 

for wigglers (K > 1) 3.31 x Vd~ 2 K 



Scalings of the relevant parameters K X u 

Betatron (n e , rp, and 7) yj^j /2k p rp ^f2r)\ v 

Conventional undulator (Bo, A u , and 7) e_BoA„/(27rmc) \ u 

Thomson backscattering (ao, Al, and 7) ao Al/2 



deserves to be developed as it is on the path toward a 
free-electron laser based on laser-plasma accelerators. 

In the soft x-ray range, we discussed nonlinear Thom- 
son scattering. This radiation is produced by electrons 
oscillating in an intense laser field. While this source has 
a spectrum extending in the few hundreds of eV range 
using tens of TW-class lasers, it may become an inter- 
esting mechanism to generate keV x rays with PW-class 
lasers. 

In the range of a few keV, the betatron mechanism has 
been demonstrated to be an efficient and simple method 
to produce short x-ray pulses. This source can produce 
up to 10 9 photons per shot, collimated within less than 
50 mrad (FWHM) . On the one hand, controlling the elec- 
tron orbits within the cavity in order to ensure a better 
energy transfer from the electron to the radiation could 
allow one to extend the spectral range to the tens of keV 
range, while still working with 50 TW-class lasers. On 
the other hand, petawatt-class lasers and/or capillaries 
would be suitable to create larger ion cavities with lower 
plasma density and to accelerate electrons to the GeV 
range. This should allow one to produce x-ray beams 
whose divergence is in the milliradian range and photon 
energies in the tens and hundreds of keV range, with a 
higher number of photons. 

Above 10 keV and up to a few MeV, the mechanism 
used is Thomson backscattering. Here high-energy ra- 
diation is produced thanks to the short period of the 
electron motion. This scheme is similar to the case of a 
laser-plasma accelerator coupled with a conventional un- 



dulator, except that it benefits from the short period of 
the electron motion resulting in emitted photons of high 
energies and it does not need transport. The photon en- 
ergy can be chosen by tuning the electron energy. This 
scheme has been demonstrated recently and it is the most 
promising for generating tunable incoherent bright hard 
x-ray and 7-ray radiation. 

The different schemes reviewed here produce incoher- 
ent radiation whose energy is proportional to the number 
of electrons contributing to the emission. These sources 
can be compared, in some ways, to conventional syn- 
chrotron radiation. The main advantage provided by 
the laser-plasma approach is, of course, the reduced size 
(from hundreds of meters or kilometers to the laboratory 
room) and the reduced cost (from multihundrcd million 
dollars or multibillion dollars to multimillion dollars). 
Moreover, the pulse duration of laser-based sources is 
femtosecond when it is commonly picosecond for conven- 
tional synchrotrons even if, with the latest technology, 
it can be reduced down to femtosecond (and this has 
a certain cost). Another advantage is the natural syn- 
chronization of these sources with exciting laser pulses 
needed in pump-probe experiments, since both the pump 
and the probe come from the same laser system. Again, 
external femtosecond synchronization in synchrotron fa- 
cilities is feasible with recent methods, but it remains 
difficult. The femtosecond x-ray sources based on the 
relativistic laser-plasma interaction presented in this re- 
view are attractive mainly because of their low cost, and 
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this is why perspectives cannot be extended indefinitely 
to the use of higher intensity, since the cost will be- 
come comparable with conventional accelerators and syn- 
chrotron facilities. However, compared to synchrotrons, 
laser- and plasma-based sources present significant shot- 
to-shot fluctuations in terms of photon energy, number 
of photons, and divergence. They can therefore not cover 
a range of applications as large as synchrotrons. These 
sources remain at the research stage with limited appli- 
cations for the moment. Important challenging develop- 
ments are required to produced stable and robust x-ray 
beams appropriate for applications. In addition, the rep- 
etition rate is problematic: while synchrotron facilities 
work at MHz repetition rate allowing one to consider- 
ably increase the signal-to-noise ratio, laser-based sources 
work at 10 Hz on paper, but in practice they work in the 
single-shot regime in most 50 TW-class laser facilities to 
prevent rapid damage on optics such as in compressors. 
This means that theses sources are efficient for a single- 
shot experiment or experiments requiring not too much 
data accumulation. 

In Fig. 41, the x-ray sources are compared to other 
existing x-ray sources in terms of their peak brightness. 
This includes x-ray sources based on the laser-plasma 
interaction such as high harmonic generation from rela- 
tivistic laser and overdense plasma interaction (Dromcy 
et al, 2007, 2006; Tarasevitch et al, 2000; Teubner and 
Gibbon, 2009; Thaury et al, 2007), and K a radiation 
from laser-produced plasmas (Kieffer et al, 1993; Mur- 
nane et al., 1991; Rousse et at, 1994), as well as x-ray 
sources based on conventional accelerators: the undu- 
lator radiation source from the European Synchrotron 
Radiation Facility (ESRF), the femto-slicing-based un- 
dulator radiation source from the femto-laser station at 
the Paul Scherrer Institut (PSI), and the Thomson 7-ray 
source from the Lawrence Livermore National Labora- 
tory (LLNL) (Albert et al, 2010). Figure 41 shows that 
laser-plasma-based x-ray sources have competitive peak 
brightness, in particular, the betatron radiation and the 
all-optical Thomson source, because of their femtosecond 
duration and their micron source size. 

As free-electron lasers are the natural evolution of syn- 
chrotrons, FEL based on laser-plasma interaction could 
be the next step in the development of laser- and plasma- 
based x-ray sources. In the free-electron laser process, 
the emitted radiation is amplified and becomes coher- 
ent with the number of photons increased by orders of 
magnitude. The most promising scenario uses a con- 
ventional undulator coupled to a laser-plasma acceler- 
ator. Because stringent conditions apply to the electron 
beam parameters such as energy spread and transverse 
omittance, challenging developments are required in the 
laser-plasma accelerator domain, including ultracompact 
transport system. Thanks to the high brightness of FEL 
radiation, many experiments could be performed in single 
shot or with small data accumulation such that the low 
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FIG. 41 Peak brightness for different types of x-ray sources: 
high harmonic generation from relativistic laser and overdense 
plasma interaction (triangle and yellow line) [experiment con- 
ducted with the Vulcan Petawatt laser system at Ruther- 
ford Appleton Laboratory (RAL) by Dromey et al. (2006)], 
K a radiation from laser-produced plasmas (black star) [ex- 
periment conducted at the Laboratoire d'Optique Appliquee 
(LOA) by Rousse et al. (1994)], nonlinear Thomson scat- 
tering (NTS) from relativistic laser and underdense plasma 
interaction (diamond and cyan line) [experiment conducted 
at LOA by Ta Phuoc et al. (2003b)], conventional undula- 
tor radiation from laser-plasma accelerators (green square) 
[experiment conducted with the ATLAS laser system at the 
Max-Planck-Institut fur Quantenoptik (MPQ) by Fuchs et al. 
(2009)], betatron radiation from laser-plasma accelerators 
(pentagon and blue line) [experiment conducted with the Her- 
cules laser system at the University of Michigan (U-Mich) by 
Kneip et al. (2010)], Thomson backscattering (TB) radiation 
from laser-plasma accelerators (circle and red line) [experi- 
ment conducted at LOA by Ta Phuoc et al. (2012)], undu- 
lator radiation from conventional accelerators (green dotted 
line) [data from the European Synchrotron Radiation Facility 
(ESRF); see the ESRF website], femto-slicing-based undula- 
tor radiation from conventional accelerators (magenta dash- 
dotted line) [data from the femto-laser station at the Paul 
Scherrer Institut (PSI); see the PSI website and Beaud et al. 
(2007)] and Thomson backscattering (TB) from conventional 
accelerators (black dashed line) [experiment conducted at the 
Lawrence Livermore National Laboratory (LLNL) by Albert 
et al. (2010)]. The repetition rate and the pulse duration used 
in the peak brightness estimation are given in the legend, and 
permit one to convert peak brightness into average brightness. 
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disadvantage than for incoherent sources. 

In conclusion, this article gives a comprehensive review 
of the different mechanisms proposed for the production 
of femtosecond x-ray radiation from laser-plasma accel- 
erators. These sources have been described within the 
same formalism, considering different configurations of 
accelerator and undulator or wiggler and extracting the 
relevant parameters such as K, X u , and 7. Finally, the 
possibility to implement a free-electron laser using these 
configurations has been discussed and the electron beam 
quality required for such realization has been presented. 
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